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Scattering  from  Inclusions  in  Marine  Sediments: 
SAX04  Data/Model  Comparisons 


Anatoliy  N.  Ivakin 


Abstract — The  role  of  discrete  scatterers  in  marine  sediments  is  evaluated  based  on  acoustic  and  environmental  measurements  at  the 
shallow  water  sediment  acoustics  experiment,  SAX04,  and  shown  to  be  significant.  The  sediment  at  SAX04  site  was  complicated  and 
characterized  as  a  mostly  medium  sand  /  mud  mixture.  Analysis  of  the  sediment  samples  showed  presence  of  a  small  volume  portion  of 
larger  particles,  such  as  coarse  sand  fraction  and  shell  fragments,  which  were  considered  as  inclusions,  or  sparsely  distributed  discrete 
scatterers,  embedded  in  a  homogeneous  effective  fluid  with  parameters  corresponding  to  medium  sand  or  mud  sediment  frame.  This 
analysis  provided  also  the  size  distributions  for  both  types  of  inclusions,  coarse  sand  particles  and  shells,  which  were  used  as  inputs  to  a 
model  of  incoherent  discrete  scattering  in  the  sediment.  It  is  demonstrated  that  this  approach  in  general  is  capable  to  provide  a 
reasonable  explanation  of  both  frequency  and  angular  dependencies  of  the  SAX04  bottom  backscattering  strength,  and  discrete 
scattering  from  the  sediment  inclusions  can  be  a  significant  mechanism  of  the  SAX04  high  frequency  reverberation  in  a  wide  range  of 
grazing  angles  (15  to  50  degrees)  and  frequencies  (30  kHz  to  500  kHz). 

Index  Terms — Backscattering  strength,  medium  sand  and  mud  sediments,  shell  fragments,  coarse  sand  inclusions,  particle  size 
distribution. 


I.  Introduction 

Testing  various  models  and  mechanisms  of  high  frequency  seabed  scattering  was  one  of  primary  goals  of  recent  major  shallow 
water  experiments  in  sediment  acoustics,  SAX99  [Williams  et  al  2002]  and  SAX04  [Thorsos  et  al  2005].  The  scattering 
mechanisms  are  related  with  different  types  of  seabed  irregularities:  volume  heterogeneity,  continuous  (smooth  spatial 
fluctuations  of  the  sediment  acoustic  parameters)  and  discrete  (large  particles  in  the  sediment),  roughness  of  the  seabed 
interfaces,  and  volume -roughness  interactions.  Given  different  environmental  conditions,  frequencies  and  angles  of  scattering, 
relative  contributions  of  these  mechanisms  can  be  also  very  different.  For  example,  in  the  SAX99  situation,  sediment  rough 
surface  scattering  was  demonstrated  as  a  dominant  mechanism  at  frequencies  below  130  kHz  and  sub-critical  grazing  angles 
[Williams  et  al2002]. 

Environmental  conditions  at  SAX04  were  quite  different  from  those  at  SAX99  due  to  known  weather  events  preceding  and 
during  the  experiment,  and,  as  a  result,  the  SAX04  sediment  was  much  more  complicated  and  generally  characterized  as  a 
complex  mixture  of  mostly  medium  sand  and  mud  [Richardson  et  al  2005,  Williams  et  al  2008].  Acoustic  backscattering  from  the 
sediments  at  SAX04  was  also  very  different  from  SAX99.  For  example,  the  seabed  roughness  has  been  shown  to  be  insignificant 
contributor  to  the  sediment  backscatter  [Williams  et  al  2008].  It  could  be  expected  that  other  scattering  mechanisms  can  be 
important  in  these  complicated  conditions.  The  objective  of  this  paper  is  to  evaluate,  using  a  theoretical  modeling  approach,  the 
role  of  discrete  scattering  mechanism  for  the  SAX04  environment.  This  mechanism  of  scattering  can  be  due  to  various  discrete 
“inclusions”  in  the  sediment,  which  can  be  defined  as  relatively  large  particles  with  the  size  larger  than  the  mean  grain  size  of  the 
sediment,  such  as  rocks,  coarse  sand  and  gravel  particles,  shells  and  shell  fragments,  shelled  animals,  etc.  It  is  assumed  that 
acoustic  measurements  of  the  seabed  backscattering  strength  made  at  SAX04  [Williams  et  al  2008]  along  with  environmental 
characterization  of  the  sediment  at  this  site  and,  particularly,  analysis  of  the  sediment  inclusions  [Ivakin  2008]  can  provide  a 
necessary  ground  truth  for  such  an  evaluation.  If  this  role  is  proved  to  be  significant,  this  approach  can  be  used  to  provide  a 
physics-based  model  for  analysis  of  bottom  reverberation  at  SAX04  and  other  shallow  water  regions  with  similar  conditions. 

The  SAX04  acoustic  experiments  were  accompanied  by  extensive  measurements  of  the  sediment  acoustic  parameters,  the 
density,  sound  speed  and  attenuation,  and  others  [Richardson  et  al  2005].  These  parameters  are  necessary  inputs  for  description 
of  sound  propagation  in  the  sediment,  reflection  from  and  transmission  through  the  sediment  interfaces.  To  describe  the  process 
of  scattering  in  the  sediment,  another  set  of  sediment  parameters  is  required.  In  a  companion  paper  [Ivakin  2008],  we  present 
results  of  granulometric  analysis  of  the  SAX04  sediment  samples  for  particles  of  coarse  fractions  with  the  grain  size  (equivalent 
diameter)  greater  than  1  mm.  These  particles  comprise  only  a  small  volume  portion  in  the  SAX04  medium  sand  /  mud  sediment, 
and  can  be  considered  as  sparse  discrete  “inclusions”  in  the  sediment.  In  this  paper  (Section  II),  a  model  of  incoherent  discrete 
scattering  in  the  sediment  is  presented  and  discussed.  The  model  considers  inclusions  as  discrete  scatterers  embedded  in  a 
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sediment  half-space  with  given  effective  fluid  parameters,  density,  sound  speed  and  attenuation.  Acoustic  scattering  from  the 
sediment  inclusions  is  described  in  terms  of  the  individual  scattering  functions  of  discrete  targets  and  statistical  distributions  of 
their  parameters  (size,  shape,  orientation,  material,  etc).  Parameters  of  these  distributions  (or  at  least  possible  limits  of  their 
values)  can  be  obtained  from  an  analysis  of  the  sediment  samples.  In  this  paper  (Section  III),  the  results  of  analysis  of  the  SAX04 
sediment  samples  [Ivakin  2008]  are  used  to  provide  the  size  distributions  of  the  two  types  of  inclusions,  shell  fragments  and 
coarse  sand  particles,  as  the  input  parameters  to  the  “inclusion  scattering”  model.  In  Section  IV,  this  model  is  used  to  explain 
frequency-angular  dependencies  of  seabed  scattering  strength  measured  at  SAX04  experiment  over  a  wide  frequency  range  (30 
kHz  to  500  kHz)  at  grazing  angles  both  above  and  below  critical  (about  30  degrees  for  medium  sand).  In  Section  V,  some 
discrepancies  in  model/data  comparisons  and  assumptions  made  in  the  model  of  scattering  are  discussed  and  possible  future 
modifications  of  the  model  to  eliminate  these  discrepancies  and  assumptions  are  outlined.  Section  VI  gives  the  paper  summary 
and  conclusions. 


II.  Model  of  incoherent  scattering 

The  “inclusion  scattering”  model  is  based  on  a  theoretical  approach  described  in  [Ivakin  2000,  2004,  2005].  It  assumes,  in 
particular,  that  the  inclusions  are  sparsely  distributed  in  space  and  comprise  only  a  small  portion  of  the  scattering  volume.  The 
sparse  distribution  of  inclusions  allows  considering  their  positions  as  independent  and  mutually  uncorrelated  and, 
correspondingly,  the  incoherent  summation  of  their  scattered  intensities.  This  approach  is  rather  general  and  is  often  applied  to 
various  types  of  scatterers  in  the  ocean,  such  as  fish  in  water  column,  clouds  of  gas  bubbles  near  the  sea  surface,  naturally 
occurring  plumes  in  the  ocean,  and  others,  see,  e.g.,  [Medwin  and  Clay  1998,  and  Refs  therein].  A  similar  approach  has  been 
widely  used  in  acoustic  sensing  of  suspended  sediments  [Sheng  and  Hay  1988,  Thorne  et  al  1995,  Thorne  and  Backingham 
2004],  which  assumes  that  the  particles  are  small  and  light  enough  so  that  they  can  be  suspended  in  water  above  the  seabed.  This 
allows  particularly  estimating  the  particle  sizes  from  multi-frequency  acoustic  backscatter  at  very  high  frequencies  (about  1  MHz 
and  higher)  corresponding  to  suspended  sediment  grain  sizes. 

The  case  of  non-suspended  sediment  (or  large  particles  and  correspondingly  lower  frequencies)  is  more  complicated  in  the 
sense  that  such  sediment  is  generally  a  densely  packed  granular  medium,  and  a  complex  process  of  interactions  between  these 
large  particles  and  surrounding  smaller  grains,  presence  of  the  sediment  interfaces  and  related  effects  should  be  somehow  taken 
into  account.  It  seems  reasonable  however  to  assume  that  acoustic  behavior  of  a  large  particle  (or  “inclusion”)  surrounded  by  a 
big  number  of  small  grains  is  similar  to  the  case  where  this  inclusion  is  submerged  to  an  effective  continuous  medium  with 
acoustic  parameters  corresponding  to  the  surrounding  sediment  frame.  The  effective  parameters  of  this  frame  and, 
correspondingly,  sound  propagation  in  the  surrounding  medium,  where  grains  are  much  smaller  than  the  wavelength,  can  be 
defined  by  various  “effective  medium”  models,  e.g.,  the  “effective  density  fluid”  model  [Williams  2001],  The  “inclusion 
scattering”  model  presented  and  discussed  in  this  paper  represents,  in  a  way,  represents  a  hybrid  approach  to  discrete  scattering  in 
non-suspended  sediments  using  advantages  of  both  approaches,  an  “effective  medium”  approach,  in  part  of  description  of 
propagation  in  the  sediment,  and  the  suspended  particle  approach,  in  part  of  description  of  local  process  of  discrete  scattering  in 
the  sediment. 

An  essential  question  in  the  case  of  non-suspended  sediments  remains  when  it  is  hard  to  distinguish  the  inclusions  from 
surrounding  granular  sediment  frame.  In  this  case,  the  shape  of  the  sediment  grain  size  distribution  becomes  important  and  should 
be  considered  in  more  details.  Usually,  the  sediment,  as  a  granular  medium,  is  being  considered  as  having  a  rather  narrow  size 
distribution  of  particles  and  its  mean  grain  size  is  assumed  as  a  major  parameter  defining  the  sediment  geotechnical  and  acoustic 
properties.  For  example,  traditionally  it  defines  the  sediment  type  and  related  “word-based”  characteristics  (e.g.,  fine,  medium,  or 
coarse  sand).  This  however  has  much  less  sense  if  the  sediment  has  two  or  more  “modes”  with  significantly  different  grain  sizes, 
or  has  a  wide  range  of  sizes.  In  this  case,  standard  granulometric  description  of  the  grain  size  distribution,  even  rather  extended 
and  involving  other  parameters  in  addition  to  the  mean  grain  size,  such  as  sorting,  skewness,  kurtosis,  etc,  see  [Jackson  and 
Richardson  2007],  with  respect  to  effect  on  the  sediment  scattering  properties,  is  simply  not  relevant.  From  other  hand,  presence 
of  a  small  volume  portion  of  particles  with  the  size  much  larger  than  in  the  main  mode  can  contribute  significantly  to  scattering 
properties  of  the  sediment,  although  it  does  not  affect  the  mean  grain  size  of  the  sediment.  A  more  accurate  analysis  of  the 
sediment  grain  size  distributions  shows  that  they  normally  do  have  particles  in  a  wide  range  of  sizes  larger  than  the  mean  grain 
size.  Particles  of  this  “tail”  in  the  size  distribution  are  usually  sparse  in  space  and  satisfy  the  requirement  for  their  total  volume 
concentration  to  be  small.  Important  is  that  they  do  not  affect  the  value  of  the  sediment  mean  grain  size.  Therefore,  in  the 
scattering  model,  these  particles  can  be  considered  as  inclusions  submerged  in  an  effective  continuous  medium  with  the  effective 
acoustic  parameters  corresponding  to  this  mean  grain  size.  Important  is  that,  in  this  case,  the  processes  of  propagation  in  this 
volume  and  local  scattering  by  a  single  inclusion  can  be  treated  separately. 

The  incoherent  scattering  approach  for  inclusions  in  underwater  sediments  can  be  described  as  follows.  Consider  a  sediment 
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volume  with  a  number  of  discrete  scatterers,  or  "inclusions”.  The  intensity  I  j  of  the  field  scattered  by  the  j- th  inclusion  located 

in  this  volume  is  proportional  to  the  intensity  incident  from  water  on  the  sediment  surface,  I () .  Note  that  it  is  not  required  here  to 

specify  the  intensity  of  the  field  incident  on  an  individual  inclusion.  At  far  enough  distances  ro  from  the  sediment  surface,  taking 

into  account  the  spherical  spread  in  water,  the  scattered  intensity  can  be  written  as 


!j  =  (JJ/o  lr0‘ 


(1) 


where  <Jj  is  a  factor  that  can  be  defined  as  the  efficient  target  scattering  cross-section  for  this  inclusion.  Regardless  of 

complexity  of  the  sediment  structure,  if  inclusions  are  sparse  enough  and  distributed  randomly,  their  scattered  intensities  can  be 
summed  incoherently.  This  summation  determines  the  seabed  scattering  coefficient  (or  the  scattering  cross  section  per  unit  seabed 
area,  or,  if  expressed  in  dB,  the  seabed  scattering  strength) 


m 


s 


A  - /  J 


(2) 


where  A  is  the  horizontal  cross  section  of  the  ensonified  sediment  volume.  It  can  be  assumed  here  that  the  scattering  cross 
section  of  each  inclusion  essentially  depends  on  its  depth  and  size.  Then  summation  in  (2)  can  be  performed  using  a  number- 
size-depth  distribution  function 

A N(a,z)  (3) 

AAa  Az 


where  AN  is  defined  as  the  number  of  particles  having  the  size  and  depth  in  small  enough  intervals,  Cl  j  e  \a,  U  +  Aa]  and 
Zj  e  [z,  Z  +  Az] .  Then  summation  in  (2)  can  be  replaced  by  integration  and  the  seabed  scattering  coefficient  at  given  frequency 
and  directions  of  the  incident  and  scattered  waves  can  be  written  as  follows 


m 


(ks,ki)  =  II  y/N  {a,  z)  &{ks  ,kt,a,z)  da  dz 


(4) 


where  k  and  k  are  the  wave  vectors  of  the  scattered  and  incident  waves  in  water  near  the  seabed.  The  expression  (4)  is  rather 

general.  It  is  assumed  here  only  that  all  the  particles  with  the  same  size  and  depth  intervals  belong  to  the  same  class  and  therefore 
their  scattering  intensities  are  averaged  within  this  class,  i.e.,  not  only  over  these  intervals,  but  also  over  their  horizontal  positions 
within  the  ensonified  area  and  all  other  possible  parameters,  such  as  their  material  properties,  orientation  (for  non-spherical 
particles),  and  others.  Alternatively,  the  number  of  essential  parameters  of  inclusions  should  be  increased,  which  means  that,  e.g., 
the  size  parameter,  a  ,  should  be  understood  as  a  vector,  or  a  set  of  parameters,  and  the  size  distribution  (3)  and  integration  in  (4) 
should  be  extended  on  these  parameters  as  well. 

The  expression  (4)  can  be  used  to  introduce  a  generalized  volume  scattering  coefficient  of  the  medium,  m,  or  the  scattering 
cross  section  per  unit  sediment  volume  element  located  at  given  depth,  as  follows 

ms  =  |  m(z)dz 


m(z)  =  J  y  N  {a,  z)a{a,  z)  da 


(6) 


Using  the  size-depth  distribution  function,  all  practically  important  statistical  characteristics  of  the  sediment  inclusions  can  be 
found.  Then,  using  (4)-(6),  the  seabed  scattering  coefficient  and  the  sediment  volume  scattering  coefficient  can  be  expressed 
through  such  parameters  as  the  average  cross  section  and  average  volume  of  the  inclusions,  their  average  volume  concentration, 
the  surface  "coverage”  parameter,  and  others  (see  Appendix  I). 

It  is  important  to  note  that  the  individual  scattering  cross-section  <J  in  the  expressions  above  is  defined  rather  generally.  It  is 
different  from  one  in  water  or  in  the  sediment  unbounded  (free)  space,  particularly,  because  it  is  affected  not  only  by  the 
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scattering  process  itself,  but  also  by  interactions  with  the  sediment  interfaces.  For  example,  in  the  case  of  one  interface  (water- 
sediment),  the  interaction  processes  include  transmission  of  the  incident  and  scattered  waves  through  the  interface,  their 
propagation  and  attenuation  in  the  sediment,  and  possible  multiple  reflections  between  the  interface  and  the  target.  All  these 
processes  are  included  to  the  individual  scattering  function  as  it  is  defined  by  expression  (1),  which  can  be  used  for  direct 
measurements  of  this  function.  Such  measurements,  however,  might  be  difficult  for  not  large  targets  as  their  scattering  amplitude 
can  be  too  small  compared  to  the  level  of  background  scattering  from  surrounding  granular  medium  (this  background  can  be 
rather  strong  at  high  frequencies  [Ivakin  and  Sessarego  2007]).  Other  approaches  can  be  related  with  computer  simulations  and 
theoretical  modeling,  which,  using  reasonably  simplifying  assumptions,  can  provide  approximate  solutions  for  this  function. 

In  this  paper,  we  consider  the  case  of  particles  which  are  buried  in  the  sediment  and  assume  that  all  possible  multiple 
reflections  between  the  water-sediment  interface  and  particle  can  be  neglected.  In  this  case,  the  process  of  scattering  from  the 

particle  in  the  sediment  can  be  described  by  a  “free-field”  scattering  cross  section  <J0  ( ks ,  ,  a) ,  i.e.,  by  that  in  unbounded 

homogeneous  sediment  free  space.  The  transmission  through  the  water-sediment  interface  and  propagation  in  the  sediment,  at 
given  directions  of  incident  and  scattered  waves,  can  be  taken  into  account  by  a  bistatic  propagation  factor,  T (k  s ,  k{ ,  z) ,  so  that 

<y(ks ,ki,a,z)  =  T (ks , kt , z)  <j0 (ks ,kt,a,z ) 

Correspondingly,  Eq.  (6)  now  becomes  of  the  form 

m(z )  =  T(z)m0(z) 


(7) 


(8) 


where  m  is  the  “free  field”  volume  scattering  coefficient  of  the  sediment  defined  as  follows 


m 


,(z)  =  \vn  (a’z)cr0  ( a >  z)  da 


(9) 


Assume  now  that  particles  are  uniformly  distributed  through  the  layer  of  thickness  h  ,  so  that  ZQ  <  Z  <  Z i  =  ZQ  +  h  .  Then 
for  the  seabed  scattering  coefficient,  using  Eq.(5),  one  obtains 


ms  =m0Z(z0,zl ) 

zi 

Z(zj,z2)  =  jl(z)dz 


(10) 


Consider  the  case  of  the  sediment  inclusions  embedded  in  a  homogeneous  fluid  half-space.  The  propagation  factor  for  this  case  is 
described  in  details  in  Appendix  II.  It  is  shown  that  this  factor  can  be  presented  in  an  exponential  form,  see  Eq.  (35).  Then  for  the 
bistatic  seabed  scattering  coefficient,  using  Eq.(5),  one  obtains 

exp(-  kzg  )  -  cxp(-  KZ]  )  (11) 

Tomo 

K 


For  the  backward  direction,  assuming  also  that  the  top  boundary  of  the  scattering  layer  is  close  enough  to  the  water-sediment 
interface,  so  that  Kz  «  1 ,  Eq.  (11)  results  in  the  following  expression  used  for  calculations  later  in  this  paper 


ms  = 


\Wt\4 in 0  sin %t 


2ft2  0 


(l  -  exp(-  Ifih!  sin  %t )) 


(12) 


where  //  is  the  sediment-to-water  density  ratio,  W-  is  the  water-to-sediment  transmission  coefficient  of  the  incident  wave,  /)  is 
the  sound  attenuation  coefficient  in  the  sediment,  and  %t  is  the  grazing  angle  of  the  incident  wave  in  the  sediment  (defined  in 

Appendix  II).  This  expression  is  similar  to  given  previously  in  [Ivakin  2001]  and  generalizes  the  one  for  volume  scattering  from  a 
half-space 
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\Wj\4 m0  sin Xt 
m ,  = - - - 

*  2  ju2/3  (13) 

see,  e.g.,  [Ivakin  and  Lysanov  1981,  Ivakin  2000,  and  Jackson  and  Richardson  2007],  to  the  case  of  a  scattering  layer  of  a  finite 
thickness. 

To  define  the  sediment  volume  scattering  coefficient,  m  ,  in  according  with  Eq.  (9),  we  need  to  specify  the  form  of  the  size 

distribution  of  the  particles  and  their  individual  scattering  cross  section.  The  size  distribution  for  the  SAX04  sediment  discrete 
scatterers  was  measured  [Ivakin  2008]  and  will  be  discussed  in  the  following  section.  There  were  no  measurements  of  the 
individual  scattering  from  inclusions  in  the  sediment,  and  there  is  no  exact  solution  for  non-spherical  particles  of  irregular  shape, 
which  would  be  most  relevant  for  natural  sediment  particles.  In  this  paper,  we  will  use  an  approximate  “high-pass”  expression  for 
this  function,  given  in  Appendix  III.  Following  other  approaches  used  to  define  “high-pass”  expressions,  we  assume  that  after 
average  over  random  orientations  of  the  particle,  it  can  be  approximated  by  that  for  spherical  particle  with  the  same  equivalent 
radius  (the  radius  of  a  sphere  having  the  same  volume).  With  this  modification,  the  sediment  particle  size  distribution  as  a 
function  of  the  equivalent  radius  is  required  to  provide  inputs  to  the  inclusion  scattering  model  to  provide  comparison  with  the 
SAX04  acoustic  backscatter  data. 


III.  Some  results  of  analysis  of  sax04  sediment  samples 


In  a  companion  paper  [Ivakin  2008],  we  present  detailed  results  of  granulometric  analysis  of  thirteen  SAX04  sediment  samples 
for  particles  of  coarse  fractions  with  equivalent  diameter  greater  than  1  mm.  The  two  visually  different  classes  of  grains  in  these 
fractions,  coarse  sand  (mainly  quartz  particles)  and  shell  fragments,  were  segregated  from  each  other  and  analyzed  separately. 
The  size  distributions  of  particles  as  functions  of  the  equivalent  diameter  D,  D=2a,  were  obtained  and  the  results  are  shown  in 
Fig.  1  for  sand  and  shell  particles  in  all  the  sediment  samples.  An  important  result  of  this  analysis  is  that  the  size  distributions  of 
these  two  classes  of  particles  closely  follow  the  power  law  in  all  the  range  of  sizes  presented  in  these  samples,  1mm  <  D  <  4  mm 
for  coarse  sand  particles,  and  1mm  <  D  <  20  mm  for  shell  fragments.  It  is  seen  also  that  size  distributions  are  much  flatter  for 
shells  while  decreasing  steeply  as  the  size  increases  for  coarse  sand  particles,  and  this  results  in  that  shells  become  a  dominating 
coarse  fraction  of  the  sediment  for  particles  larger  than  certain  “transition”  size  found  here  to  be  about  2  mm. 


Fig.  1.  Number-size  distributions  for  the  SAX04  sediment  inclusions  obtained  from  thirteen  SAX04  sediment  samples. 
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These  results  are  used  here  in  this  paper  to  provide  a  set  of  input  parameters  for  the  seabed  scattering  model.  First,  however,  it 
is  important  to  notice  that  presented  here  size  distributions,  as  usual  in  practice,  result  from  an  average  over  the  sediment  sample 
thickness  and,  therefore,  are  related  with  actual  size-depth  distributions  through  an  integral  relationship 

VSN(a)  =  H~l^y/N(a,z)dz  (14) 


where  H  is  the  sediment  sample  thickness.  Here,  for  the  SAX04  samples,  H  —  6  cm.  Correspondingly,  there  is  a  resolution 
uncertainty  in  defining  the  local  volume  concentration  of  the  particles  within  the  sediment  sample.  Also,  one  needs  to  assume  if 
there  are  particles  beyond  the  sediment  sample.  If  particles  are  concentrated  actually  in  a  layer  of  thickness  h  ,  we  use  following 
assumptions  to  resolve  this  uncertainties 

y/N  =  (H / h)y/N,  h<H 

y/N=y/N’h>H  (i5) 

This  means,  in  particular,  that,  if  particles  are  concentrated  in  a  thinner  layer,  h  <  H  ,  then  the  volume  concentration  in  this  layer 
(see  Appendix  I)  is  larger  than  average  in  whole  sample  and  can  be  defined  as  follows 

Cv  =  ( H/h)Cv ,  Cv  =  J v{aWN{a)da  (16) 


For  the  “effective  surface”  concentration  of  the  particles,  if  they  are  concentrated  in  a  layer  thinner  than  sample,  there  is  no  such 
uncertainty  and  it  can  be  defined  as  follows 


Cs  =  //J  s(a)y/N(a)da 


(17) 


In  following  model/data  comparisons,  we  use  two  models  for  the  depth  distribution  of  particles.  First  one  corresponds  to  a  half¬ 
space  with  uniformly  distributed  particles  having  the  same  volume  concentration  as  the  average.  In  the  second  model,  the 

thickness  of  the  scattering  layer  is  taken  to  be  smaller  than  for  the  sample,  h  =  1.5  cm,  and  therefore,  Cv  =  4  Cv . 

To  provide  inputs  for  the  theoretical  model  of  scattering,  we  use  power  law  approximations  for  the  size  distribution  function  of 
the  form 


yr  N(a)  «  Ba0  4  (a  /  a0)  7 ,  al<a<a2, 

where  ao  being  introduced  as  an  intercept  scale.  This  also  makes  a  power  law  amplitude,  Bf) ,  dimensionless.  Here  we  use  a 
scale  ao  =  0.5  mm.  For  shells  and  coarse  sand  particles,  Eq.(18)  is  used  with  two  different  sets  of  parameters.  In  Fig.l,  they  are 
represented  by  two  dashed  lines.  One  line  corresponds  to  coarse  sand  particles  and  was  calculated  with  the  following  values  of 
parameters:  Ba  =  0.018,  y  =  8,  =  0.4  mm,  a2  =  2  mm.  Using  Eqs.  (16),  (17),  the  particle  volume  concentration 

averaged  over  the  sample  volume  and  the  effective  surface  concentration  were  estimated  as  well  and  obtained  as  follows 
Cv  =  0.0058,  Cs  =  0.52  .  Corresponding  values  for  another  dashed  line  provide  inputs  to  the  model  of  scattering  for 
shells:  B  =  0.00045,  y  =  4,  fifj  =  0.5  mm,  a7  =  45  mm.  These  values  give  estimates  for  concentrations  of  shells: 

Cv  =  0.00 1 1,  C s  =  0.02 1 ,  which  are  smaller  than  for  coarse  sand  particles.  In  both  cases,  therefore,  including  the  4-times 
increase  of  volume  concentration  in  the  thinner  layer  (1.5  cm)  model,  we  have  Cv  «  1 ,  and  the  requirement  for  scatterers  to  be 
sparse  enough  is  not  violated.  Also,  in  both  cases,  we  have  C s  <  1  that  allows  considering  inclusions  as  a  thin  layer  of  particles. 

Based  on  these  estimates,  we  can  consider  shells  and  coarse  sand  particles  as  two  separate  modes  in  the  sediment  grain  size 
distribution,  or  as  sparse  inclusions  of  these  two  types,  suspended  in  a  homogeneous  and  uniform  substrate  of  densely  packed 
smaller  particles  (medium  sand  or  mud). 


7 


joe  858,  A.  Ivakin,  Scattering  from  Inclusions 

Submitted  by  author  to  IEEE  Journal  of  Ocean  Engineering  14  Jul  2008 


IV.  SAX04  DATA/MODEL  COMPARISONS 

The  obtained  sediment  particle  size  distributions  were  used  as  inputs  to  the  “inclusion  scattering”  model  to  provide  comparison 
with  the  frequency-angular  dependencies  of  acoustic  backscattering  from  the  seabed  observed  at  the  SAX04  site  in  wide  ranges  of 
frequency  (30  to  500  kHz)  and  grazing  angles  (15  to  50  degrees)  by  [Williams  et  al  2008].  The  first  and  main  question  to  answer 
using  this  data/model  comparisons  in  this  paper  is  if  the  level  of  scattering  corresponding  to  SAX04  sediment  inclusions  is 
comparable  with  the  observed  level  at  SAX04  and,  therefore,  if  the  inclusions  can  play  a  noticeable  role  in  the  SAX04  seabed 
scattering.  Next  question  is  if  a  simple  model  of  scattering  described  in  previous  sections  and  based  on  a  very  simplified 
presentation  of  the  complicated  SAX04  environment  is  able  to  provide  general  features  of  the  frequency-angular  dependences. 
Another  set  of  questions  for  discussion  is  which  features  of  the  observed  dependencies  do  not  appear  in  the  model  simulations, 
what  can  be  the  reason  of  the  discrepancies,  and  which  features  of  the  SAX04  environment  should  be  taken  into  account  in  future 
work  to  eliminate  such  discrepancies. 

Examples  of  typical  angular  dependencies  of  the  seabed  backscattering  strength,  10  log(/Ms )  ,  are  given  in  Fig.  2a  for  two 

frequencies,  70  and  400  kHz.  It  is  seen  that  they  are  quite  different.  The  angular  dependence  at  the  lower  frequency  shows  a 
distinct  “critical  angle”  effect  for  volume  scattering,  see,  e.g.,  [Ivakin  1981,  2001,  Jackson  2007],  with  the  critical  angle  being 
around  30°  which  corresponds  to  a  medium  sand  sediment.  At  the  higher  frequency,  this  effect  is  not  pronounced  or  at  least 
shifted  to  angles  smaller  than  20°  ,  where,  at  these  frequencies,  there  were  no  measurements.  Examples  of  frequency 
dependencies  are  shown  in  Fig.  2b  for  two  fixed  grazing  angles,  %w  =  35°  (above  critical  grazing  angle)  and  25°  (below 

critical).  It  is  seen  that,  first,  that  they  are  significantly  different  at  medium  frequencies,  70  to  130  kHz,  and,  secondly,  for  both 
angles,  they  have  different  slopes  at  low  and  high  frequencies.  A  general  idea  of  analysis  in  this  paper  is  that  such  frequency 
dependencies  of  the  scattering  strength  can  be  due  to  two  different  types  (or  modes)  of  inclusions  in  the  SAX04  sediment,  shells 
and  coarse  sand  particles. 

In  the  previous  section  of  this  paper,  it  was  demonstrated  that  the  size  distributions  of  the  two  modes  of  inclusions  have 
different  slopes  and  this  results  in  the  shell  mode  being  dominating  at  large  sizes  (greater  than  2  mm)  while  the  coarse  sand  mode 
dominating  at  smaller  sizes.  One  can  show  (see  Appendix  IV)  that  the  power  law  size  distribution  of  inclusions  results  in  the 
power  law  frequency  dependence  of  the  scattering  strength.  This  is  illustrated  in  Fig. 3,  where  the  normalized  volume  scattering 
strength  (described  in  details  in  Appendix  IV)  is  shown  as  a  function  of  dimensionless  frequency,  kao  ,  for  different  values  of  the 

power  exponents  y ,  at  fixed  cut  parameters,  ao  /  =  2,  ao  /  a2  =  0.02  .  The  following  analysis  shows  that  that  the 

difference  in  frequency  dependencies  of  the  scattering  strength  at  low  and  high  frequency  ranges  can  be  due  to  two  different  types 
(or  modes)  of  inclusions  in  the  SAX04  sediment,  shells  dominating  at  lower  frequencies,  and  coarse  sand  particles  dominating  at 
high  frequency  range.  This  explanation  holds  in  both  ranges  of  grazing  angles,  above  and  below  30°  . 

First,  calculations  of  the  backscattering  strength  were  made  using  Eq.  (13)  for  scatterers  uniformly  distributed  in  a  fluid 
sediment  half-space.  The  density,  sound  speed  and  attenuation  in  water  above  the  sediment  were  taken  as  follows:  pw  =1.0 

g/cm3 ,  Cw  =1.5  km/s,  8  =  0 .  Correspondent  parameters  for  the  inclusions  in  the  sediment  were  taken  to  be:  pa  =  2.7 

g/cm3 ,  Ca  =5.7  km/s,  8  =  0.  Values  of  effective  acoustic  parameters  of  the  sediment  half-space  in  the  model  were  taken  to 

be:  p  =  2.0  g/cm3 ,  C  =  1 .75  km/s,  8  =  0.006  ,  which  is  within  reasonable  uncertainty  of  environmental  measurements  at 
the  SAX04  site.  In  Figs  4a  and  4b,  frequency  and  angular  dependencies  for  the  backscattering  strength  for  the  two  modes  of 
particles,  shells  and  coarse  sand  (as  inclusions  in  medium  sand  frame  half-space),  are  presented  by  two  solid  lines 
correspondingly.  The  size  distributions  are  taken  for  both  cases  as  shown  in  Fig.l  by  two  dashed  lines  for  shells  and  coarse  sand 
particles  and  as  specified  in  previous  Section.  The  main  result  of  this  first  data/model  comparison  is  that  the  level  of  scattering  by 
the  SAX04  sediment  inclusions  is  comparable  with  the  observed  level  and,  therefore,  the  inclusions  indeed  can  play  a  noticeable 
role  in  the  SAX04  seabed  scattering.  Moreover,  this  simple  model  provides  a  good  explanation  of  both  frequency  and  angular 
dependencies  of  the  SAX04  backscatter  at  grazing  angles  starting  from  the  critical  angle  (about  30° )  and  above. 
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SAXD4  Backscatter  data 


Fig.  2a. 


SAX04  backscattering  at  various  grazing  angles 


Fig.  2.  Angular  and  frequency  dependences  of  the  bottom  backscattering  strength  measured  at  SAX04  by  [Williams  et  al  2008]. 


9 


joe  858,  A.  Ivakin,  Scattering  from  Inclusions 

Submitted  by  author  to  IEEE  Journal  of  Ocean  Engineering  14  Jul  2008 


ka 

O 


Fig.  3.  The  dependence  of  the  normalized  volume  scattering  coefficient  M ()  as  a  function  of  the  dimensionless  frequency  kao 
for  different  values  of  the  power  exponents  y ,  see  Appendix  IV. 


From  Fig.  4b,  it  is  easy  to  see  however  that  there  is  a  significant  discrepancy  of  the  observed  and  calculated  angular 
dependencies  at  grazing  angles  below  critical  for  both  low  and  high  frequencies.  To  understand  the  reason  for  this  discrepancy, 
one  should  take  into  account  that  actually  the  SAX04  sediment  was  not  purely  medium  sand,  but  rather  a  complex  medium  sand  / 
mud  mixture,  with  a  thin  (up  to  a  few  cm  thick)  mud  layer  covering  the  sand  basement.  A  consistent  consideration  of  the  effects 
of  the  sediment  layering  is  beyond  the  scope  of  this  simple  model  and  will  be  a  subject  of  future  work.  However,  some  estimates 
can  be  done  even  without  consideration  of  the  sediment  layering.  At  this  stage,  it  can  be  enough  to  consider  the  stratification  of 
scatterers  only,  that  is  to  take  into  account  a  finite  thickness  of  the  scattering  layer  within  a  mud  environment.  Physics  of  such 
simplification  is  that  in  this  case  all  possible  reflections  from  the  mud/sand  interface  are  neglected.  This,  however,  still  can 
provide  an  approximate  average  level  of  scattering,  but  will  neglect  all  possible  interference  patterns  in  frequency  angular 
dependencies  of  the  scattering  strength  usually  appearing  in  the  case  of  stratified  sediments  [Ivakin  1986  and  1998,  Mourad  and 
Jackson  1993]. 

For  such  rough  estimates,  calculations  were  made  for  the  same  model  of  “inclusion  scattering”,  but  using  Eq.  (12)  for  a  thin 
scattering  layer,  instead  ofEq.  (13)  for  the  scattering  half-space,  and  using  effective  acoustic  parameters  of  mud  sediments, 
instead  of  those  for  medium  sand.  The  scattering  layer  thickness  was  taken  to  be  1.5  cm,  and  parameters  of  the  mud  sediment 

were  taken  to  be  as  follows:  p  —  1.8  g/cm3 ,  C  =  1.61  km/s,  8  =  0.006  .  In  Figs  5a  and  5b,  frequency  and  angular 
dependencies  for  the  backscattering  strength  are  presented  by  two  dashed  lines  for  two  modes  of  particles,  shells  and  coarse  sand, 
correspondingly,  and  the  size  distributions  were  taken  again  for  modes  of  inclusions  as  shown  in  Fig.  1  by  two  dashed  lines. 

Figure  5a  shows  that  estimates  of  the  scattering  level  for  the  SAX04  sediment  inclusions  are  fitting  well  with  the  observed  at 
grazing  angles  below  30°  .  Figure  5b  also  shows  a  reasonable  agreement,  with  the  most  significant  discrepancy  at  low 
frequencies,  30  to  50  kHz.  Generally,  results  shown  in  Fig.5  demonstrate  that  inclusions  in  the  top  sediment  layer  of  mud  can 
indeed  be  considered  as  a  possible  mechanism  responsible  for  the  enhanced  level  of  scattering  observed  at  the  grazing  angles 
below  30°  in  both  low  and  high  frequency  range.  This  good  fit  in  data/model  comparison  however  is  paid  by  an  appearing 
inconsistency  in  the  sediment  model,  where  the  mud  inputs  have  to  be  used  at  small  grazing  angles,  below  30°  ,  while  medium 
sand  inputs  are  used  for  the  grazing  angles  above  30°  .  This  however  gives  a  hint  for  a  reasonable  modification  of  the  model  to 
eliminate  such  inconsistency.  It  seems  that  a  systematic  consideration  of  the  sediment  stratification  can  be  a  natural  next  step  in 
improving  the  model.  This,  along  with  other  ways  to  improve  the  model  of  scattering,  will  be  discussed  in  the  next  Section. 
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SAXD4  backscattering  at  grazing  angle  35° 


Fig.  4a. 


Fig.  4b. 


Fig.  4.  SAX04  data/model  comparisons  for  frequency  and  angular  dependencies  of  the  bottom  backscattering  strength.  The  solid 
lines  correspond  to  scattering  from  inclusions  with  the  power  law  grain  size  distribution  (see  dashed  lines  in  Fig.  1).  It  is  assumed 
that  inclusions  are  embedded  to  an  effective  fluid  half-space  with  parameters  corresponding  to  medium  sand. 
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Fig.  5a. 


SAXQ4  backscattering  at  grazing  angle  25° 
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Fig.  5b. 


Fig.  5.  SAX04  data/model  comparisons  for  angular  and  frequency  dependencies  of  the  bottom  backscattering  strength.  The 
dashed  lines  correspond  to  the  “inclusion  scattering”  model  with  the  power  law  grain  size  distribution  as  was  used  for 
calculations  in  Fig. 4  (shown  here  by  solid  lines  as  well  for  comparison).  The  difference  is  due  to  assumption  that  here  inclusions 
are  located  in  a  thin  top  layer  of  the  sediment  and  the  sediment  parameters  correspond  to  mud. 
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V.  Discussion 

The  SAX04  data/model  comparisons  were  made  for  the  model  of  scattering  using  certain  assumptions  and  simplifications, 
which  can  be  responsible  for  some  quantitative  discrepancies  and  an  inconsistency  noticed  in  previous  Section.  In  particular,  the 
following  assumptions  can  be  important: 

Inclusions  were  considered  as  fully  buried  in  the  sediment  assuming  that  contribution  of  particles  located  near  the 
interface  (including  those  partially  buried)  is  insignificant; 

Individual  scattering  function  was  taken  of  a  simplified  form  with  low-  and  high-frequency  limits  as  for  equivalent  fluid 
spheres  assuming  that  for  a  non-spherical  particle,  upon  average  over  their  random  orientation,  it  can  be  approximated 
by  that  for  spherical  one  with  the  same  equivalent  radius; 

An  incoherent  summation  of  intensities  scattered  by  different  inclusions  was  assumed,  neglecting  possible  effects  of 
their  correlations; 

Possible  multiple  scattering  effects  were  neglected; 

The  sediment  surface  was  assumed  flat,  so  that  possible  volume -roughness  interactions  were  ignored; 

The  sediment  was  considered  as  homogeneous  fluid  substrate,  i.e.,  its  stratification  was  not  taken  into  account. 

In  order  to  eliminate  these  assumptions,  corresponding  modifications  are  required.  For  example,  elimination  of  the  first 
assumption  is  necessary  to  take  into  account  scatterers  located  on  the  bottom  surface  or  not  buried  completely  but  rather  exposed 
in  the  water.  In  this  case,  analytical  solutions  are  complicated,  but  can  be  analyzed  numerically  [Lim  et  al  1993  and  Lim  1998, 
Fawcett  1996,  Fawcett  and  Lim  2003,  Zampolli  et  al  2008].  There  is  also  a  simpler  semi-phenomenological  approach,  using  ray 
considerations,  that  can  be  useful  at  high  enough  frequencies.  In  this  case,  the  incident  and  scattered  waves  can  be  considered 
taking  into  account  possible  reflections  on  the  water-sediment  interface  and  corresponding  bistatic  channels  of  scattering 
[Williams  et  al  2001].  Although  in  the  SAX04  conditions,  exposed  shells  and  other  discrete  scatterers  on  the  sediment  surface 
were  not  observed,  such  observations  are  common  in  many  other  shallow  water  regions.  This  can  be  especially  important  at 
higher  frequencies  and  sub-critical  grazing  angles  because  of  smaller  sound  penetration  of  the  sediment  and  respectively 
decreased  role  of  buried  scatterers.  Another  correction  should  be  made  for  scatterers  buried  near  the  water-sediment  interface,  so 
that  their  depth  is  smaller  or  comparable  with  their  size.  In  this  case,  their  individual  scattering  functions  can  be  affected  by 
possible  multiple  reflections  between  the  interface  and  the  particle.  However,  this  effect  can  be  neglected  if  the  contrast  of 
parameters  at  the  water-sediment  interface  is  small  enough  so  that  reflections  are  weak.  Note  also  that  these  reflections  can  be 
substantially  weakened  by  presence  of  a  thin  transition  layer  known  to  be  usual  near  the  water-sediment  interface. 

Elimination  of  the  second  assumption  is  particularly  important  for  better  description  of  the  angular  dependence  of  scattering 
which  is  known  to  be  sensitive  to  the  shape  of  the  scatterers  (the  so-called  aspect-effect)  that  is  commonly  used  in  remote  sensing 
(see,  e.g.,  [Ivakin  1981,  1985,  2001]).  Third  and  forth  assumptions  can  become  too  strong  at  higher  frequencies  (above  500  kHz), 
where  the  size  of  the  most  of  the  sediment  particles  in  the  medium  sand  sediment  becomes  comparable  with  the  wavelength  and, 
therefore,  “inclusions”  become  not  sparse.  In  this  case  corrections  should  be  considered  to  take  into  account  correlations  of  the 
particle  positions,  possible  “collective”  effects  and  effects  of  multiple  scatter,  see  discussions  in  [Bunchuk  and  Ivakin  1989, 

Lyons  2005,  Isimaru  1997].  Elimination  of  the  fifth  assumption  would  make  it  possible  to  give  a  consistent  consideration  of 
scattering  at  sub-critical  grazing  angles,  taking  into  account,  e.g.,  possible  effects  of  the  volume  scattering  enhancement  due  to 
better  sound  penetration  into  the  sediment  [Ivakin  1981,  2004,  2005]. 

The  sixth  assumption  neglects  possible  strong  effects  of  the  sediment  stratification  or  layering,  gradient  effects  in  the  top 
transition  layer,  related  effects  of  focusing,  appearing  new  channels  of  scattering  and  their  interference,  etc.  According  to  a 
results  of  the  data/model  comparisons  presented  in  previous  Section  of  this  paper,  the  sediment  stratification  seems  to  be  a  natural 
candidate  to  focus  in  future  work  to  better  understand  results  of  the  SAX04  acoustic  measurements.  This  corresponds  also  to  real 
situation  of  the  complex  SAX04  sand/mud  environment,  appearing  usually  as  a  thin  but  distinctive  mud  layer  covering  a  sand 
basement.  The  thickness  of  the  mud  layer  varied  from  a  few  mm  (forming  a  thin  transition  layer  between  water  and  sand)  to  about 
5  cm  (filling  depressions  in  the  sand  rough  surface  [Richardson  2005,  Thorsos  2005]).  One  of  possible  effects  of  the  transition 
layer  near  the  sediment  surface  is  smoothing  of  an  impedance  contrast  at  the  interface  which  can  result  in  an  additional 
enhancement  of  acoustic  penetration  and,  therefore,  enhances  volume  scattering.  It  is  known  that  the  sediment  stratification  can 
make  a  significant  impact  on  the  seabed  scattering  properties,  see  [Ivakin  1986  and  1998,  Mourad  and  Jackson  1993],  where  it 
was  shown,  for  example,  that  at  high  frequencies  and  small  grazing  angles,  because  of  the  smaller  sound  penetration  to  the 
sediment,  only  its  near  surface  properties  are  important  (and  the  SAX04  sediment  therefore  can  be  considered  as  having 
parameters  of  mud).  At  low  enough  frequencies  and  larger  grazing  angles,  the  thin  top  layer  parameters  become  unimportant  (and 
the  SAX04  sediment  can  be  treated  as  a  medium  sand  half-space).  In  this  paper,  we  consider  only  these  two  limit  cases,  leaving  a 
more  consistent  consideration  of  the  sediment  stratification  effects  to  a  future  work. 
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VI.  Summary  and  conclusions 

In  this  paper,  an  “inclusion  scattering'’  model  is  presented  and  used  for  analysis  of  backscattering  data  obtained  during  the 
SAX04  sediment  acoustics  experiment.  The  model  considers  discrete  scatterers  larger  than  mean  grain  size  as  inclusions 
embedded  in  a  sediment  half-space  with  given  effective  fluid  parameters,  density,  sound  speed  and  attenuation.  Acoustic 
scattering  from  the  sediment  inclusions  is  described  in  terms  of  the  individual  scattering  functions  of  discrete  targets.  It  was 
assumed  that  the  scattering  function  for  non-spherical  particles  upon  average  over  their  random  orientation  can  be  approximated 
by  that  for  spherical  one  with  the  same  equivalent  radius  (same  volume).  With  this  modification,  the  sediment  particle  size 
distributions  as  function  of  the  equivalent  radius  are  required  to  provide  input  parameters  to  this  model.  They  were  obtained  from 
granulometric  measurements  for  coarse  fractions  of  the  SAX04  sediment  particles  with  the  grain  size  larger  than  1  mm.  These 
particles  comprise  only  a  small  volume  portion  of  the  sediment  and  therefore  do  not  affect  the  mean  grain  size.  However,  they 
can  contribute  significantly  to  the  seabed  acoustical  properties  affected  by  processes  of  discrete  scattering  in  the  sediment. 

Generally,  the  results  of  data/model  comparison  demonstrate  that  scattering  from  the  sediment  inclusions,  such  as  shell 
fragments  and  coarse  sand  particles,  can  be  a  significant  or  even  a  dominant  mechanism  of  the  SAX04  seabed  scattering  in  a  wide 
frequency  range,  30  kHz  to  500  kHz,  and  in  a  wide  range  of  grazing  angles,  15  to  50  degrees.  A  possible  impact  on  improving 
existing  algorithms  for  geoacoustic  characterization  of  the  sediments  is  that  this  work  has  demonstrated  that  the  sediment  grain 
size  distribution  “tails”,  corresponding  to  particles  with  sizes  greater  than  the  mean  grain  size,  or  the  sediment  inclusions,  are 
critical  characteristics  and  important  factors  affecting  seabed  scattering  properties.  This  can  be  used,  for  example,  to  upgrade 
existing  grain-size  based  algorithms  for  predicting  bottom  reverberation  [Jackson  and  Richardson  2007,  Briggs  et  al  1989  and 
2002], 

The  analysis  has  shown  however  that,  although  the  current  version  of  the  model  is  able  to  provide  reasonable  estimates  for 
interpretation  of  SAX04  scattering  data,  some  quantitative  discrepancy  in  model/data  comparisons  appear.  This  can  be  due  to 
using  certain  assumptions  made  in  the  model  of  scattering.  It  is  clear  now  that  the  model  needs  some  additional  modifications  to 
improve  understanding  of  the  role  of  the  discrete  scattering  mechanism  in  the  sound-bottom  interactions,  and,  therefore,  enhance 
its  capabilities  for  potential  acoustic  data  predictions  and  inversions.  For  example,  taking  into  account  the  sediment  stratification 
might  be  an  important  step  in  improving  the  model.  This  and  other  modifications  can  be  interesting  subjects  for  future  work. 


Acknowledgment 

This  research  was  conducted  in  collaboration  with  colleagues  at  the  APL-UW  and  NRL  and  their  contribution  is  gratefully 
acknowledged.  Dr.  M.  Richardson  of  NRL  arranged  the  sediment  samples  collection  by  APL  and  NRL  divers.  Dr.  K.  Williams  of 
APL  provided  results  of  his  acoustic  measurements  at  SAX04  used  for  data/model  comparisons  in  this  paper.  This  work  was 
supported  by  the  U.S.  Office  of  Navel  Research. 


APPENDIX  I 

STATISTICAL  CHARACTERISTICS 

Using  the  size-depth  distribution  function,  all  practically  important  statistical  characteristics  of  the  sediment  inclusions  can  be 
found.  For  example,  from  the  definition  of  this  function  (3),  the  total  number  of  inclusions  in  the  given  sediment  volume  (or  in  the 
sediment  sample)  is 

N  =  A  JJ  y/N  (a,  z)  da  dz 

Normally,  the  scattering  volume  depends  on  the  directivity,  frequency  and  duration  of  the  signal.  The  number  of  scatterers  in  this 
volume  (or  in  the  “resolution  cell”)  plays  a  critical  role  in  statistics  of  scattering  [Abraham  2003,  Abraham  and  Lyons  2004].  In 
particular,  it  is  required  for  definition  of  the  probability  density  function 

(PN(a,z)  =(A/ N)y/  N(a,z)  (20) 

using  which  the  average  of  an  arbitrary  function  appearing  in  the  problem,  <  U  > ,  over  the  particles  size  and  depth  can  be 
defined  as  follows 
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<u  >=  ju(a,z)<pN(a,z)da dz 


(21) 


Consider  a  dimensionless  volume  concentration  of  particles  at  a  given  depth 

Cv  (z)  =  J  v(a,  z)y/N  ( a ,  z)  da 


Here  v{a,  z)  is  the  volume  of  the  particle  having  the  size  a  and  located  at  the  depth  z  .  This  dependence  generally  can  be 
different  for  particles  with  different  shape.  To  establish  such  a  relationship,  one  can  introduce  the  equivalent  radius  of  the 
particle,  i.e.,  the  radius  of  the  sphere  with  the  same  volume.  The  equation  (22)  is  important  for  controlling  the  requirement  of  not 
interpenetrating  of  particles,  for  which  at  any  depth  an  inequality  Cv(z)  <  1  should  hold.  Also,  comparing  (22)  with  (6),  the 
volume  scattering  coefficient  in  an  arbitrarily  thin  layer  of  sediments  can  be  defined  as  follows 

m  =  Cv  <  <7  >  /  <V  >  (23) 


For  the  average  volume  concentration  of  the  inclusions  in  the  sediment  sample,  assuming  that  the  sample  has  a  vertical  size  H , 
one  obtains 


c„ 


=  Cv(z)dz  = 


N  <  v  > 
AH 


(24) 


If  the  sediment  sample  is  thin  enough,  it  can  be  more  convenient  to  introduce  a  dimensionless  “coverage”  parameter,  or  an 
“effective  surface”  concentration  of  the  particles  in  this  sample 


y/N  {a,  z )  s(a,  z)  da  dz 


N  <s> 
A 


(25) 


Here  s(a,  z)  is  an  “effective”  cross-section  of  the  particle,  which  defines  a  geometric  projection  of  the  particle  onto  a  certain 
surface  in  the  sediment.  It  can  be  a  sediment  interface,  or  a  certain  plane  defined  in  the  problem,  e.g.,  oriented  normally  with 
respect  to  direction  of  sound  scattering.  This  dependence,  however,  can  be  different  for  non-spherical  particles  with  different 
orientation.  A  way  to  specify  such  a  relationship  can  be  averaging  over  possible  orientations  and  defining  an  equivalent  radius  of 
the  average  projection.  The  equation  (25)  can  be  also  important  in  the  case  where  a  monolayer  of  particles  is  considered.  In  this 
case,  an  inequality  Cs  <  1  should  hold. 

Using  (19)-(21),  it  is  easy  to  see  that  the  scattering  coefficient  in  (4)  can  be  obtained  from  the  equation 

ms  =CN(a)  (26) 


where  the  factor  CN 


can  be  defined  by  any  of  the  following  expressions 


N  _CS  _  CVH 

a  =  A=  <v> 


(27) 


and  the  choice  of  a  particular  expression  depends  on  convenience  for  a  specific  scattering  problem.  Consider,  e.g.,  the  case  where 
scatterers  are  distributed  within  a  layer  in  a  free-space.  In  a  particular  case  of  spherical  particles  with  the  radius  a  large  in 
comparison  with  the  wavelength,  the  individual  scattering  function  can  be  approximated  by  its  geometric  acoustics  expression 

( 7  =  a~R~  /  4  ,  where  R  is  the  reflection  coefficient  for  the  material  of  the  sphere  (see  Appendix  III).  Then  from  (27)  with 

S  =  na~  one  obtains  a  known  expression,  see,  e.g.,  [Bunchuk  1989,  Stanton  2000]  ms  =  CSR~  K/\n ) ,  where  however,  the 
coverage  parameter,  if  it  is  more  convenient,  can  be  replaced  by  the  average  volume  concentration  as  follows 

C,=C,H(s)/(v). 
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APPENDIX  II 
PROPAGATION  FACTOR 

If  scatterers  are  buried  in  the  sediment,  the  incident  and  scattered  waves  in  water  are  related  with  those  in  the  sediment  through 


corresponding  transmission  coefficients,  Wj  and  Ws  ,  and  plane  wave  propagation/attenuation  factors.  Such  approach  is  widely 

used  in  models  of  volume  scattering  in  sediments  [Ivakin  1981,  2000,  2001,  2004,  Jackson  2007]  and  in  detection  of  buried 
objects  in  the  seabed  [Hovem  1998].  Then  the  individual  scattering  cross-section  can  be  presented  in  the  form  (7)  with  the 
bistatic  propagation  factor  T  defined  as  follows 

T  =  \T?S  |V2 

(28) 

Ti,s  =  wi,s  cxp(iku  Pl  sz) 

(29) 

Wi  s  =  sin  Xi,s  /Osin  xUs  +  P,  J 

(30) 

Pi,s  =  (n2  -  cos2  Xi,s)112 

(31) 

Here  //  =  /)/  pw  is  the  sediment  to  water  density  ratio,  n  =  Cw  /  C  =  tl0  (1  +  iS)  is  the  sediment  complex  refraction  index,  or 
the  water  to  sediment  sound  speed  ratio,  Xi  s  are  grazing  angles  of  the  incident  and  scattered  waves  near  the  seabed,  kw  is  the 
sound  wave  number  in  water,  and  Wt  is  the  water-to-sediment  transmission  coefficient  in  the  direction  of  the  incident  and 
scattered  wave  (with  corresponding  subscript  index). 

Propagation  of  the  incident  and  scattered  waves  in  the  sediment  are  controlled  by  the  function  Pi  s,  =  P(Xi  s )  •  One  can  show 


P  =  ( n 2  -  cos?  z)W2  =  no(q  +  id /  q) 


q  =  q{%)  =  -J=  ,  V  =  1  - - d2 

V2  n~ 


The  real  part  of  equation  (32)  represents  in  fact  a  generalized  Snell’s  law  for  the  practically  important  case  where  the  refraction 
index  is  complex,  corresponding  to  absorbing  (lossy)  sediments.  It  defines  a  relationship  between  directions  of  propagation  of 
plane  waves  in  an  ideal  (lossless)  water  and  lossy  sediment,  i.e.,  directions  of  the  normal  to  the  phase  fronts  of  the  incident  and 
transmitted  (refracted)  plane  waves,  and  can  be  presented  in  the  form 


sin  Zt  =q(Xi)  =  —  RQ(Pi) 

n„ 


Correspondingly,  the  imagery  part  ofEq.(32),  through  the  Eq.  (29),  describe  exponential  attenuation  of  the  refracted  wave  along 
this  direction  and,  at  the  same  time,  gives  a  uniform  description  of  the  amplitude  of  both  refracted  and  evanescent  waves  in  the 
sediment  including  their  smooth  mutual  transformation  near  the  critical  grazing  angle  where  Xi  ~  Xc  ’  C°S  /K(  =  no  .  In 

particular,  at  the  critical  grazing  angle,  q{%c)  =  yfS  ,  and  at  small  subcritical  grazing  angles,  at^;  «  %c  ,  q  ~  8 / yjl  —  n~ 
Using  Eqs.  (28),  (32)  and  (34),  the  propagation  factor  becomes  of  the  form 


T(z)  =  T0e-*z 


(35) 
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I  .  2  _ 2 

where  Tg  =  | Wt  Ws  |  fl  ,  K  =  f:R\  / sin  y_tj  +  1  / sin  %ts  )  ,  and  [)  =  2 5kwn0  is  the  sound  attenuation  coefficient  in  the 
sediment. 


APPENDIX  III 

INDIVIDUAL  SCATTERING  FUNCTION 


Examples  of  exact  expression  for  the  individual  scattering  cross  section  in  the  case  of  an  elastic  solid  and  fluid  spheres  can  be 
found,  e.g.,  in  [Hickling  1962,  Morse  and  Ingard  1968,  Medwin  and  Clay  1998].  There  is  no  exact  solution  for  non-spherical 
particles  of  irregular  shape.  However,  there  are  approximate  solutions  for  low-frequency  (Rayleigh  scattering)  and  high- 
frequency  (geometry  acoustics)  limits  [Palmer  1996,  Bass  and  Fuks  1979].  There  are  also  wide  frequency  range  approaches  that 
can  be  numerically  efficient  in  the  case  of  smooth  enough  bodies  [Reeder  and  Stanton  2004].  For  many  practical  applications, 
there  have  been  developed  simple  empirical  formulae,  so-called  “high-pass”  expressions,  defined  in  such  way  that  they  satisfy 
known  low-frequency  (Rayleigh  scattering)  and  high-frequency  (geometry  acoustics)  limits  and  have  a  reasonable  continuous 
behavior  for  the  intermediate  frequencies  (see,  e.g.,  [Sheng  and  Hay  1988,  Stanton  1989,  Thorne  et  al  1995,  Thorne  and 
Backingham  2004]).  In  this  paper,  for  numerical  illustrations,  we  will  use  a  simplest  of  possible  “high-pass”  approximations  for 
the  scattering  function  in  the  form 

aQ(a)  =  a2F(ka) , 


x<xQ  =  (ax  / a 0) 
x>x„ 


1/4 


(37) 


where  a  is  a  “size”-parameter,  CC  and  GCm  are  dimensionless  parameters  defined  by  material  of  the  particle  and  their 

orientation.  Generally,  the  values  of  this  parameters,  as  well  as  the  meaning  of  the  “size”  of  a  non-spherical  particle,  should  result 
from  a  statistical  averaging  over  possible  orientation  of  the  particle.  In  this  paper,  following  other  “high-pass”  approaches,  we 
assume  that,  after  average  over  random  orientations,  the  parameters  in  (37)  can  be  approximated  by  those  for  spherical  particles 
with  the  same  equivalent  radius  (the  radius  of  a  sphere  having  the  same  volume).  As  such,  for  illustration  purposes  of  this  paper, 
we  will  use  parameters  of  a  fluid  sphere  [Morse  1968,  Medwin  1998,  Stanton  1989]. 


r; 


a0  =- 


(38) 


with  parameters  R0  and  R  defined  as  follows 


R 


o 


pc2  V  Pa  ~P  R  =  PaCa  ~  PC 
PaCl)  Pa+P/I'  PaCa+PC 


(39) 


where  Ca  and  pa  are  the  sound  speed  and  the  density  of  the  sphere  (representing  inclusions)  ,  C  and  p  are  the  sound  speed 
and  the  density  of  a  surrounding  fluid  medium  (representing  sediment). 


APPENDIX  IV 

SCATTERING  FROM  INCLUSIONS  WITH  POWER-LAW  SIZE  DISTRIBUTION 


Consider  the  inclusion  size  distribution  of  the  power  law  form  (or,  more  accurately,  “truncated”  power  law)  given  by  Eq.  (18). 
From  Eqs  (9)  and  (36),  one  obtains 


a  2 

M Q  =  m0a0  t B  -  a70 ~3  J a2~'' F (ka)da 

a\ 


(40) 
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This  function  represents  a  normalized  volume  scattering  coefficient  of  the  sediment.  Consider  its  frequency  dependence  using  a 
"high-pass”  approximation  for  the  individual  scattering  function  given  by  (37).  For  high  frequencies,  where 

kax  >  X0  =  /  0Co  )*/4 ,  it  can  be  presented  in  the  form 


M0=CCaoJ  3_r 


CL  i  Cl  n 


\ao  aoJ 


(41) 


where 

J  (x,  y )  =  f  f~xdt  =  - - — ,  J0  (x,  y )  =  In {y  /  x) 

{  7  (42) 


Note  that  in  this  case,  it  becomes  frequency  independent  regardless  the  power  exponent  of  the  size  distribution.  For  low 
frequencies,  where  ka 2  <  X() ,  we  have 


M0  =  a0(ka0)4J7 


ax 

Kao  ’  aoJ 


(43) 


In  this  case,  it  follows  the  Rayleigh  forth  power  frequency  dependence  regardless  the  power  exponent  of  the  size  distribution. 
At  intermediate  frequencies,  where  XQ  /  a-,  <  k  <  XQ  /  Ct{,  one  obtains 


(44) 


Mo  =a0(ka0)4J1, 


Vao  kaoJ 


+  a 


v  kao  ao  j 


For  power  exponents  in  the  range  3  <  y  <  7  and  wide  enough  size  distributions,  so  that  X0  /  Cl  2  «  k  «  XQ  /  ax, 

°°  a(7-y)/4  (y-3)/4 

M0  =  ( kaoy~3G(r X  G(r)=  f  x2^F(x)dx  *  - 

1  (7  -/)(/-  3)  (45) 

It  is  important  to  note  that  the  integral  in  (45)  has  little  sensitivity  to  a  particular  behavior  of  the  individual  scattering  function 
F  at  intermediate  x=ka,  which  makes  the  approximation  (45)  very  robust  and  stable  with  respect  to  variations  of  the  scatterer 
type,  its  shape,  material,  etc.  Therefore,  using  a  simple  high-pass  expression  (37),  for  which  the  approximate  equality  in  (45) 
becomes  exact,  seems  to  be  reasonable  and  won’t  affect  the  result  significantly.  If  more  realistic  individual  scattering  functions 
are  used,  the  reasonable  accuracy  of  the  approximate  equality  in  (45)  will  still  hold.  All  the  material  parameters  of  the  inclusion 
are  included  in  the  factors  (X0  and  Cl.f  .  To  define  these  parameters,  for  illustration  purposes,  Eqs.  (38)  and  (39),  corresponding 

to  a  fluid  sphere  case,  are  used  in  this  paper.  The  result,  the  dependence  of  the  normalized  volume  scattering  coefficient  M  as  a 
function  of  the  dimensionless  frequency  kao  for  different  values  of  the  power  exponents  y  is  shown  in  Fig.3.  Another  general 

result  is  important  to  emphasize  here.  According  to  (45),  in  the  case  of  the  wide  power  law  size  distribution  of  inclusions,  and  for 
the  intermediate  frequency  range,  the  frequency  dependence  of  the  volume  scattering  coefficient  is  very  sensitive  with  respect  to 
parameters  of  the  size  distribution.  In  particular,  it  also  obeys  power  law  and  its  power  exponent  is  linearly  related  to  the  one  for 
the  size  distribution.  This  can  be  an  important,  particularly,  for  remote  acoustic  sensing  of  the  sediment. 
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Frequency  Dependence  of  Phase  Speed,  Group 
Speed,  and  Attenuation  in  Water-Saturated  Sand: 

Laboratory  Experiments 

Jean-Pierre  Sessarego,  Anatoly  N.  Ivakin,  and  Didier  Ferrand 


Abstract — Sound  propagation  in  water-saturated  sandy  sedi¬ 
ments  was  studied  under  controlled  laboratory  conditions  in  the 
0.1-1.3-MHz  frequency  range.  In  the  “low-frequency  domain” 
(100-200  kHz),  the  results  obtained  were  consistent  with  classical 
Biot  theory,  which  predicts  a  positive  phase-speed  dispersion, 
but  at  higher  frequencies  (/  >  500  kHz),  a  strong  negative 

dispersion  of  the  phase  speed  was  observed,  which  is  in  apparent 
contradiction  with  the  Biot  theory  predictions.  In  addition,  the 
attenuation  coefficient  was  found  to  increase  nonlinearly  with 
frequency  from  0.5  to  1.3  MHz,  whereas  Biot  theory  predicted  an 
almost  linear  increase  of  this  coefficient.  The  unexpected  behavior 
of  the  phase  speed  and  the  attenuation  coefficient  observed  in 
these  sandy  sediments  in  the  very  high-frequency  domain  may  be 
attributable  to  scattering  by  the  sand  particles,  which  becomes 
the  predominant  mechanism  involved  when  the  size  of  particles  is 
comparable  to  the  wavelength. 

Index  Terms — Attenuation,  Biot  model,  dispersion,  group  speed, 
phase  speed,  sediment  acoustics. 


I.  Introduction 

THE  understanding  of  sound-seafloor  interactions  is  a 
challenging  issue  that  requires  a  thorough  knowledge  of 
the  intrinsic  properties  of  the  sediments.  In  particular,  the  phase 
speed  and  the  attenuation  coefficient  for  the  longitudinal  waves 
in  sandy  sediments  are  parameters  that  need  to  be  accurately 
determined  because  they  are  used  as  input  data  to  most  acoustic 
models  for  sea-bottom  reflection  and  backscattering.  These  pa¬ 
rameters  have  been  measured  in  marine  sand  by  many  scientists 
during  the  last  20  years  (see,  e.g.,  [  1  ]— [8]).  The  measurements 
have  been  performed  under  both  laboratory  and  natural  (at  sea) 
conditions,  but  the  frequency  range  studied  has  usually  been  rel¬ 
atively  low  (<300  kHz).  At  these  low  frequencies,  sediment  can 
be  modeled  as  a  porous  medium,  and  corresponding  theories 
(such  as  that  developed  by  Biot  and  other  related  approaches, 
e.g.,  “effective  density  fluid  model”  developed  by  Williams  [7]) 
effectively  explain  the  main  features  of  sound  propagation  in 
these  media.  One  of  the  main  assumptions  adopted  in  “low-fre¬ 
quency  effective  medium”  models  is  that  the  particle  size  in 
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the  medium  is  much  smaller  than  the  acoustic  wavelength.  At 
these  low  frequencies  (<300  kHz),  the  p- wave  in  the  sediment 
is  known  to  show  positive  phase-speed  dispersion.  In  this  study, 
the  negative  dispersion  of  the  phase  speed  was  found  to  occur 
at  higher  frequencies  (>500  kHz).  This  finding  suggests  that 
in  the  high-frequency  domain,  scattering  by  sand  particles  be¬ 
comes  the  predominant  mechanism  responsible  for  the  negative 
phase-speed  dispersion. 

II.  Preliminary  Measurements 

Before  the  acoustical  experiments,  preliminary  measure¬ 
ments  were  carried  out  to  determine  the  geotechnical  parame¬ 
ters  of  the  sand  used  in  this  study.  These  measurements  were 
performed  at  The  Laboratoire  de  Mecanique  et  d’Acoustique 
(LMA,  Marseille,  France)  with  a  well-sorted  sand  having  a 
mean  grain  diameter  of  245  pm  and  a  relatively  narrow  size 
distribution  with  d.{ 0.1)  =  165  pm  and  d{ 0.9)  =  364  pm. 
The  parameter  d.(W )  is  defined  such  that  W  is  the  ratio  of  the 
volume  (weight)  of  the  grains  with  diameters  less  than  d(W ) 
to  the  total  volume  of  all  the  grains  (dry  weight  of  the  sediment 
sample).  To  determine  the  wet  density  of  the  sediment,  two 
different  measurements  were  performed.  First,  the  density 
was  measured  without  any  settling  of  the  sediment  and  the 
value  obtained  was  p  —  1977  kg/m3±  10  kg/m3.  Second,  the 
density  was  measured  once  the  sediment  had  settled  (after 
manual  shaking).  The  value  obtained  in  this  case  was  p  = 
2004  kg/m3±  10  kg/m3.  The  porosity,  which  is  defined  as 
the  ratio  of  the  sediment  pore  volume  to  its  total  volume  was 
found  to  be  /?  =  0.365  ±  0.02.  This  result  is  consistent  with 
many  measurements  on  the  porosity  of  sand  with  a  similar 
granulometry. 

The  permeability  and  the  tortuosity  of  the  sediment  were  also 
measured  at  the  LMA.  Tortuosity  was  measured  by  placing  elec¬ 
trical  conductivity  probes  in  a  small  box  containing  sand  satu¬ 
rated  with  a  saline  solution  (30%o)-  The  ratio  between  the  re¬ 
sistivity  of  the  saline  solution  and  the  resistivity  of  the  satu¬ 
rated  sand  gives  the  formation  factor  F,  which  appears  in  the 
equation  a  =  Ffd  relating  the  tortuosity  a  to  the  porosity  (3. 
The  value  of  the  tortuosity  obtained  with  this  technique  was 
a  =  1.31  ±0.13.  To  measure  the  permeability,  a  similar  method 
to  that  described  by  Darcy  was  used,  which  gave  the  following: 
k  =  1.65  10_11m2±1.2  10-n.  It  can  be  seen  here  that  the  un¬ 
certainty  is  very  high,  but  this  is  quite  usual  with  such  measure¬ 
ments;  see,  for  example,  [7]  and  [9], 

To  compute  the  phase  speed  and  the  attenuation  of  the  longi¬ 
tudinal  wave  using  the  Biot  model  [10],  [1 1],  these  experimental 
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TABLE  I 

Parameters  Used  to  Compute  the  Phase  Speed  and  the  Attenuation  With  the  Biot  Model 


Input  parameters  for  the  Biot  model 

Bulk  modulus  of  grains 

3.2  10'°Pa 

Permeability  of  sediment 

*1.6510'"  m2+/-1.2  10'" 

Tortuosity  of  sediment 

*1.31  +/-0.13 

Porosity  of  sediment 

*0.365+/-0.02 

Viscosity  of  water 

0.001  kg.m'.s'1 

Mass  density  of  sand  grains 

2690  kg/m3 

Bulk  modulus  of  water 

2.19  109  Pa 

Mass  density  of  water 

998  kg/m3 

Shear  modulus  of  sediment  frame 
Real  part 

Imaginary  part 

2.92  107Pa 
-1.8  106  Pa 

Bulk  modulus  of  sediment  frame 
Real  part 

Imaginary  part 

4.36  107  Pa 
-2.08  1 06  Pa 

The  symbol  (*)  indicates  the  measured  values. 


TABLE  II 

Characteristics  of  the  “Panametrics”  Transducers  Used  in  the  Experiments 


Central  frequency 
and  type  of 
transducer 

Frequency 
band  with 
@-3dB 

Diameter  D 

2A./D 

Far  Field 
Distance 
N=D2Aa 

Angular 
Beamwidth  203 
@-3dB 

250  kHz  (V3536) 

120  kHz 

40  mm 

0.3 

6.5  cm 

1.5° 

250  kHz 

25  mm 

0.24 

5.2  cm 

6° 

1MHz  (V302) 

480  kHz 

25  mm 

0.12 

10.5  cm 

3° 

phase-speed,  group-speed,  and  attenuation  measurements  ob¬ 
tained  with  the  three  transducers. 

III.  Sound-Speed  Measurements 

Fig.  4  gives  a  schematic  diagram  of  the  experimental  setup. 
This  setup  consisted  of  a  rectangular  PVC  box[AU Pie  CIS  e 
define  PVC]  (l  =  1 5  cm,  L  =  20  cm,  H  =  15  cm)  in  which 
ten  transducers  (four  working  at  500  kHz,  four  at  1  MHz,  and 
two  at  250  kHz)  were  fixed  so  that  they  were  perfectly  aligned 
face  to  face  in  pairs.  The  acoustic  method  used  to  measure  the 
distances  between  transducers  as  accurately  as  possible  will  be 
described  below. 

Based  on  the  far-field  distances  given  in  Table  II,  it  can  be 
seen  in  Table  III  that  far-field  conditions  were  fully  satisfied  in 
the  case  of  the  two  lower  frequency  sources,  whereas  they  were 
only  partly  fulfilled  in  that  of  the  1-MHz  source.  This  situation 
did  not  affect  the  phase-speed  measurements.  It  could  only  in¬ 
troduce  minor  errors  in  the  attenuation  measurements,  but  the 
values  measured  at  1  MHz  were  so  high  that  this  effect  was  taken 
to  be  negligible. 

Transmission  experiments  were  carried  out  separately  in 
fresh  water  and  sand.  The  transmission  in  fresh  water  was 
used  to  measure  the  distance  between  the  transducers.  For 
this  measurement,  it  was  assumed  that  the  phase  speed  can  be 
exactly  determined  as  long  as  the  water  temperature  is  known 
and  also  that  sound  absorption  in  fresh  water  can  be  neglected. 
It  is  worth  noting  here  that  the  sound  attenuation  in  fresh  water 


data  were  introduced  into  the  Biot  equations  and  standard  values 
available  in  the  literature  [7],  [12]  were  used  for  the  other  model 
parameters,  such  as  the  grain  density,  the  bulk  and  shear  mod¬ 
ulus  of  the  sediment  frame,  the  bulk  modulus  and  mass  density 
of  the  sand  grains,  and  the  viscosity  and  bulk  modulus  of  water. 
All  the  parameters  used  as  inputs  in  the  Biot  model  are  summa¬ 
rized  in  Table  I. 

In  our  acoustical  measurements  of  the  phase  speed  and  the 
attenuation  coefficient,  we  used  wet  and  carefully  degassed 
sand.  Because  the  presence  of  small  air  bubbles  in  sands  is 
known  to  considerably  affect  the  results,  the  sand  used  in  our 
experiments  was  immersed  in  fresh  water  for  several  months 
before  the  experiments  and  careful  degassing  operations  were 
performed  for  several  weeks.  The  acoustical  measurements 
were  then  performed  on  the  water-saturated  sand  mixture  after 
pouring  it  into  a  rectangular  box,  making  sure  that  the  sand  was 
never  exposed  to  air. 

To  fully  cover  the  wide  frequency  range  chosen,  three  dif¬ 
ferent  “panametrics”  transducers  were  used.  Their  characteris¬ 
tics  (the  bandwidth,  the  diameter  of  the  piezoelectric  ceramic, 
the  far-field  distance,  and  the  beamwidth)  are  given  in  Table  II. 
The  frequency  spectra  of  the  transducers  are  shown  in  Figs.  1-3. 
It  can  be  seen  from  these  figures  that  each  of  the  three  trans¬ 
ducers  used  covered  a  very  broad  frequency  band.  Their  fre¬ 
quency  spectra  overlapped,  and  together,  they  covered  the  whole 
0.1-1.5-MHz  frequency  range.  It  can  be  noted  that  the  over¬ 
lapping  frequencies  were  used  to  check  the  consistency  of  the 
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Fig.  1 .  Time  response  of  the  250-kHz  transducer  and  the  corresponding  spectrum. 


Fig.  2.  Time  response  of  the  500-kHz  transducer  and  the  corresponding  spectrum. 


Fig.  3.  Time  response  of  the  1-MFlz  transducer  and  the  corresponding  spectrum. 


is  about  0.1  dB/m  at  1  MHz,  which  is  negligible  compared  with 
the  sound  attenuation  in  sediment  (>  100  dB/m  at  1  MHz). 

After  immersing  the  PVC  box  in  a  large  water  tank,  an 
acoustic  method  was  used  to  measure  the  distances  between 
the  transducers.  This  method  consisted  first  in  measuring  the 
delay  between  the  transmitted  and  received  short  pulses  in  the 


transmitting  mode.  The  distance  d  between  the  transducers  was 
then  calculated  using  the  simple  equation:  d  —  Cwt ,  where 
t  is  the  delay  between  the  transmitted  impulse  and  received 
signal,  and  Cw  is  the  sound  speed  in  water.  In  the  case  of  sound 
propagation  in  water,  there  is  no  dispersion  and  the  beginning  of 
the  signal  received  can  be  measured  with  good  accuracy,  which 
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Fig.  4.  Schematic  diagram  of  the  box  equipped  with  transducers. 


TABLE  III 

Distances  Between  the  Transducers  During  Measurements 


Transducers 

Min  distance  (cm) 

Max  distance  (cm) 

250  kHz 

12.78+/- 0.02 

500  kHz 

10.89 +/-0.02 

15.88 +/-0.02 

1  MHz 

1 1.76  +/-0.02 

15.88 +/-0.02 

gives  a  satisfactory  estimate  of  the  propagation  time.  The  value 
of  Cw  was  obtained  using  the  relation  given  by  Grosso  [14]  for 
the  sound  speed  in  fresh  water  as  a  function  of  temperature 

CW(T)  =  1402.388  +  5.03711T  -  0.0580852T2 

+3.342  10~4T3  -  1.478  10~6T4  +  3.15  10-9T5  (1) 

where  T  is  the  temperature  in  Celsius  degrees.  This  relation  is 
valid  for  0°  <  T  <  95°.  It  gives  the  sound  speed  in  the  water 
with  an  uncertainty  of  ±0.015  m/s. 

Some  minor  corrections  were  made  to  take  into  account  the 
transition  time  in  the  electrical  circuit  and  the  propagation  time 
through  the  protection  membrane  of  the  transducers.  This  addi¬ 
tional  delay  was  measured  by  bringing  the  transducers  into  con¬ 
tact  (using  Panametrics  gel  between  the  two  transducer  faces  to 
ensure  perfect  coupling).  With  each  pair  of  transducers  in  con¬ 
tact,  the  delay  rc  was  determined  and  used  to  correct  the  distance 
measurements  as  follows: 


d  =  Cw(t  -  tc).  (2) 

Because  the  error  in  the  time  measurements  was  less  than 
0.15  fis,  this  method  gives  an  acoustical  estimate  of  the  dis¬ 
tance  between  the  transducers  with  an  accuracy  of  more  than 
0.2  mm.  It  is  worth  noting  here  that  the  overall  accuracy  of  the 
phase-speed  measurements  in  the  sand  was  mainly  limited  by 
the  lack  of  precision  in  the  distance  measurements  [12].  The 
results  obtained  are  summarized  in  Table  III. 

During  the  calibration  process,  the  amplitudes  Aw\(f)  and 
71,1,2 (/)  of  the  signals  transmitted  over  distances  d\  and  d 2  in 
water  were  also  measured  as  a  function  of  the  frequency.  When 
the  calibration  measurements  in  fresh  water  were  completed,  the 
PVC  box  was  filled  with  water-saturated  sand. 


A.  Phase-Speed  Measurements 

1 )  Broadband  Measurements:  In  this  wideband  experiment, 
very  short  pulses  were  used  in  the  transmission  mode.  The  elec¬ 
tric  signals  sent  to  the  transducer  were  very  short  pulses  of  high 
voltage  (>300  V)  delivered  by  a  “panametrics”  pulser/receiver 
5055PR.  The  corresponding  acoustical  signals  were  recorded  on 
a  sampling  oscilloscope  and  transmitted  to  a  PC  for  processing. 
To  obtain  the  phase  speed,  the  following  relation  was  used  [15]: 

2itvd  ... 

CpiAt')  =  - y. - j ~ (3) 

( imag  (5water(*')) 

2ttut  +  arg  - - 77; - -t-tt 

\  imag  (osand(;z)) 

where  5sanri  and  SwatGr  are  the  Fourier  transform  of  the  signal 
received  after  propagating  through  sand  and  water,  respectively, 
d  is  the  distance  between  the  transmitter  and  the  receiver,  and  t 
is  the  time  of  arrival  of  the  signal,  which  has  been  defined  above. 
In  this  case,  only  one  distance  was  required  for  phase-speed 
measurements. 

The  phase  speed  is  given  in  Fig.  5  as  a  function  of  frequency. 
The  curves  in  this  figure  give  the  results  of  measurements  ob¬ 
tained  with  the  three  wideband  transducers.  In  the  region  where 
the  curves  overlap,  reasonable  agreement  can  be  seen  to  exist  be¬ 
tween  the  data  obtained  with  the  three  transducers.  The  phase- 
speed  measurement  error  was  about  +/— 10  m/s.  The  most  im¬ 
portant  result  obtained  here  is  the  noticeable  decrease  in  the 
phase  speed  observed  in  the  0.5-1. 3-MHz  frequency  range.  In 
this  experiment,  the  temperature  of  water  was  20.9  °C  giving  a 
sound  velocity  [AU:  SOUtld  Speed? ]of  1485  m/s. 

2)  Tone  Bursts:  It  was  then  proposed  to  confirm  the  results 
obtained  in  the  above  wideband  measurements  using  narrower 
band  signals.  For  this  purpose,  another  set  of  measurements 
were  performed  using  tone  bursts  with  various  frequencies.  By 
measuring  the  time  of  arrival  of  the  signals  transmitted  in  sand 
along  two  distances  d\  and  d2,  the  phase  speed  was  obtained  as  a 
function  of  the  frequency.  In  this  case,  the  electrical  signal  deliv¬ 
ered  to  the  transducer  was  a  sinusoidal  wave  train  of  1 0, 1 5,  or  20 
periods  depending  on  the  frequency  used.  At  “low  frequencies” 
(<200  kHz),  ten  periods  were  used  for  the  sinusoidal  signal, 
while  at  “high  frequencies”  (>500  kHz),  the  number  of  periods 
was  set  at  20.  The  working  frequency  was  increased  in  20-kHz 
increments  to  cover  the  whole  frequency  range  of  interest.  All 
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Phase  speed  as  a  function  of  frequency 


Fig.  5.  Phase-speed  measurements  in  sand. 


Fig.  6.  Phase-speed  measurements  obtained  at  two  different  depths  on  different 
days  (500-kHz  panametrics  transducer). 


the  signals  transmitted  over  the  two  distances  d\  and  d.2  were 
sent  to  an  analog-to-digital  converter  before  been  recorded  on 
a  PC.  To  estimate  the  phase  speed,  the  classical  zero-crossing 
method  of  measurement  was  used  and  the  sound  speed  was  then 
obtained  from  the  relation 


where  T\  and  T2  are  the  delays  associated  with  a  given  number 
of  zero  crossings  for  distances  d\  and  d. 2.  It  was  checked, 
for  example,  that  with  high-frequency  signals,  the  tenth  zero 
crossing  was  associated  with  a  portion  of  the  signal  where 
the  instantaneous  frequency  was  constant.  The  velocity  mea¬ 
surement  [AU:  Can  "velocity"  be  changed 
to  "speed"  throughout? /error  obtained  using  this 
method  was  estimated  to  be  approximately  ±5  m/s. 

Curves  giving  the  phase  speed  as  a  function  of  the  frequency 
are  shown  in  Fig.  6  in  the  180-700-kHz  frequency  range.  Mea¬ 
surements  were  performed  on  various  days  at  various  depths 
to  determine  whether  the  settlement  of  the  sediment  had  any 
effects.  The  data  given  in  Fig.  6  show  that  at  the  same  trans¬ 
ducer-embedding  depth  (10  cm),  the  results  obtained  before  and 
after  a  period  of  two  weeks  were  very  similar  (February  17  to 
March  2).  It  can  also  be  seen  that  the  phase  speed  measured  at  a 
depth  of  3  cm  was  practically  parallel  to  that  obtained  at  a  depth 
of  10  cm.  This  finding  proves  that  the  same  dispersion  law  ap¬ 
plied  at  two  different  depths.  In  addition,  the  phase  speed  at  3  cm 
below  the  sand-water  interface  was  20  m/s  lower  than  that  mea¬ 
sured  at  a  depth  of  10  cm.  This  difference  might  be  attributable 
to  the  fact  that  the  near-surface  sediment  was  more  saturated 
by  water  than  the  deeper  layer,  and  that  there  was  probably  a 
sound-speed  gradient  in  the  sand  bottom.  It  was  also  checked 
that  from  February  15  to  March  6,  the  temperature  measured 
both  in  the  water  tank  and  in  the  sediment  was  constant  and 
equal  to  22  °C  giving  a  sound  speed  of  1488.3  m/s  in  the  water. 

The  results  obtained  with  the  500-  and  250-kHz  transducers 
were  compared  in  the  same  way  (see  Fig.  7).  Although  the 
measurements  were  performed  on  different  days,  the  values 
obtained  with  the  two  transducers  were  consistent.  Positive 


phase-speed  dispersion  can  be  observed  in  the  low-frequency 
domain  (100-300  kHz),  which  is  in  agreement  with  the  Biot 
predictions,  but  at  higher  frequencies,  negative  phase-speed 
dispersion  begins  to  occur. 

To  confirm  the  decrease  in  the  phase  speed  observed  here, 
another  set  of  measurements  were  performed  at  higher  fre¬ 
quencies  (0.5-1 .2  MHz).  These  measurements  were  carried  out 
using  long  sinusoidal  pulses  (20  periods).  The  results  are  given 
in  Fig.  8  and  compared  with  those  obtained  with  the  500-kHz 
transducer.  Both  curves  show  the  occurrence  of  a  similar  de¬ 
crease  in  the  phase  speed  with  the  frequency.  Some  differences 
observed  between  the  results  obtained  with  the  two  different 
transducers  (at  the  overlap  frequencies)  are  within  the  range  of 
uncertainty  defined  above.  It  is  worth  noting  that  the  negative 
dispersion  behavior  of  the  phase  speed  at  high  frequencies 
(>400  kHz)  was  also  observed  by  Moussatov  et  al.  [16]  in  the 
case  of  sediments  and  glass  beads  of  about  the  same  sizes. 

The  measurements  were  performed  with  the  transducers 
placed  3  cm  below  the  water-sand  interface;  during  these 
measurements,  the  temperature  in  both  the  water  tank  and  the 
sand  box  was  around  17.5  °C. 
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Fig.  8.  Phase-speed  measurements  (500  and  1  MHz). 


Fig.  9.  Phase-speed  measurements  as  a  function  of  the  frequency  (comparison 
with  Biot  model  [10],  [1 1]). 


To  compare  the  experimental  data  with  model  predictions,  it 
has  often  been  proposed  to  normalize  the  sound  speed  in  sed¬ 
iment  with  respect  to  the  sound  speed  in  water  at  the  time  of 
the  measurements.  The  results  obtained  in  the  0. 1-1 ,5-MHz  fre¬ 
quency  range  are  given  in  Fig.  9.  In  this  figure,  the  phase  speed 
computed  using  the  Biot  model  [10],  [11]  for  water-saturated 
sand  is  given  by  the  dashed  line.  The  values  of  the  parameters 
used  with  the  Biot  model  are  given  in  Table  I. 

The  following  comments  can  be  made:  At  high  frequencies 
(>500  kHz),  the  results  of  our  experiments  show  that  the  phase 
speed  no  longer  obeys  the  Biot  model.  The  phase  speed  was 
found  to  decrease  strongly  as  the  frequency  increased.  This  de¬ 
crease  was  of  the  order  of  40  m/s  in  the  0.5-1 -MHz  frequency 
range.  This  is  in  apparent  contradiction  with  the  Biot  model  pre¬ 
dictions,  but  it  can  be  noted  that  at  these  high  frequencies,  the 
size  of  the  sand  particles  becomes  comparable  to  that  of  the 
wavelength,  and  it  is,  therefore,  not  surprising  that  some  scat¬ 
tering  processes  begin  to  occur. 


B.  Group-Speed  Measurements 

We  also  measured  the  group  speed  by  making  a  correlation 
between  the  signals  received  in  the  sand  at  distances  d\  and 
d2.  The  maximum  value  of  the  correlation  function  gives  the 
time  of  arrival  of  the  energy,  from  which  the  group  speed  can  be 
calculated  by  using  the  relation 


(5) 


where  r  is  the  delay  obtained  from  the  position  of  the  max¬ 
imum  value  of  the  correlation  function.  Fig.  10  gives  the  group 
speed  as  a  function  of  frequency  from  0.2  to  1 .2  MHz.  The  group 
speed  curve  shows  a  strong  dispersion,  which  is  even  more  pro¬ 
nounced  than  in  the  case  of  the  phase  speed.  The  mean  rate  of 
decrease  in  the  group  speed  was  60  m/s  between  frequencies  of 
0.4  and  1.2  MHz. 

The  phase  speed  and  the  group  speed  are  shown  in  Fig.  1 1 
on  the  same  scale.  It  can  be  seen  from  this  figure  that  the  group 
speed  was  always  lower  than  the  phase  speed.  The  higher  the 
frequency,  the  greater  is  the  difference  between  the  two  curves. 
This  result  can  be  explained  using  the  classical  relationship  be¬ 
tween  the  group  speed  and  the  phase  speed 


6/yr  —  Cpli 


dCph 

dk 


(6) 


where  k  is  the  wave  number. 


IV.  Attenuation  Measurements 

To  measure  the  sound  attenuation  in  the  sediment,  the  same 
arrangement  as  that  described  above  for  measuring  the  phase 
speed  was  used.  The  amplitude  of  the  wideband  signals  trans¬ 
mitted  in  both  water  and  sand  was  measured.  A  Fourier  trans¬ 
form  was  applied  to  the  signals  received  to  determine  the  am¬ 
plitude  spectrum  in  a  very  wide  frequency  range.  The  attenua¬ 
tion  in  the  sediment  was  then  determined  as  a  function  of  the 
frequency  from  the  ratio  between  the  magnitude  of  the  Fourier 
transforms  obtained  in  the  case  of  sediment  and  water.  It  was 
again  assumed  that  there  was  no  attenuation  in  fresh  water  (the 
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Fig.  1 1 .  Comparison  between  the  phase  speed  and  the  group  speed  in  the 
0.2-1.2-MHz  frequency  band. 


reference  medium).  This  method  was  used  previously  by  Buck¬ 
ingham  [17],  The  attenuation  is  given  by  the  following  relation: 


a(i /)  —  8.686-Log10 


(  Nw{u)\ 
\Ns(u)J 


where 


(7) 


Nw{  v)  = 


Ay)l  (/Z) 

Aw2(v)' 


N.(v)  = 


As\  jv) 

As  2{v)' 


Aw \(v),  AW2(v),  Asi(u),  and  As2{v)  are  the  amplitude  spectra 
of  the  wideband  signals  received  in  water  and  sand  at  distances 
d\  and  <i2-  The  attenuation  coefficient  a  is  given  in  decibels  per 
unit  distance. 

The  results  of  the  measurements  obtained  with  the  1-MHz 
transducer  are  shown  in  Fig.  12.  The  maximum  frequency  used 
to  calculate  (7)  was  1 .3  MHz  because  at  higher  frequencies, 
no  signals  could  be  measured,  because  the  attenuation  was  too 
strong.  In  this  figure,  each  data  point  is  the  mean  value  of  ap¬ 
proximately  1 5  measurements  and  the  error  bars  denote  the  stan¬ 
dard  deviations  of  these  measurements. 

The  results  of  the  attenuation  obtained  using  the  Biot  model 
are  also  given  in  Fig.  12.  As  was  to  be  expected,  this  model 
does  not  accurately  describe  the  attenuation  in  the  very  high 
frequency  range.  We  have  also  added  some  results  obtained  by 
Seifert  et  al.  [18]  on  sand  particles  with  a  mean  diameter  of 
230pm.  Our  measurements  gave  a  lower  attenuation,  but  the 
size  distribution  of  the  sand  particles  was  not  exactly  the  same 
in  both  studies.  In  addition,  while  correcting  this  manuscript, 
we  received  the  paper  by  Lee  et  al.  [19]  and  we  have  also  added 
these  authors’  results  to  Fig.  12.  The  latter  authors  obtained  sim¬ 
ilar  results  to  those  obtained  here,  but  Lee  et  al.  observed  higher 
sound  attenuation  levels,  probably  because  the  mean  particle  di¬ 


Fig.  12.  Sound  attenuation  measurements  versus  frequency.  Comparison  with 
the  Biot  model  computed  in  the  case  of  fresh  water  and  with  data  by  Seifert  [18] 
and  Lee  et  al.  [19].  (The  stars  correspond  to  the  scattering  attenuation  obtained 
by  Seifert  and  the  circles  to  data  by  Lee  et  al.) 


ameter  of  the  sand  they  used  was  425  pm,  which  is  much  larger 
than  that  of  our  sand  (245  pm).  This  difference  suggests  that 
the  attenuation  may  depend  not  only  on  the  frequency  but  also 
on  the  particle  size.  From  all  these  results,  it  can  be  concluded 
that  at  high  frequencies  (>500  kHz),  the  attenuation  coefficient 
increases  nonlinearly  with  frequency  (the  attenuation  seems  to 
be  proportional  to  /2),  which  is  probably  due  to  both  classical 
viscous  losses  (Biot  theory)  and  scattering  by  sand  particles. 
Ohkawa  reached  similar  conclusions  [20],  He  mentioned  that 
the  scattering  effect  is  the  main  process  involved  at  frequencies 
above  200  kHz  in  sand  with  a  mean  particle  diameter  of  230  pm. 

V.  Conclusion 

It  can  be  concluded  that  although  the  Biot  model  can  be  used 
to  predict  the  phase  speed  and  the  attenuation  in  the  low-fre¬ 
quency  domain  in  the  case  of  a  medium  sand,  this  method  is 
not  appropriate  for  dealing  with  the  high-frequency  part  of  the 
spectrum  (>500  kHz)  where  the  sound  scattering  by  sand  par¬ 
ticles  becomes  the  main  contribution.  At  the  higher  frequencies 
tested,  a  strong  decrease  in  both  phase  speed  and  group  speed 
and  a  strong  increase  in  the  attenuation  were  observed.  Upon 
comparing  our  results  with  those  of  Lee  et  al.  [19],  it  is  sug¬ 
gested  that  the  scattering  by  sand  particles  contributes  to  gen¬ 
erating  both  attenuation  and  sound-speed  “anomalies.”  These 
scattering  effects  can  be  very  strong  and  can  considerably  affect 
sound  reflection  and  backscattering  by  sediments. 

Finally,  it  can  be  mentioned  that  for  larger  sand  particles 
(coarse  sand,  for  example)  scattering  can  occur  at  lower  fre¬ 
quencies  and  fall,  in  particular,  to  the  frequency  range  of  most 
high-resolution  sonars  used  for  sea-bottom  imaging. 
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High-Frequency  Sound  Reflection  by  Water- Saturated 

Sediment  Interfaces 

Jean-Pierre  Sessarego,  R.  Guillermin,  and  Anatoly  N.  Ivakin 


Abstract — Sound  reflection  by  water-saturated  sands  and  glass 
beads  with  a  flattened  surface  was  studied  under  controlled  labora¬ 
tory  conditions  in  a  wide  frequency  range,  from  200  kHz  to  7  MHz. 
In  the  “low-frequency”  domain  and  in  the  case  of  medium  sand, 
the  reflected  sound  level  was  found  to  be  in  good  agreement  with 
both  classical  sonar  measurements  and  classical  theories  of  reflec¬ 
tion  developed  for  fluid-porous  media  (this  reflected  level  is  prac¬ 
tically  independent  of  the  frequency);  as  the  frequency  increases, 
a  large  decrease  in  the  reflected  level  occurs,  possibly  due  to  inco¬ 
herent  scattering.  In  the  very  high-frequency  regime  (>3  MHz), 
the  sound  level  measured  was  more  than  20  dB  below  the  classical 
level,  and  it  remained  constant  at  higher  frequencies.  Similar  ex¬ 
periments  were  carried  out  with  coarse  sand  to  study  the  effect 
of  grain  size  on  the  reflection  loss.  The  same  behavior  with  only  a 
frequency  shift  was  observed.  These  effects  were  confirmed  by  re¬ 
peating  the  experiments  with  glass  beads  of  two  sizes.  The  anom¬ 
alies  observed  in  the  reflected  levels  measured  seem  to  be  directly 
connected  to  the  ratio  between  the  grain  size  and  the  wavelength. 
One  of  the  main  conclusions  reached  in  this  study  was  that  for  very 
coarse  sand  and  gravels,  the  effects  of  the  granular  structure  of  the 
bottom  should  not  be  neglected,  even  at  the  usual  sonar  frequen¬ 
cies. 

Index  Terms — Granular  media,  reflection  loss,  scattering  by 
grains,  sediment  acoustics,  sound  reflection. 


I.  Introduction 

MEASUREMENTS  of  sound  reflection  from  sandy  bot¬ 
toms  have  been  carried  out  for  years  in  a  relatively  low- 
frequency  regime  (ranging  between  a  few  kilohertz  to  few  hun¬ 
dreds  of  kilohertz)  by  many  authors  [  1  ]— [3]  and  [4].  This  interest 
in  the  low-frequency  domain  is  mainly  because  most  sonars 
work  at  relatively  low  frequencies  (<300  kHz).  Therefore,  it 
has  often  been  assumed  that  sound  reflection  is  affected  only 
by  the  surface  roughness  and,  in  the  case  of  very  low-frequency 
sonars,  by  volume  inhomogeneities;  the  various  models  which 
have  been  developed  up  to  now  for  acoustical  backscattering, 
therefore,  do  not  include  scattering  by  sand  grains  because  it 
is  generally  assumed  that  at  these  usual  sonar  frequencies,  the 
grain  size  does  not  affect  the  sound  reflection  levels.  The  ap¬ 
parent  lack  of  marine  applications  in  the  high-frequency  domain 
may  explain  why  very  few  sound  reflection  measurements  have 
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TABLE  I 

Statistical  Characteristics  of  Sands  and  Glass  Beads 


Medium  Sand 

Coarse  Sand 

Small  Beads 

Large  Beads 

Mean  grain  size 

245  |um 

1.50  mm 

200  |rm 

3  mm 

d(0. 1  )=I  65  |um 

d(0.9)=364|im 

d(0.1)=1.3mm 

d(0.9)=2.6mm 

been  carried  out  above  500  kHz.  However,  some  interesting  re¬ 
sults  have  been  obtained  under  laboratory  conditions  by  Nolle 
[5]  and  Williams  [6],  and  more  recently,  by  Greenlaw  et  al.  [7] 
and  Ivakin  et  al.  [8]— [10], 

To  better  understand  the  complex  mechanisms  involved  in 
sea  bottom  reflection,  experiments  were  conducted  at  Labora¬ 
toire  de  Mecanique  et  d’Acoustique  (LMA),  Centre  National 
de  la  Recherche  Scientifique  (CNRS,  Marseille,  France)  with 
water-saturated  sands  and  glass  beads,  where  all  the  experi¬ 
mental  parameters  were  perfectly  controlled.  The  aim  of  this 
study  was  to  measure  the  reflection  loss  as  a  function  of  the 
frequency  over  a  very  wide  frequency  range,  from  200  kHz  to 
7  MHz,  to  see  whether  a  transition  zone  between  continuous 
and  granular  behavior  could  be  detected  in  the  reflection  levels 
measured.  Four  wideband  Panametrics  transducers  were  used  to 
cover  this  very  wide  frequency  range. 

n.  The  Experiment 

Because  the  main  goal  of  this  study  was  to  determine  the  re¬ 
flection  loss  for  different  types  of  “sediments”  as  a  function  of 
the  frequency,  several  experiments  were  performed  on  two  types 
of  sand  (medium  and  coarse  sand)  and  two  types  of  glass  beads 
(small  and  large  beads).  The  granulometry  of  the  sand  and  beads 
is  given  in  Table  I.  In  the  case  of  medium  sand,  the  size  distribu¬ 
tion  of  the  particles  was  found  to  be  Gaussian,  and  the  mean  di¬ 
ameter  was  245  //m  (Table  I).  This  measurement  was  performed 
with  a  Mastersizer  2000.  The  results  are  given  in  Fig.  1. 

In  the  case  of  coarse  sand,  the  size  distribution  of  a  series  of 
244  particles  was  measured  manually  with  a  microscope.  Fig.  2 
is  a  photo  of  the  particles  analyzed:  the  ruler  at  the  top  of  the 
figure  gives  a  good  idea  of  the  particle  sizes.  As  can  be  seen  from 
the  histogram,  the  size  distribution  was  not  exactly  Gaussian, 
but  the  mean  diameter  of  the  particles  was  found  to  be  1 .5  mm. 

Reflection  coefficient  measurements  were  alternately  per¬ 
formed  on  sand  and  beads  of  about  the  same  mean  size  (Table  I) 
to  determine  whether  a  monosize  distribution  of  beads  might 
explain  the  main  features  obtained  with  a  Gaussian  distribution 
of  sand  particles.  The  results  of  this  comparison  should  be 
extremely  useful  for  future  modeling  purposes. 
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Fig.  1.  Histogram  of  grain  size  in  the  case  of  medium  sand  (particle  size  is  expressed  in  microns). [AU,  Pl€ClS€  J)  I'OV  1(1  C  tl€W  figUV€  With 

English  description] 


Histogram  of  sizes 


0  ^  N  -A  ^  N%  *  <v>  ^  ^  03 


Grain  diameter  (in  mm) 


Fig.  2.  Photograph  of  coarse  sand  particles  and  corresponding  size  distribution. 


Water  tank 


Fig.  3.  Experimental  arrangement. 


In  Table  I,  the  parameter  d(w)  is  defined  so  that  w  is  the  ratio 
between  the  volume  (weight)  of  the  grains  with  diameters  less 
than  d(w)  and  the  total  volume  of  all  the  grains  (dry  weight  of 
the  sediment  sample). 


Sand  and  glass  beads  were  placed  in  four  different  plastic 
boxes  of  the  same  size:  24  cm  x  18  cm  x  10  cm.  These  boxes 
were  then  immersed  in  a  water  tank  (2  m  long,  1.2  m  wide,  and 
1  m  deep);  see  Fig.  3.  The  water  had  been  in  the  tank  for  sev¬ 
eral  months  to  make  sure  that  it  contained  no  bubbles.  In  addi¬ 
tion,  careful  degassing  operations  were  performed  in  both  the 
sand  and  glass  bead  boxes  to  prevent  the  presence  of  bubbles 
in  the  medium.  After  these  degassing  operations,  both  the  sand 
and  bead  surfaces  were  carefully  flattened  to  remove  any  unde¬ 
sirable  ripples  on  the  water/sediment  and  water/bead  interfaces. 
Any  surface  roughness  could,  therefore,  be  said  to  be  due  solely 
to  the  granular  properties  of  the  medium. 

Four  Panametrics  transducers  were  used  to  cover  the  fre¬ 
quency  range  from  200  kHz  to  7  MHz.  These  transducers 
were  selected  so  that  their  spectra  overlapped  over  a  wide 
frequency  band.  This  frequency  overlapping  was  subsequently 
used  to  confirm  the  validity  of  the  measurements  by  comparing 
the  results  obtained  with  various  transducers.  The  physical 
characteristics  of  these  transducers  are  given  in  Table  II. 
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TABLE  II 

Physical  Characteristics  of  the  Panametrics  Transducers 


Transducers 

Bandwidth  at  -3dB 

Far-field  distance 

Beam  width 

Diameter  D 

Tl  (250  kHz) 

120  kHz 

6.5  cm 

7.5° 

40  mm 

T2  (500  kHz) 

250  kHz 

5.2  cm 

6° 

25  mm 

T3  (2.25  MHz) 

1  MHz 

7  cm 

2.5° 

13  mm 

T4  (5  MHz) 

2  MHz 

14  cm 

1.2° 

13  mm 

This  large  frequency  range  was  chosen  so  that  with  the  var¬ 
ious  sands  and  beads,  the  dimensionless  parameter  ka  would 
also  cover  a  very  large  range,  10-1  to  40,  where  k  is  the  wave 
number  and  a  denotes  the  radius  of  the  beads  and  the  mean  ra¬ 
dius  of  the  sand  particles. 

Before  the  beginning  of  the  experiments,  the  transducers  were 
calibrated  using  the  air/water  interface  as  a  perfect  reflector.  The 
distance  between  transducer  and  air/water  interface  was  about 
25  cm  in  all  the  calibration  experiments  except  for  that  in  which 
the  2.25-MHz  transducer  was  used,  where  for  the  sake  of  con¬ 
venience,  the  distance  was  set  at  only  16  cm.  It  is  important 
to  stress  that  these  distances  were  chosen  to  satisfy  transducer 
far-field  conditions  (Table  II). 

As  an  example.  Fig.  4  gives  the  time  signal  reflected  by  the 
air/water  interface  with  transducers  Tl,  T2,  T3,  and  T4.  This 
figure  also  gives  the  corresponding  spectra,  which  show  that  the 
frequency  bands  of  the  various  transducers  overlapped  consid¬ 
erably. 

A.  Experiments  on  Sand 

The  reflection  coefficient  of  a  sand  bottom  was  measured  at 
normal  incidence  and  at  the  same  distance  as  that  chosen  for 
the  calibration  at  the  air/water  interface.  In  this  case,  we  did 
not  have  to  apply  any  divergence  correction  to  the  signal  am¬ 
plitude.  The  backscattered  echoes  were  recorded  with  the  trans¬ 
ducers  in  various  horizontal  positions.  These  transducers  were 
mounted  on  a  carriage  moved  by  stepping  motors  (Fig.  3).  The 
increment  chosen  for  the  displacement  was  0.5  cm,  which  was 
much  smaller  than  the  footprint  of  the  beam  on  the  sand  bottom. 
The  footprint  size  of  the  beam  depends  on  the  size  of  the  trans¬ 
ducer,  the  working  frequency,  and  the  distance  from  the  trans¬ 
ducer  to  the  reflecting  surface.  In  our  experiments,  the  footprint 
diameter  was  about  2  cm  with  transducers  T3  and  T4  and  about 
7  cm  with  transducers  Tl  and  T2.  Examples  of  echograms  are 
shown  in  Fig.  5  in  the  case  of  medium  and  coarse  sand.  In  this 
experiment,  the  two  boxes  filled  with  medium  and  coarse  sand 
were  successively  ensonified  with  transducers  Tl,  T2,  T3,  and 
T4.  The  1 15  traces  shown  in  Fig.  5  correspond  to  a  total  trans¬ 
ducer  displacement  of  57  cm.  In  the  lower  part  of  the  figure,  the 
signal  reflected  by  the  medium  sand  is  shown,  and  in  the  upper 
part,  the  signal  reflected  by  the  coarse  sand  is  shown.  The  ef¬ 
fects  of  the  walls  of  the  boxes  (ripples  in  the  lower  part  of  the 
figure)  are  also  visible.  These  spurious  echoes  were  due  to  the 
diffraction  of  the  signal  by  the  plastic  edges  of  the  boxes.  This 
was  checked  by  making  simple  ray  travel  time  calculations.  This 
effect  was  particularly  pronounced  with  low-frequency  trans¬ 
ducers  because  of  their  large  aperture.  It  can  also  be  seen  that 
with  the  medium  sand  at  the  lowest  frequencies  (transducers  Tl 
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Fig.  5.  Received  echoes  from  medium  and  coarse  sand.  Transducers  T1  (250  kHz),  T2  (500  kHz),  T3  (2.25  MHz),  and  T4  (5  MHz). 


and  T2),  a  very  weak  echo  was  reflected  by  the  bottom  of  the 
plastic  box  (echoes  labeled  bottom  echo  in  Fig.  5).  At  relatively 
low  frequencies  (<500  kHz),  the  attenuation  of  the  signal  by 
the  sand  was  not  sufficiently  strong  to  absorb  the  energy  of  the 
transmitted  signal  over  a  distance  of  10  cm,  and  a  weak  echo  was 
reflected  by  the  plastic  structure  of  the  box.  All  these  undesirable 
echoes  were  eliminated  in  the  data  processing,  using  a  suitable 
time  windowing  procedure.  It  can  be  noted  that  these  spurious 
echoes  were  not  visible  in  the  case  of  the  coarse  sand  because  the 
signal  had  a  highly  scattered  structure  just  after  the  first  reflec¬ 
tion.  This  complex  structure  is  often  called  a  coda  in  the  seismic 
field.  In  our  experiments,  the  coda  was  due  to  the  multiple  scat¬ 
tering  by  the  grains.  This  process  also  began  to  appear  in  the 
case  of  medium  sand  at  very  high  frequencies  (Fig.  5). 

To  determine  the  reflection  loss,  we  measured  successively 
the  total,  coherent,  and  incoherent  intensities  defined  as  follows. 
If  pi(t )  is  a  time  series  of  the  signal  reflected  by  the  bottom 
when  the  transducer  is  in  a  position  denoted  by  i,  based  on  the 
definition  of  the  total  intensity  (It),  coherent  intensity  (ICoh)> 
and  incoherent  intensity  (7incoh),  we  can  write 

It  =  <I^(/)|2>1 
icoh  =  i  (Pi(m2 


-lincoh  —  J-t  icoh* 

To  obtain  the  reflection  loss  corresponding  to  the  total  intensity, 
the  following  procedure  was  applied. 

—  Take  the  Fourier  transform  of  the  time  signal  Pi(f)  = 

TF{pi(t)}. 

—  Take  the  square  root  of  the  total  averaged  intensity  7t(/). 

—  Compute  the  following  ratio:  square  root  of  the  total  av¬ 
eraged  intensity/absolute  value  of  the  Fourier  transform 
of  the  calibration  signal  at  the  air/water  interface.  The  re¬ 
sult  gives  the  reflection  coefficient  as  a  function  of  the  fre¬ 
quency.  The  reflection  loss  is  then  obtained  by  taking  the 
decimal  logarithm  of  this  expression. 

To  obtain  the  coherent  part  of  the  reflection  coefficient,  we 
used  the  coherent  intensity  instead  of  the  total  intensity.  The  in¬ 
coherent  part  of  the  reflection  coefficient  was  then  obtained  by 
calculating  the  difference  between  total  and  coherent  intensi¬ 
ties.  The  same  procedure  was  successively  applied  with  all  the 
transducers.  Both  total  and  coherent  intensities  were  obtained 
by  averaging  about  50  signals.  A  time  window  was  applied  to 
remove  the  signals  affected  by  spurious  echoes  before  the  aver¬ 
aging  process. 

The  reflection  loss  for  this  medium  sand  is  given  in  Fig.  6, 
where  also  the  contributions  of  the  incoherent,  coherent,  and 
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Fig.  6.  Reflection  losses  in  the  case  of  medium  sand  at  various  frequencies,  (a)  Transducer  Tl.  (b)  Transducer  T2.  (c)  Transducer  T3.  (d)  Transducer  T4.  (e)  All 
transducers. 


total  intensity  are  shown.  From  Fig.  6(a)  and  (b),  it  can  be  seen 
that  the  incoherent  intensity  level  was  very  low  (around  —35  dB) 
in  comparison  with  both  total  and  coherent  intensity  levels;  this 
means  that  the  reflection  was  mainly  coherent.  In  other  words, 
no  scattering  was  introduced  by  the  granular  structure  of  the  in¬ 
terface,  and  the  interface  acted  just  like  a  perfectly  flat,  homo¬ 
geneous  interface. 

Fig.  6(e)  shows  that  at  frequencies  less  than  1.4  MHz,  the  re¬ 
flected  level  was  constant  and  the  amplitude  was  around  —10 
dB,  which  is  in  agreement  with  classical  laboratory  measure¬ 
ments  on  sand  at  these  frequencies,  as  well  as  with  the  classical 
theory  of  reflection  on  a  water/sand  interface.  At  very  high  fre¬ 
quencies  (>3  MHz)  the  level  dropped  to  —30  dB,  which  corre¬ 


sponds  to  a  loss  of  amplitude  of  20  dB  between  the  lowest  fre¬ 
quency  (<1  MHz)  and  the  highest  frequency  (>3  MHz).  This 
was  certainly  due  to  a  complex  scattering  process  introduced 
by  the  sand  grains.  Note  that  in  this  region  of  the  spectrum,  the 
wavelength  began  to  be  of  the  same  order  of  magnitude  as  the 
mean  particle  size  (a/X  =  0.25  at  3  MHz).  Another  important 
feature  that  can  be  observed  here  is  the  transition  zone  between 
the  “low-frequency”  and  the  “high-frequency”  part  of  the  spec¬ 
trum.  In  this  transition  zone  [Fig.  6(c)],  the  reflection  loss  can 
be  seen  to  have  decreased  quasi-linearly  with  the  frequency. 

The  particular  point  at  which  a  break  occurred  in  the  slope 
corresponds  to  a  frequency  of  1.4  MHz  (wavelength:  1  mm), 
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Reflection  Loss  at  250  kHz  -  coarse  sand. 
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Fig.  7.  Reflection  losses  in  the  case  of  coarse  sand  at  various  frequencies,  (a)  Transducer  Tl.  (b)  Transducer  T2.  (c)  Transducer  T3.  (d)  Transducer  T4.  (e)  All 
transducers. 

which  corresponds  to  ka^ 0.76  (where  k  is  the  wave  number 
and  a  is  the  mean  particle  radius). 

Similar  measurements  were  performed  with  the  coarse  sand, 
the  granulometry  of  which  is  given  in  Table  I. 

The  reflection  loss  for  this  coarse  sand  is  given  in 
Fig.  7(a)-(d)  in  the  case  of  transducers  Tl,  T2,  T3,  and 
T4.  Fig.  7(e)  gives  the  overall  result  obtained  in  the  frequency 
range  from  200  kHz  to  6  MHz. 

It  can  be  observed  in  Fig.  7(a)  and  (b)  that  the  reflection  loss 
occurring  at  low  frequencies  (250  kHz  <  /  <  500  kHz)  and 
computed  with  total  intensity  was  around  —15  dB,  which  is  5 
dB  less  than  the  level  observed  with  medium  sand.  The  inco¬ 


herent  intensity  was  also  plotted  in  the  same  figure.  It  can  be 
seen  that  this  incoherent  intensity  was  relatively  high  in  compar¬ 
ison  with  the  coherent  part,  which  proves  that  some  scattering 
process  was  certainly  involved.  When  the  frequency  increases 
(/  >  3  MHz),  the  level  of  the  reflection  loss  seems  to  tend 
asymptotically  to  —30  dB.  Note  that  this  value  is  the  same  as 
that  obtained  with  medium  sand.  Moreover,  the  breaking  point 
seems  to  have  shifted  to  the  low-frequency  region.  It  is  difficult 
to  determine  the  exact  frequency  corresponding  to  this  point, 
but  200  kHz  seems  to  be  a  reasonably  accurate  estimate.  At 
this  breaking  point,  the  dimensionless  parameter  ka  was  about 
0.65,  which  is  very  similar  to  the  value  obtained  with  medium 
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Fig.  8.  Received  echoes  from  small  and  large  glass  beads.  Transducers  T1  (250 


sand.  This  finding  seems  to  indicate  that  the  main  parameter 
governing  the  reflection  loss  is  the  dimensionless  parameter  ka, 
which  depends  directly  on  the  ratio  between  the  particle  size  and 
the  wavelength. 

B.  Experiments  on  Glass  Beads 

Similar  experiments  to  those  performed  on  the  two  types  of 
sand  were  carried  out  with  glass  beads  of  two  sizes.  Fig.  8  gives 
the  echogram  obtained  with  transducers  Tl,  T2,  T3,  and  T4.  As 
in  the  experiments  with  sand,  the  lower  part  of  the  image  corre¬ 
sponds  to  the  small  beads  and  the  upper  part  to  the  large  beads. 
The  edge  effects  (ripples)  that  can  be  observed  here  in  the  lower 
part  of  the  image  were  again  due  to  the  scattering  by  the  edges  of 
the  plastic  box.  It  can  also  be  seen  in  this  case  that  the  water/sed¬ 
iment  interfaces  were  not  perfectly  horizontal,  which  resulted  in 
the  “bottom  trace”  being  tilted.  The  slope  of  the  bottom  was  cal¬ 
culated  and  found  to  be  around  1 .7°  with  small  beads  and  less 
than  —  1°  with  large  beads.  The  main  objective  in  preparing  the 
sample  was  to  obtain  a  perfectly  flattened  bottom  (with  no  rip¬ 
ples),  rather  than  obtaining  a  perfectly  horizontal  interface.  An 
error  of  +/— 1.5°  in  the  horizontality  was  taken  to  be  accept¬ 
able.  It  was,  therefore,  necessary  to  shift  the  signals  before  car¬ 
rying  out  any  data  processing  to  compensate  for  the  time  shift 
observed  between  signals  due  to  the  bottom  slope.  A  weak  echo 
generated  by  the  bottom  of  the  plastic  boxes  can  also  be  ob¬ 
served  at  low  frequencies  (transducers  Tl  and  T2).  This  weak 
echo  was  not  visible  with  the  large  beads  because  the  signals 
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were  strongly  scattered  just  after  the  reflection  by  the  interface 
(coda). 

The  reflection  loss  for  small  beads  (200  /rm)  is  given  in  Fig.  9. 
In  this  figure,  as  with  the  medium  sand,  we  can  see  that  the  re¬ 
flection  level  was  around  — 10  dB  and  that  this  amplitude  was 
practically  constant  over  the  whole  frequency  range  from  200 
kHz  to  1 .3  MHz.  The  incoherent  reflection  level  was  very  low  in 
comparison  with  the  coherent  one,  which  means  that  no  strong 
scattering  was  generated  by  the  interface.  It  can  also  be  observed 
that  the  reflection  levels  computed  with  the  total  intensity  and 
the  coherent  intensity  were  the  same.  This  means  that  the  re¬ 
flection  was  mainly  coherent.  As  the  frequency  increased,  the 
reflection  loss  computed  with  the  total  intensity  decreased  to  a 
value  of  —25  to  —30  dB,  as  previously  observed  with  medium 
sand. 

Fig.  9(e)  shows  the  reflection  loss  over  the  whole  frequency 
range  from  250  kHz  to  7  MHz  .  Between  300  kHz  and  1 .3  MHz, 
the  reflection  loss  remains  constant  before  showing  [Fig.  9(c)]  a 
linear  decrease  as  the  frequency  increases.  It  reaches  a  constant 
level  again  at  very  high  frequencies,  but  other  experiments  are 
now  required  at  higher  frequencies  to  confirm  this  pattern. 

The  reflection  loss  for  large  beads  showed  different  patterns, 
depending  on  the  frequency  of  the  ensonification.  Fig.  10  shows 
the  reflection  loss  obtained  with  transducers  Tl,  T2,  T3,  and  T4. 
The  reflection  behavior  observed  was  actually  more  chaotic  than 
with  small  beads,  but  the  mean  level  seems  to  have  been  stabi¬ 
lized  in  the  low-frequency  part  of  the  spectrum  (250  kHz  <  /  < 
500  kHz)  at  —15  dB,  which  is  5  dB  lower  than  with  the  small 
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Reflection  Loss  at  250  kHz  -  small  beads. 
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Fig.  9.  Reflection  loss  in  the  case  of  small  beads  at  various  frequencies,  (a)  Transducer  Tl.  (b)  Transducer  T2.  (c)  Transducer  T3.  (d)  Transducer  T4.  (e)  All 
transducers. 


beads.  In  the  high-frequency  regime  (1  MHz  <  /  <  3  MHz), 
the  reflection  loss  subsequently  decreased  with  the  frequency 
before  reaching  (at  /  >  3  MHz)  a  constant  level  (at  around  —30 
dB). 

In  this  case,  at  all  the  frequencies,  the  reflection  loss  was 
found  to  have  a  significantly  large  incoherent  component. 

If  we  compare  now  the  reflection  losses  obtained  with  large 
beads  and  coarse  sand,  we  can  see  some  similarities  (Fig.  11). 
This  observation  confirms  that  a  monosize  distribution  of  beads 
can  be  used  to  simulate  the  reflection  loss  in  a  complex  system 
of  sand  particles.  The  level  of  the  reflection  loss  at  very  high 


frequencies  also  seems  to  be  independent  of  the  particle  size  as 
well  as  being  independent  of  the  frequency. 

If  we  now  compare  the  reflection  losses  for  both  medium  sand 
and  small  beads  (Fig.  12),  we  can  see  that  the  two  patterns  were 
very  similar  in  the  low-frequency  range.  The  only  difference 
was  the  rate  of  decrease  of  the  reflection  loss  as  the  frequency  in¬ 
creased.  This  observation  again  confirms  that  small  beads  have 
a  very  similar  pattern  of  behavior  to  that  of  sand,  as  long  as  the 
bead  diameter  is  of  the  same  order  as  the  mean  size  of  the  sand 
particles.  Further  research  is  now  required  to  explain  why  a  dif¬ 
ference  in  the  amplitude  of  the  reflection  loss  was  observed  in 
the  transition  zone,  but  it  can  already  be  assumed  that  in  this  fre- 
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Fig.  10.  Reflection  loss  in  the  case  of  large  beads  at  various  frequencies,  (a)  Transducer  Tl.  (b)  Transducer  T2.  (c)  Transducer  T3.  (d)  Transducer  T4.  (e)  All 
transducers. 


quency  range,  where  large  scattering  effects  occur,  both  the  size 
and  the  shape  of  the  grains  have  greatly  affected  the  reflection 
loss. 

Last,  comparisons  between  small  and  large  beads  are  made  in 
Fig.  1 3  and  between  medium  and  coarse  sand  in  Fig.  14.  The  two 
curves  in  these  figures  both  point  to  the  same  conclusion.  With 
medium  sand  and  small  beads,  the  reflection  level  decreased 
conspicuously  in  the  high-frequency  domain;  it  can,  therefore, 
be  concluded  that  in  the  case  of  medium  sand  at  the  usual  sonar 
frequencies,  the  granular  structure  of  the  sediment  does  not  play 
any  role,  whereas  in  the  case  of  coarse  sand,  the  decrease  in 
the  reflected  level  begins  at  the  usual  sonar  frequencies.  When 
working  with  systems  operating  at  300  kHz  at  sea  on  coarse  sand 


bottom,  the  influence  of  the  granular  structure  of  the  sediment 
cannot,  therefore,  be  ignored. 

III.  Conclusion 

The  results  of  this  study  show  that  the  sound  reflected  by  a 
sand  bottom  depends  strongly  on  the  ratio  between  the  wave¬ 
length  and  the  particle  size.  This  effect  was  also  observed  with 
a  monosize  distribution  of  glass  particles.  The  main  result  ob¬ 
tained  in  this  paper  is  that  when  the  ka  parameter  is  less  than 
1,  the  reflection  loss  is  constant  and  the  sand  can  be  regarded 
as  a  continuous  homogeneous  medium.  When  the  ka  param¬ 
eter  is  around  1,  a  strong  decrease  in  this  reflection  loss  occurs, 
which  means  that  the  medium  behaves  like  a  scattering  medium. 
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Fig.  11.  Comparison  between  large  beads  and  coarse  sand. 
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Fig.  12.  Comparison  between  medium  sand  and  small  beads. 
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Fig.  13.  Comparison  between  small  and  large  beads. 

This  finding  shows  the  role  played  by  the  size  of  particles  in  the 
process  of  reflection  by  granular  media.  Although  we  can  ex¬ 
plain  at  this  stage  why  for  medium  sand  and  small  glass  beads 
the  reflected  level  is  frequency  independent  (this  level  depends 
mainly  on  the  sound  speed  and  the  density  contrast  between 
water  and  sediment  assumed  to  be  a  homogeneous  medium),  it  is 
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Fig.  14.  Comparison  between  medium  sand  and  coarse  sand. 

not  yet  possible  to  explain  why  the  level  dropped  to  around  —30 
dB  when  very  high  frequencies  were  used.  This  level  of  —30 
dB  was  obtained  with  both  medium  and  coarse  sand,  as  well  as 
with  both  small  and  large  glass  beads.  These  data,  therefore,  in¬ 
dicate  that  at  very  high  frequencies,  sand  behaves  mainly  like  a 
granular  medium  and  that  the  transition  between  continuous  to 
granular  behavior  is  linearly  dependent  on  the  frequency. 
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Abstract:  Sound  backscattering  from  water-saturated  granular  sediments 
at  frequencies  from  150  kHz  to  8  MHz  at  oblique  incidence  was  studied  in 
controlled  laboratory  conditions.  Two  kinds  of  sediments,  medium  and 
coarse  sands,  were  degassed,  and  their  surface  was  flattened.  In  these  condi¬ 
tions,  the  sediment  granular  structure  can  be  considered  as  a  controlling 
mechanism  of  backscattering.  Comparison  of  frequency  dependencies  of 
backscatter  for  the  two  sediments  with  different  mean  grain  size  shows  the 
existence  of  a  persistent  scaling  effect  that  allows  description  of  the  back- 
scattering  strength  as  a  function  of  one  parameter,  the  mean  grain  size/ 
wavelength  ratio. 
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1.  Introduction 

Reported  measurements  of  seabed  scattering  are  largely  confined  to  relatively  low  frequencies, 
up  to  300  kHz  (see,  e.g.,  Jackson  and  Richardson1  and  references  therein).  Still  higher  frequen¬ 
cies,  up  to  a  few  megahertz,  are  used  in  seafloor  imagery  (e.g.,  in  pencil-beam  sonars),  particu¬ 
larly  to  observe  dynamic  processes  at  the  seafloor.  From  the  physics  of  sound-sediment  inter¬ 
action  standpoint,  however,  interesting  effects  can  be  anticipated  at  these  high  frequencies  as 
the  wavelength  can  become  comparable  with  the  sediment  typical  grain  sizes.  Unfortunately, 
existing  observations  in  this  case  are  usually  noncalibrated  and,  therefore,  do  not  allow  mea¬ 
surements  of  system  independent  scattering  characteristics  of  the  seabed,  such  as  the  bottom 
scattering  strength.  A  recent  exception  is  presented  in  Greenlaw  et  al.~  where  the  backscatter¬ 
ing  strength  for  sand  sediments  in  shallow  water  was  measured  at  frequencies  from  265  kHz  to 
1 .85  MHz.  However,  physical  interpretation  of  results  obtained  in  this  work  is  quite  ambiguous 
due  to  a  significant  uncertainty  of  the  sediment  parameters.  As  the  authors  noted,2  the  sediment 
surface  was  only  relatively  smooth  and  had,  e.g.,  pockmarks  created  by  fish  foraging  on  the 
bottom  fauna.  A  small  but  detectable  sediment  gas  content  was  noticed  as  well.  Therefore, 
acoustic  observations  could  be  affected  by  the  small-scale  roughness,  the  presence  of  gas 
bubbles  in  the  sediment,  and  other  mechanisms,  e.g.,  those  related  to  possible  near-surface 
stratification  of  the  sediment.3'4 

Unfortunately  such  ambiguities  are  quite  usual  for  experimental  studies  of  sediment 
acoustics  in  shallow  water  where  the  dynamical  complexity  and  unpredictability  of  environ¬ 
mental  conditions  can  be  so  great  that  even  very  extensive  time-  and  labor-consuming  environ¬ 
mental  measurements  (such  as  those  at  the  recent  major  sediment  acoustics  experiment, 
SAX045'6)  may  not  be  enough  to  sufficiently  reduce  the  uncertainty  in  interpretation  of  acoustic 
data.  In  this  connection,  conducting  experiments  in  well-controlled  laboratory  conditions  can 
become  a  valuable  (but  much  less  expensive)  supplement  to  the  experiments  at  sea.  An  impor- 
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tant  advantage  of  laboratory  studies  of  sediment  scattering  is  that  they  make  it  possible  to  sim¬ 
plify  the  problem  by  reduction  of  the  number  of  controlling  parameters  and  allow  observation 
of  the  effects  of  different  scattering  mechanisms  separately.  For  example,  the  sediment  can  be 
degassed  to  eliminate  effects  of  micro-bubbles,  and  its  surface  can  be  flattened  to  exclude  con¬ 
tributions  of  the  sediment  roughness  at  scales  greater  than  the  grain  size.  This  approach  was 
used  by  Nolle  et  al.1  in  their  study  of  acoustical  properties  of  sand  sediments  at  two  frequencies, 
500  kHz  and  1050  kHz.  Williams  et  al.&  studied  broadband  backscattering  from  several  types 
of  the  sediments  in  a  wider  frequency  range,  about  200  kHz  to  2.5  MHz.  Unfortunately,  the 
bottom  scattering  strength  was  not  obtained,  as  the  measurements  were  still  not  calibrated. 

In  this  paper  (Sec.  2),  we  describe  experiments  conducted  at  the  CNRS/LMA  water 
tank  facility  on  backscatter  from  sands  over  a  wide  frequency  range,  150  kHz  to  8  MHz.  Note 
that  existing  methods  for  measuring  the  bottom  scattering  strength  normally  assume  using  nar¬ 
row  band  signals  (see,  e.g.,  Jackson  and  Richardson,1  Ch.  G.2,  pp.  497-502).  To  be  applied  to 
broadband  transducers,  significant  modifications  are  required.  One  such  modification  is  consid¬ 
ered  in  Sec.  3  of  this  paper.  It  made  possible,  to  our  knowledge  for  the  first  time,  measurement 
and  analysis  of  the  frequency  dependence  of  the  sediment  backscattering  strength  over  such  a 
wide  frequency  range.  The  results  are  given  in  Sec.  4.  The  analysis  assumes  that,  for  the  given 
conditions,  only  one  mechanism  of  scattering  can  be  dominating,  which  is  due  to  the  sediment 
granular  structure,  and  that  the  controlling  sediment  parameter  is  the  mean  grain  size.  In  this 
case  a  scaling  effect  is  possible:  Given  grazing  angle,  the  backscattering  strength,  which  is 
system  independent  and  depends  only  on  frequency  and  sediment  properties,  can  be  presented 
as  a  unique  function  of  only  one  parameter,  the  mean  grain  size/wavelength  ratio.  This  effect  is 
demonstrated  and  the  scaling  function  is  presented  in  Sec.  4  as  well. 

2.  Experiments 

Two  kinds  of  water-saturated  sediments  with  different  grain  sizes,  moderately  well  sorted  me¬ 
dium  and  coarse  sands,  were  chosen  for  the  study.  Sediment  properties  such  as  the  density, 
porosity,  and  grain  size,  were  measured  by  nonacoustical  methods.  The  medium  sand  had  the 
mean  grain  diameter  0.245  mm,  sediment/water  density  ratio  1.98,  and  porosity  36.5%.  The 
coarse  sand  had  the  mean  grain  diameter  1.55  mm,  sediment/water  density  ratio  2.02,  and  po¬ 
rosity  33%.  Prior  to  the  experiments,  the  two  sediments  were  stored  in  containers  filled  with 
water  treated  with  chlorine  to  exclude  the  presence  of  living  organisms  which  might  generate 
bubbles.  The  sediments  were  transferred  in  smaller  containers  (without  exposure  to  air)  to  a 
large  water  tank  (chlorine  treated  as  well)  where,  in  addition,  they  were  sieved  and  agitated  to 
eliminate  remaining  bubbles.  After  such  preparations,  the  two  sediments  were  placed  in  differ¬ 
ent  plastic  rectangular  boxes  of  the  same  size  with  horizontal  dimensions  17  X  23  cm  and  9  cm 
in  the  vertical.  The  sediment  surface  was  carefully  flattened  by  scraping  even  with  the  box  edges 
to  eliminate  roughness  at  scales  larger  than  the  sediment  grains.  Therefore,  the  necessary  mea¬ 
sures  were  taken  to  ensure  that  only  the  sediment  granular  structure  be  considered  as  a  domi¬ 
nating  mechanism  (rather  than  large  scale  roughness  and/or  gas  bubbles)  controlling  total  scat¬ 
tering  in  and  from  the  sediment. 

Experiments  on  backscattering  from  the  sediments  were  conducted  using  a  typical 
monostatic  geometry  with  the  transducers  acting  as  both  source  and  receiver  (see,  e.g.,  Jackson 
and  Richardson,1  Fig.  G3).  Two  broadband  Parametrics  transducers  with  nominal  center  fre¬ 
quencies  500  kHz  and  5  MHz  were  used  which,  altogether,  allowed  covering  continuously  the 
wide  frequency  range  of  150  kHz  to  8  MHz.  The  transducers  are  circular  pistons  with  radii 
12.5  mm  (500  kHz)  and  6.5  mm  (5  MHz).  Their  directivities  are  documented  over  all  used 
frequencies  and  well  described  by  the  corresponding  theoretical  function  l)=2J](ka  sin  (1) / 
(, ka  sin  0) ,  where  a  is  the  piston  radius,  k=2ir/k  is  the  wave  number  in  water,  /i(...)  is  the 
Bessel  function  of  the  first  order,  and  ft  is  the  angle  between  the  direction  of  observation  and  the 
maximum  response  axis  (MRA)  of  the  transducer. 

The  position  of  the  transducers  was  controlled  by  a  system  allowing  their  automatic 
vertical  and  horizontal  translation  with  0. 1  mm  accuracy  and  their  rotation  with  0.1°  accuracy. 
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For  calibration  of  each  transducer,  first  the  time  series  of  the  echo  signal  reflected  from  the 
water-air  interface  at  normal  incidence  were  measured  and  Fourier  analyzed.  Then,  the  trans¬ 
ducer  was  rotated  toward  the  sediment  surface  and  set  in  a  position  with  a  fixed  direction.  Hori- 
zontal  translations  of  the  transducer  (parallel  to  the  sediment  surface)  were  used  to  develop  a 
statistically  uniform  ensemble  of  realizations  for  the  echo  signals.  To  provide  the  necessary 
statistics,  a  number  of  horizontal  positions  of  the  transducer  was  set  automatically  with  a  con¬ 
secutive  horizontal  shift  of  0.5  cm.  At  each  position,  time  series  of  the  echo  signal  scattered 
from  the  sediment  (at  a  fixed  incidence  angle)  were  measured  for  a  number  of  pings  (up  to  64) 
and  coherently  averaged  to  reduce  possible  effects  of  noise. 

Fourier  spectral  analysis  was  performed  on  the  echo  signals  with  a  flexible  time¬ 
windowing.  The  windowing  was  used,  first,  to  reduce  the  effects  of  reverberation  from  the  sedi¬ 
ment  box  walls  and,  second,  to  accommodate  corresponding  distances  and  the  scattering  areas 
of  the  sediment  surface  with  the  frequency  dependent  beam  pattern  of  the  transducer  (as  dis¬ 
cussed  below).  The  windowed  echo  frequency  spectra,  F(f),  were  normalized  by  the  spectra  of 
the  calibration  signal,  F0(f).  These  normalized  spectra  provided  a  data  set  for  the  statistical 
processing.  Particularly,  their  second  moments  (the  squared  magnitudes  averaged  over  realiza¬ 
tions),  were  used  to  estimate  the  normalized  scattered  intensity 

ism0(f)=(\F(f)\2)/\FQ(f)\\  (i) 

and  to  obtain  the  frequency  dependencies  of  the  sediment  backscattering  coefficient,  ms,  and  its 
decibel  equivalent,  the  scattering  strength,  10  log  ms,  as  given  in  the  next  section. 

3.  Extracting  scattering  strength 

The  intensity  of  calibration  and  scattered  signals  in  the  frequency  domain,  Ia  and  Is  [see  Eq. 
(1)],  can  be  calculated  from  the  corresponding  signal  time  series  taken  in  the  observation  time 
window.  Specifically,  the  scattering  intensity  is  defined  as  the  squared  magnitude  of  the  Fourier 
transform  of  the  scattered  signal  time  series  taken  in  the  time  window,  t  e  [r, ,  t2],  related  to 
time -delays  occurring  in  the  process  of  scattering  and  defined  by  corresponding  range-time 
inequality 


ch  <  ri  +  rs<  ct2, 


(2) 


with  rt  and  rs  being  the  distances  from  the  source  and  receiver  to  the  ensonified  area  (or,  more 
accurately,  the  observation  area)  and  c  being  the  sound  speed  in  water. 

Consider  the  case  of  a  monostatic  geometry,  where  ri=rs=r=ct/2.  The  incidence 
angle  (measured  from  the  vertical),  0,  is  defined  as  follows:  cos  0=tfl t,  t^tf=2H/c,  where  /y is 
the  arrival  time  in  the  vertical  direction  (fathometer  return)  and  FI  is  the  vertical  distance  from 
the  source  to  the  sediment  surface.  For  the  angle  window  0  e  \_6X ,  02]>  taking  into  account  that 
tan  9dQ=d  In  /,  we  have  02-0\~  ( h~t])ltw  tan  0W  with  fu,=  (f1  +  /2)/  2  and  0n  =  (02+  9\)I2  be¬ 
ing  centers  of  the  corresponding  time  and  angle  windows. 

The  scattering  cross  section  per  unit  area  per  unit  solid  angle  (see,  e.g..  Ref.  1,  p.  23) 
or,  for  brevity,  the  scattering  coefficient,  ms,  can  be  defined  from  the  equation 


2 

_/  'ji-A 

>  —  lO  A^i 

r 


>fns 


(3) 


where  Ia  is  the  intensity  at  the  distance  r0=ct0  in  the  direction  of  maximum  directivity  (on  the 
transducer  MRA)  with  t0  being  the  arrival  time  of  the  calibration  signal,  and  Aef  is  an  effective 
observed  scattering  surface  obtained  from  the  directivity  D,  defined  with  respect  to  the  field 
amplitude,  by  integration  over  the  observation  area,  defined  by  Eq.  (2),  as  follows: 

Aef=  D4dA.  (4) 

JA(t) 
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The  integral  (4)  can  be  written  using  spherical  variables  (0,  ip),  which  are  the  incidence 
and  azimuth  angles  of  the  incident  wave.  However,  for  directional  sources  not  vertically  ori¬ 
ented,  it  is  more  convenient  to  describe  their  directivity  and  analyze  the  integral  (4)  in  angular 
coordinates  measured  not  from  vertical  direction,  but  from  the  transducer  MRA,  i.e.,  from  its 
depression  incidence  angle,  90.  Consider  new  variables  (not  spherical),  6'  =  0-  90  and  c p' 
=  <p  sin  0,  which  represent  angular  coordinates  measured  in  two  mutually  orthogonal  planes  so 
that  their  intersection  is  the  transducer  MRA,  where  6'  =  0  and  ip'  =  ip=Q.  One  of  these  two 
planes  is  always  vertical  (the  incidence  plane).  The  second  is  an  auxiliary  one  perpendicular  to 
the  incidence  plane  and  referred  to  as  the  slant  azimuth  plane. 

Using  the  new  variables,  the  integral  (4)  can  be  written  in  the  form 

r2 

^ef=  (5) 

COS  0W 

where  Hef  is  an  effective  solid  angle  of  the  observation  area 

r  6 2  r  tt  sin  6 

«er=  \~do\  D\e-e0,ip')dip' .  (6) 

J  J  -it  sin  0 


Note  that  the  integral  (6)  at  high  enough  frequencies  (and  correspondingly  narrow  directivities) 
has  a  maximal  value  at  9W=  90,  where  the  observation  angle  window  matches  with  the  maxi¬ 
mum  of  the  transducer  directivity. 

Assume  that  the  transducer  has  directivity  with  an  axial  symmetry,  so  that  for  the 
directivity  factor  we  have  D4(9' ,  p')  =  Q({)/ ()„),  where  0  is  the  angle  between  the  direction  of 
observation  and  the  transducer  MRA,  and  i)0  is  a  parameter  defining  the  directivity  angular 
width.  Note  that  at  i)n  -C  9,  limits  in  the  <p' -integral  in  Eq.  (6)  can  be  replaced  by  infinities,  and 
integration  can  be  performed  using  an  approximation  (0'2+  ip'2) 1/2.  Then,  in  the  matching 
case,  9W=  0q,  we  obtain 

nef«  p(a)  =  f  dxj  Q(^x2  +  y2)dy.  (7) 

The  integral  in  Eq.  (7)  has  two  simple  asymptotic  approximations: 

r  oo 

P  ~  Mxa ,  1,  Mi  =  I  Q(s)ds, 


P  ~  M2,  a^$>  1,  M2  =  2ir\  sQ(s)ds.  (8) 

J  o 


A  uniform  approximation  (for  arbitrary  a)  can  be  sought  in  the  form  of  a  scaling  expression 

P(a)^M2Pv(aMl/M2),  Pv{x)  =  {  77  >0.  (9) 

At  any  value  of  77,  this  approximation  has  two  correct  asymptotes  [Eq.  (8)]  for  small  and  large 
a.  This  allows  use  of  77  as  a  free  parameter  that  can  be  chosen  so  as  to  provide  the  best  fit  of  the 
scaling  function  [Eq.  (9)]  to  the  integral  (7)  at  intermediate  a,  taking  into  account  a  specific 
transducer  directivity  function 

In  our  case  of  the  circular  piston  transducers,  we  can  assume  that  at  high  frequencies 
#o=(^a)_1  <C  1  and  perform  the  integration  in  Eqs.  (7)  and  (8)  numerically  with  Q ( .v ) 
=  (27|(,v)  / ,v ) 4 .  For  the  integrals  in  Eq.  (8)  we  obtain  Mx  ~  2.4  and  M2  ~  5.8  ~  M2.  The  result  for 
the  integration  in  Eq.  (7)  is  shown  in  Fig.  1  by  the  solid  curve.  The  uniform  scaling  approxima¬ 
tion  given  in  Eq.  (9)  is  shown  for  different  77  by  dashed  curves.  It  is  seen  that  the  best  agreement 
is  provided  by  77—  4.  Therefore,  the  set  of  Eqs.  (1),  (3),  (5),  (7),  and  (9)  provides  a  practical 
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Fig.  1.  (Color  online)  Choosing  the  parameter  97  for  a  scaling  power  law  approximation  (9):  numerical  integration  for 
P(a)/M2  (solid  curve)  compared  to  the  scaling  function  Py(aM ]  / M2)  for  different  97  (dashed  curves).  The  best  fit  is 
provided  by  7/= 4. 


algorithm  for  extracting  the  backscattering  coefficient  from  analysis  of  the  normalized  scattered 
intensity  in  the  frequency  domain  using  time-windowing  of  the  echo-signal. 

4.  Results 

In  Fig.  2,  the  frequency  dependence  of  the  backscattering  strength  is  presented  in  a  wide  fre¬ 
quency  range,  150  kHz  to  8  MHz,  for  the  two  sediments,  medium  and  coarse  sand,  at  an  inci¬ 
dence  angle  of  50°  (grazing  angle  of  40°).  It  demonstrates,  in  particular,  a  noticeable  difference 
in  the  scattering  strength  of  the  two  sediments  at  frequencies  below  2  MHz.  For  example,  this 


Fig.  2.  (Color  online)  The  backscattering  strength  at  the  incidence  angle  50°  obtained  for  medium  and  coarse  sand 
using  the  broad  band  transducers  with  500  kFlz  and  5  MFlz  nominal  center  frequencies.  The  center  frequency  of  the 
500  kHz  transducer  was  nominal  indeed,  as  the  transducer  had  actually  another,  although  weaker,  maximum  in  its 
frequency  spectrum,  at  1.5  MHz,  which  was  used  to  obtain  additional  scattering  data. 
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difference  at  frequencies  400  kHz  to  800  kHz  is  about  20  dB.  It  can  mean  practically  important 
strong  sensitivity  of  the  scattering  strength  at  these  frequencies  to  the  mean  grain  size  of  the 
sediment.  Note  also  that  at  frequencies  above  2.5  MHz  the  difference  is  small  and,  therefore, 
the  scattering  strength  may  have  too  little  sensitivity  to  the  grain  size. 

The  ratio  of  the  mean  grain  diameter  to  the  acoustic  wavelength  in  water,  d/X,  varied, 
in  our  frequency  range,  from  0.024  to  1.3  for  the  medium  sand,  and  from  0.15  to  8.3  for  the 
coarse  sand.  Therefore,  in  the  interval  0A5<d/X<  1.3,  we  have  two  data  sets  for  sediments 
with  different  mean  grain  size.  This  interval  is  especially  interesting  as  it  corresponds  to  a  “tran¬ 
sition  frequency  regime”,  where  the  acoustic  properties  of  granular  sediments  can  change  dra¬ 
matically.  This  change  is  due  to  transition  from  the  “low  frequency  regime”,  where  the  ratio  d/X 
is  very  small  and  the  continuum  media  assumptions  are  valid,  to  the  “very  high  frequency  re¬ 
gime”,  where  this  ratio  is  not  small  and  the  sediment  must  be  considered  as  an  essentially  dis¬ 
crete  granular  medium.  At  these  high  frequencies,  intrinsic  (bulk)  scattering  due  to  the  sediment 
granular  structure  becomes  important.  For  example,  the  ratio  d/X  becomes  about  0.5  at  fre¬ 
quencies  around  3  MHz  for  the  medium  sand,  and  at  500  kHz  for  the  coarse  sand.  Correspond¬ 
ing  shift  (for  the  sediments  with  different  grain  size)  related  to  the  “transition”  effects  and 
appearance  of  the  new  dominating  scattering  mechanism  can  be  anticipated.  In  Fig.  2,  such  a 
scaling  shift  is  clearly  seen. 

In  Fig.  3,  the  scaling  effect  is  demonstrated  by  plotting  scattering  data  (same  as  in  Fig. 
2)  versus  the  ratio  d/X.  It  is  seen  that  the  backscattering  strength  collapses  to  a  function  of  only 
one  parameter,  at  least  over  the  interval  QA5<d/X<  1.3,  where  both  data  sets  are  available. 
This  confirms  the  assumption  that,  for  the  given  conditions,  only  one  mechanism  of  scattering  is 
dominating,  which  is  the  intrinsic  (bulk)  scattering  due  to  the  sediment  granular  structure,  and 
that  the  controlling  sediment  parameter  is  the  mean  grain  size.  In  this  case,  the  scattering 
strength  must  be  a  unique  function  of  the  ratio  d/X.  This  explains,  e.g.,  that  the  two  different 
grain-size  cases  match  when  shifted  horizontally,  but  no  shift  in  the  vertical  is  needed  (see  Figs. 
2  and  3).  This  is  because  the  scattering  strength  is  unique,  system  independent,  and  controlled 
only  by  the  sediment  parameters.  This  result  is  practically  important,  as  it  provides  both  predic¬ 
tion  and  inversion  capabilities:  if  the  scattering  strength  is  measured  at  one  frequency  and  one 
grain  size,  it  can  be  predicted  at  others  or  inverted  with  respect  to  the  grain  sizes.  Therefore,  the 
scaling  effects  demonstrated  here  can  potentially  be  used  in  acoustic  sensing  of  marine  sedi- 


d  IX 


Fig.  3.  (Color  online)  The  same  data  as  in  Fig.  2,  but  plotted  vs  the  sediment  mean  grain  diameter  ratio  to  the 
wavelength  in  water,  d/X. 
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ments  to  remotely  estimate  their  mean  grain  size.  Further  analysis  of  these  effects  might  be  an 
interesting  and  promising  subject  of  future  theoretical  and  experimental  studies. 
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MODELING  OF  NARROW-BAND  NORMAL-MODE 
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A  model  of  reverberation  in  a  shallow-water  waveguide  caused  by  the  volume 
heterogeneity  and  rough  interfaces  is  proposed.  Normal  modes  are  used  as  a  common 
approximation  of  the  one-way  propagation.  Narrow-band  approximation  allows  defining 
a  scattering  volume  for  each  mode-to-mode  channel  of  scattering  at  given  time. 
Smoothing  of  inter-channel  interference  provides  an  estimate  for  the  intensity  of 
reverberation  as  an  incoherent  summation  over  different  channels.  Local  scattering 
matrix  elements  quantify  the  energy  exchange  between  modes.  A  simple  relationship  of 
the  temporal  dependence  of  the  reverberation  intensity  with  the  scattering  matrix, 
attenuation  and  group  velocities  of  the  normal  modes  is  obtained.  Possible  effects  of 
multiple  scattering  are  discussed  and,  in  particular,  the  effect  of  backscattering 
enhancement,  which  is  not  described  by  existing  models  of  reverberation. 


1  Introduction 

Acoustic  reverberation  in  the  underwater  waveguide  is  complicated  particularly  because 
of  appearing  different  multiple  paths  of  sound  propagation  and  related  multiple  channels 
of  local  scattering.  In  shallow  water  conditions,  there  are  additional  complications 
caused  by  difficulties  to  resolve  these  multiple  paths  and  channels.  This  is  due  to  natural 
limitations  on  vertical  apertures  of  acoustic  systems,  their  instability  in  extremely 
dynamic  shallow  water  conditions,  spatial  complexity  of  this  environment,  and  other 
reasons.  There  are  various  approaches,  in  modeling  of  shallow  water  reverberation,  to 
take  into  account  these  multiple  paths  and  complicated  conditions.  They  all  assume,  in 
various  ways,  a  summation  of  signals  arriving  through  the  different  paths.  Normal 
modes  and  rays  are  commonly  used  to  describe  these  signals  and  have  their  own 
advantages  at  different  frequency  ranges  and  in  different  environments.  An  approach 
based  on  incoherent  summation  of  intensities  results  in  various  models  of  reverberation 
in  terms  of  the  energy  flux  integral.  These  models  have  an  advantage  of  simplicity,  but 
fail  in  some  cases  where  the  various  coherence  effects  may  be  important. 

In  this  paper,  a  general  formulation  of  the  reverberation  problem  is  given,  based  on 
a  unified  approach  to  volume  and  roughness  scattering  [1,2],  in  Section  2.  A  narrow- 
band  normal-mode  reverberation  model  is  described  in  Section  3.  A  simple  relationship 
of  the  temporal  dependence  of  the  reverberation  intensity  with  the  scattering  matrix, 
attenuation  and  group  velocities  of  the  normal  modes  is  obtained.  In  Section  4,  possible 
coherence  and  multiple  scattering  effects  are  discussed.  An  approach  is  described  which 
generalizes  a  multiple-forward-single-backscatter  (MFSB)  approximation  [3].  The 
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approach,  in  particular,  allows  modeling  of  the  effect  of  backscattering  enhancement  in 
strongly  inhomogeneous  shallow  water  waveguide  (which  cannot  be  described  by 
energy  flux  theory). 

2  General  formulation 

A  general  formulation  of  the  problem  can  be  based  on  a  unified  approach  to  volume  and 
roughness  scattering  in  layered  media  described  in  [1,2].  Consider  an  inhomogeneous 
medium  with  parameters,  mass  density,  p{r )  ,  and  compressibility,  AT (r )  ,  fluctuating 
in  space  relative  to  their  unperturbed  values,  p  and  K  ,  which  can  be  given  functions  of 
r  and  called  also  the  reference  parameters,  or  parameters  of  a  reference  medium  [1]. 
Here  r  =  { R,z },  R  =  { x,y }  denote  three-  and  two-dimensional  position  vectors.  For 

given  frequency  CO  and  spectrum  U (co)  of  a  source  of  radiation  located  at  position  r' , 
the  spectral  component  of  acoustic  pressure  will  obey  the  Helmholtz  equation 


[L(r)  +  £(r)]P(co,  r,r')  =  U(co)pS(r,  r' ) 

(1) 

L  =  pV(p~W)  +  k2 

(2) 

e  =  pco2  (k  -  k)  +  pV^p^  -p”')*7] 

(3) 

where  L  is  the  unperturbed  Helmholtz  operator,  with  k  —  (o!c  and 

c  =  {Kpym 

being  wave  number  and  sound  speed,  and  £  is  the  medium  perturbation  operator. 
Assume  that  a  solution  for  the  Green  function  of  the  reference  medium,  G ,  which 
obeys  the  equation 


L(r)G(co,r,r')  =  p5{r,r')  (4) 

is  known.  Then  the  solution  of  the  Eq.  (1)  can  be  expressed  through  the  Green  function 
using  the  equivalent  to  Eq.  (1)  integral  equation  of  the  form 

P(oo,r',r")  =  U(co)G{(D,r',r")+  \s{r)G{,(o,r,r')P{co,r,r")d\  (5) 

with  the  first  term  being  the  field  in  the  reference  medium  and  the  second  one  being  the 
field  perturbation.  Then  the  first  iteration  gives  the  first-order  perturbation  field,  or 
reverberation,  as  follows 

p{co,r',r ")  =  -U(co)^£(r)G(co,r,r')G(co,r,r")d3r  (6) 
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It  is  important  to  emphasize  here  that  the  mentioned  above  reference  medium  can  be 
chosen  rather  arbitrarily,  e.g.,  as  the  same  medium  but  with  heterogeneity  and  roughness 
smoothed  over  scales  comparable  and  smaller  than  the  wavelength.  In  such  a  smoothed 
medium,  with  only  large  scale  perturbations  left,  scattering  occurs  only  in  near  forward 
directions  and  various  small-angle  approximations  can  be  used.  In  this  case,  the  Green 
function  can  be  defined  as  a  propagator  to  describe  a  one-way  propagation  process  at 
given  frequency.  It  can  include  all  multiple  forward  scattering  effects  and,  in  particular, 
the  field  fluctuations  in  the  direct  path  described  by  first  term  in  Eq.  (5).  Reverberation 
itself  appears  as  a  field  perturbation,  or  the  second  term  in  Eq.  (5),  due  to  difference 
between  real  and  reference  medium  represented  by  a  medium  perturbation  function,  or, 
more  generally,  the  perturbation  operator,  which  can  be  of  an  arbitrary  form  and 
describe  both  volume,  continuous  or  discrete  (target-like),  and  roughness  perturbations. 
In  the  case  of  rough  interfaces,  the  volume  integral  Eq.  (6)  can  be  transformed  to  an 
integral  over  the  reference  interface,  and  the  perturbation  function  describes  the 
deviation  of  the  real  interface  from  the  reference  (or  smoothed)  one,  see  details  in  [1], 
Effects  of  multiple  forward  scattering  on  reverberation  are  included,  therefore,  to  the 
Green  functions  in  Eq.  (6). 

In  the  case  of  a  waveguide,  the  sound  propagation  is  described  through  a  summation 
of  different  signals  (corresponding  to  different  rays  and  modes).  In  this  case,  the  Green 
function  can  be  presented  as  a  superposition  of  generalized  rays  or  modes  as  follows 

G(co,r,r')  =  YjSn^y  I  ty) exp(ia>Tn(r ,r'\  co))  (7) 

n 

where  gn  and  Tn  represent  mainly  geometrical  factors  in  the  amplitude  and  phase  of 

the  n- th  ray/mode,  but  can  also  be  slow  functions  of  frequency  CO  to  allow  the 
ray/mode  amplitude  and  phase  frequency  dispersion.  In  this  paper,  we  consider,  as  a 
simple  example,  the  reverberation  in  a  plane  waveguide  in  terms  of  normal  modes.  In 
this  case,  the  Green  function  can  be  presented  in  the  form 

G(co,r,r')=\R-R'\-V2  p(i^n(o))\R-R'\)  (8) 

n 

where  iff  n  represent  the  modes  with  corresponding  horizontal  wave-numbers, 
E,n(c6)  =  col  Cn{co),  and  phase  velocities,  Cn(co),  at  given  frequency.  Equation  (8) 
gives,  therefore,  a  particular  example  of  Eq.  (7)  with 

gn(r,r'\o))=\R-RTU2  i//n(co,z)y/n(6),z')  (9) 

Tn{r,r'\oo)=\R-R'\ICn(co)  (10) 
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3  Narrow-band  normal-mode  reverberation 

3. 1  Narrow  band 

For  the  radiation  signal,  we  assume 

U(t)  =  u(t,(p)exp(-ico0t)  (11) 

where  CO0  is  the  carrier  frequency,  U  is  the  signal  envelope,  U  —  0  for  t  <  0  and 
t  >  T ,  with  T  being  the  radiated  pulse  duration,  which  is  long  enough,  so  that 
T  »  1  /  C00 .  A  (p  -parameter  is  introduced  here  for  generality  to  allow  a  possibility  of 

angular  directivity  of  the  source.  For  the  radiated  signal  spectrum  and  its  frequency  band 
width,  we  have 


U(co)  =  u(a>-co0)  (12) 

|<w-ty0 1  «  Mr  «  a>Q  (13) 

The  narrow  band  assumption  allows  the  following  approximation  for  the  mode  wave 
number  within  the  radiation  frequency  band 

4  O)  *  6,  (®0  )  +  Cn  '  (®0  X®  -  ^0  )  +  •••  ( 1 4) 

with  cn  =  (d<%n  /  dco)  1  being  the  group  velocity  of  the  mode. 

3.2  Reverberation  field 
Using  Eqs.  (6)  and  (8),  one  obtains 

P(t,r',r ")  =  exp(-id)0t)^pJn(t,r\r")  (15) 

y'.« 

where  p jn  represents  the  envelope  of  the  reverberation  arrival  corresponding  to  the 
j  —  n  channel  of  scattering  in  the  waveguide,  i.e.,  scattering  from  the  j  -th  mode  to  the 
n  -th  one,  and  is  described  as  follows 

pjn=Ajn{z\z'fi\vjn{t,\R-R\\R-R'\(p)Jjn(R)d2R  (16) 

Ajn  iz'  >  2"  )  =  y/j  (®0  >  Z'  )Vn  (®0  »  Z"  )  ( 1 7) 
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vjn(t,x',x")  =  (x'x")  U2u(t-x'/Cj  -x" I cn,(p) exp(/^yx'+/^wx")  (18) 

=  P  =  K’P’Z  (19) 

p 

Here  A  jn  (z1,  z")  defines  the  excitation  coefficient  for  the  j  —  n  channel  of  scattering 
at  given  positions  of  the  source  and  receiver,  v  n  describes  the  sound  propagation  along 
this  channel  with  corresponding  time  delay  and  complex  phase  (including  attenuation) 
and  also  defines  the  scattering  volume  for  this  channel  at  given  time,  and  J Jn  defines 

the  local  scattering  due  to  different  types  of  perturbations,  spatial  fluctuations  of 
compressibility  and  density,  and  the  interface  roughness,  with  corresponding  indexes 
(3  =  K,p,£  ,  see  details  in  [1,2]. 


3.3  Bistatic  reverberation  intensity 

Using  Eqs.  (15)  and  (16),  the  ensemble-average  intensity  of  reverberation  in  a  general 
bistatic  case  can  be  obtained.  Smoothing  of  inter-channel  interference  oscillations 
provides  the  following  estimate  for  the  intensity  of  reverberation  as  an  incoherent 
summation  over  different  channels 

<  |  p(t,r',r ")  |2>  «  £<|  Pjn^,r\r")  |2>  (20) 


<|  Pjn  \2>=\AJn(z',z")\2  j|v>(/,|JR-JR'|,|/(-/?"|^)|2  M jn(R)d2 R  (21) 
M Jn(R)  =  J<  Jjn(R-a/2)J*jn (R  +  a/ 2)  >exp(//cy„  • a)d2a  (22) 

(23) 


D  £  R-R'  D  ,  R-R” 

Kin=*eZl,n  „n+Re£.- 


R-R' 


R-R"  I 


where  M  n  is  the  local  scattering  matrix  coefficient  describing  the  energy  exchange 
between  modes  for  each  channel  within  the  corresponding  scattering  volume,  and  Kjn 
is  the  respective  local  scattering  vector. 

3.4  Quasi-monostatic  reverberation 

In  the  quasi-monostatic  case,  where  r'  ~  r" ,  assuming  also  that  t  »  T  ,  one  obtains 
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<  I  P  |2>  =  AJn(z',z")  |2  Mjn(Rjn(t),<p)exp(-2ajnt)  (24) 

t  j,n 

E{cp)  =  J|  u(t,<p )  |2  dt  (25) 

Rjn  (t)  =  t/ Cjn ,  cjH  =  (cT1  +  c-1  )_1 ,  ajn  =  cjn  Im(£y  +  £, )  (26) 

where  E(tp)  is  radiation  directivity  of  the  source,  or  the  total  energy  of  the  radiated 
signal  in  given  direction,  c  Jn  and  a  •„  are  the  speed  of  the  signal  return  and 
attenuation  coefficient  for  the  j  —  n  reverberation  channel,  and  R  is  the 

7  jn 

corresponding  horizontal  range.  Assuming  that  the  parameters  of  shallow  water 
waveguide  are  known,  Eq.  (24)  can  be  used  for  inversion  of  the  reverberation  strength, 
given  by  the  scattering  matrix  M  -n ,  from  the  range-direction  measurements  of 
reverberation,  or  the  clutter  maps. 

4  Multiple  scattering  effects 

4. 1  Coherence  functions 

Recall  now  that  the  reference  medium  can  be  defined  as  the  same  medium  but  with 
smoothed  parameters  describing  all  heterogeneity  and  roughness  with  scales  larger  than 
the  wavelength.  The  Green  function  of  the  reference  medium  takes  into  account  all 
corresponding  multiple  scattering  effects  for  near-forward  directions  of  propagation  and 
can  be  calculated  using  various  small-angle  approximations.  For  other  directions  (not 
near-forward),  including  backward,  all  the  scattered  waves  are  described  by  the 
perturbation  field  Eq.  (6)  and  included  to  reverberation,  which  can  be  presented  as 
follows 


p(co,r',r")  =  U  ( co)Gh  (co,r\  r" )  (27) 

with 

Gb(r\r")  =  -  \s{r)G{r,r')G{r,r")d\  (28) 

In  Eq.  (28)  and  following  equations,  the  frequency  is  omitted.  For  the  intensity  of  the 
scattered  field,  one  obtains 

(|G,(r\r")|2)=  J5(r„r2)r4(r„r2  \r',r")d\d\  (29) 

B(r\>r2)  =  (e(rx)e*(r2i)  (30) 
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r4 (ri » r2  k,,r")  =  (G(r1,r')G(r1,r")G*(r2,r')G*(r2,r")^  (31) 

where  B(f\,r2)  is  the  spatial  correlation  function  of  the  medium  perturbations,  and 
F4  is  the  fourth  statistical  moment  (or  the  fourth-order  coherence  function)  of  the 
propagation  field.  In  particular,  for  calculation  of  the  averaged  backscattered  intensity, 
or  the  monostatic  reverberation,  at  r'=  r" ,  the  following  coherence  function  is  needed 

r4 (ri  > r2  \  r\r')  =  (G2(rx,r')G*2(r2,r')^  (32) 

4.2  Backscattering  enhancement 

Consider  now  possible  effects  of  multiple  scattering.  They  are  especially  pronounced  in 
the  case  of  a  strongly  heterogeneous  medium,  where  the  propagating  field  is  completely 
randomized,  so  that  the  mean  field  vanishes,  i.e.,<  G  >=  0.  One  can  assume  that  in 
this  case  the  two  paths  of  propagation,  that  from  the  source  to  the  scattering  point,  and 
that  from  the  scattering  point  to  the  receiver,  become  mutually  incoherent,  or 
uncorrelated,  and  therefore  the  coherence  function  Eq.  (32)  can  be  presented  in  the  form 

r4(n»r2  \r'’r")uncor  =  ^2  (>i ,  r2  \  r')T2  (r, ,  r2  |  r")  (33) 

r2(ri>G  k)  =  (g( r\ ,  r)G*  (r2 ,  rfj  (34) 


where  T2  is  the  second  order  coherence  function  of  the  propagation  field.  This  function 
plays  a  fundamental  role  in  description  of  the  energy  flux  transfer.  In  particular,  its 
spatial  spectrum  represents  the  angular  distribution  of  the  energy  flux.  Therefore,  the 
reverberation  models  based  on  the  energy  flux  theories  actually  assume  that  Eq.  (33) 
holds  in  Eq.  (29)  for  all  source-receiver  configurations,  including  monostatic,  where  Eq. 
(33)  becomes 


r4(n^2  \r',r')uncor  =(r2(r,,r2  |  r'))2  (35) 

It  is  known  however  that,  in  this  case  of  strong  fluctuations,  the  field  obeys  a 
Gaussian  random  distribution.  As  a  result,  from  Eq.  (31)  one  obtains 

r4(ri,r2  |r,,r')  =  (G2(r1,r'))(G*2(r2,r,))  +  2(r2(r1,r2  |r'))2  (36) 

Comparison  of  Eqs.  (36)  and  (35),  provides  a  simple  explanation  for  the  so-called  effect 
of  backscattering  enhancement,  or  effect  of  the  source  localization  in  a  strongly 
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heterogeneous  medium.  Indeed,  the  integration  in  Eq.  (29)  using  Eqs.  (35)  and  (36), 
results  in  inequality 


(l  G6(r',r")  |2)  >  (l  Gb{r\r") |2)  (37) 

\  I  \  /  uncor 

which  completes  the  demonstration  of  this  effect. 

Note  that  in  the  case  of  multi-paths  reverberation  in  shallow  water,  the  effect  of 
backscattering  enhancement  should  be  clearly  pronounced.  As  it  was  mentioned  above, 
models  of  incoherent  reverberation  based  on  the  energy  flux  theories  are  unable  to 
describe  this  effect.  Its  demonstration  however  is  possible  if  rays/modes  models  are 
used,  but  only  if  possible  coherence  of  different  channels  of  scattering  are  taken  into 
account.  This  means  that  the  summation  of  signals  arriving  through  different  rays  or 
modes  should  consider  their  phases.  This,  in  particular,  would  make  possible  an 
interesting  numerical  experiment.  In  the  case  of  a  homogeneous  waveguide,  this  would 
result  in  a  very  complicated  spatial  “fine  structure”  of  reverberation  comprised  of 
interferential  maximums  and  minimums,  which  would  completely  mask  the  effects  at 
near  monostatic  locations.  However,  if  heterogeneity  in  the  waveguide  are  taken  into 
account,  this  fine  structure  vanishes  except  for  the  case  where  the  receiver  becomes 
close  enough  to  the  source  and  where,  therefore,  the  effect  of  backscattering 
enhancement,  or  the  source  localization,  becomes  remarkable. 
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Environmental  data  obtained  at  SAX99  site,  including  the  sediment  particle  size-depth 
distribution  and  the  water-sediment  interface  roughness  spectra,  were  used  to  compare 
contributions  of  volume  and  roughness  components  of  the  seafloor  scattering.  It  is 
shown,  in  particular,  that  contribution  of  gravel  and  shell  inclusions  and  coarse  sand 
fraction  in  total  scattering  at  SAX99  site  can  be  dominating  (over  roughness)  at  very 
high  frequencies  (about  100  kHz  and  higher)  and  grazing  angles  above  critical  (about  30 
degrees)  while  roughness  is  likely  a  dominating  mechanism  of  bottom  scattering  at  lower 
frequencies  and  grazing  angles  below  critical. 


1  Introduction 

Models  for  seabed  scattering  based  on  realistic  assumptions  about  scattering  mechanisms 
in  the  sediment  provide  relationships  between  various  seabed  properties  and 
characteristics  of  the  scattered  field  and  are  required  for  solution  of  various  practical 
problems  such  as  prediction  of  bottom  reverberation  given  seabed  properties  or/and 
inversion  of  various  seabed  parameters  from  acoustic  scattering  data.  There  are  different 
mechanisms  of  seabed  scattering  which  are  due  to  different  types  of  seabed  medium 
irregularities:  continuous  volume  fluctuations  of  the  sediment  acoustic  parameters  and 
discrete  volume  inclusions  (rock,  shell  hash,  etc),  roughness  of  the  seabed  interfaces,  as 
well  as  volume-roughness  interactions.  Testing  various  models  and  mechanisms  of  high 
frequency  seabed  scattering  was  one  of  primary  goals  of  recent  major  experiments  in 
sediment  acoustics,  SAX99  and  SAX04  [1-4]. 

Analysis  and  model/data  comparisons  for  SAX99  show  that  scattering  from 
continuous  volume  fluctuations  and  their  contribution  to  total  scattering  at  SAX99  site  is 
much  lower  than  that  of  roughness  and,  therefore,  can  be  neglected  [3].  On  the  other 
hand,  extensive  measurements  of  roughness  spectra  allow  estimation  of  the  level  of 
roughness  scattering,  which  shows,  within  a  reasonable  uncertainty,  that  roughness 
scattering  models  provide  a  good  description  of  bottom  scattering  at  frequencies  below 
100  kHz  [3].  However,  at  higher  frequencies,  there  is  a  significant  discrepancy  between 
observed  level  of  scattering  and  predicted  by  roughness  scattering  theories  because  the 
level  of  roughness  spatial  spectrum  is  too  low  at  the  corresponding  spatial  frequencies. 

Preliminary  analysis  shows  that  SAX04  scattering  is  also  very  complicated  and  that 
there  is  a  significant  discrepancy,  even  in  a  more  wide  frequency  range  than  for  SAX99, 
between  backscattering  data  and  results  predicted  by  models  taking  into  account  only 
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roughness  scattering.  Therefore,  other  mechanisms  of  scattering  must  be  considered  and 
other  approaches  are  required  for  understanding  SAX99  and  SAX04  data. 

In  this  paper,  a  model  of  volume  scattering  in  the  sediment  is  presented,  generalizing 
results  of  previous  models  [5-9].  The  model,  in  particular,  considers  the  case  of  stratified 
discrete  inclusions  in  the  sediment.  Environmental  data  set  obtained  at  SAX99  site, 
including  the  sediment  particle  size-depth  distribution  and  the  water-sediment  interface 
roughness  spectra,  is  used  to  compare  contributions  of  volume  and  roughness 
components  of  the  seafloor  scattering. 

2  Seabed  roughness  scattering  model/data  comparison 

Extensive  measurements  of  roughness  spectra  and  other  sediment  parameters  at  SAX99 
allow,  using  various  models  of  roughness  scattering,  prediction  of  frequency-angular 
dependencies  of  seabed  roughness  scattering  strength  and  their  comparisons  with 
observed  bottom  scattering  [3].  Results  of  such  comparison  are  illustrated  in  Figure  1, 
where  various  symbols  show  frequency  dependence  of  bottom  backscattering  strength  at 
two  fixed  grazing  angles,  35  and  20  degrees,  representing  angles  above  (a)  and  below  (b) 
the  critical  grazing  angle,  which  is  about  30  degrees  for  sandy  sediments  considered  here. 
Solid  and  dash-dot  curves  show  roughness  scattering  prediction  and  bounds  of 
uncertainty  correspondingly.  The  comparison  shows,  within  a  reasonable  uncertainty, 
that  roughness  scattering  mechanism  can  be  dominating  at  frequencies  below  100  kHz. 
However,  at  higher  frequencies,  predicted  roughness  scattering  is  substantially  lower 
than  observed  level  of  scattering  because  the  level  of  roughness  spatial  spectrum  was 
found  to  be  too  low  at  the  corresponding  spatial  frequencies  at  SAX99  site  [3]. 


SAX99  Backscattering,  %=35°  SAX99  Backscattering,  %=20° 


(a)  (b) 

Figure  1 :  Roughness  scattering  model  comparison  with  data  for  frequency  dependence  of  S  AX99 
seabed  backscattering  strength  at  a  fixed  grazing  angle  above  (a)  and  below  (b)  critical. 

Preliminary  analysis  shows  that  SAX04  scattering  is  also  very  complicated  and  that 
there  is  a  significant  discrepancy,  even  in  a  more  wide  frequency  range  than  for  SAX99, 
between  backscattering  data  and  results  predicted  by  models  taking  into  account  only 
roughness  scattering.  Therefore,  other  mechanisms  of  scattering  must  be  considered  and 
other  approaches  are  required  for  understanding  SAX99  and  SAX04  data  as  well. 
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3  Seabed  volume  scattering 

Let  us  assume  that  the  seabed  scattering  is  due  to  volume  heterogeneity  of  the 
sediment  layer  with  the  average  acoustical  parameters,  sound  speed  and  density, 
independent  from  the  depth.  Between  water  and  this  heterogeneous  scattering  layer,  an 
arbitrarily  stratified  transition  layer  is  allowed.  The  volume  backscattering  coefficient  in 
the  sediment,  or  the  cross  section  per  unit  sediment  volume,  m  ,  can  be  treated  generally 

as  dependent  on  the  depth.  In  this  case,  the  seabed  backscattering  coefficient,  or 
scattering  cross  section  per  unit  area  of  seabed  surface,  m  ,  can  be  presented  in  the  form 

ms  =1 W  I4  M  2  \a  mv(z)cx  p(  -zt  hp)dz  ■  ^ 

Here,  W  is  the  sound  transmission  coefficient  of  the  water-sediment  transition  layer  (or 
interface),  p  —  plp  is  the  sediment/water  density  ratio;  h  is  the  depth  of  sound 

penetration  into  the  sediment  (see,  e.g.,  [5,6]). 

In  a  particular  case,  where  the  volume  backscattering  coefficient  of  the  sediment, 
m  ,  is  depth  independent,  we  have  a  well  known  result  [5-9] 

m  =1 W  I4  u~2m  h  (2) 

s  r  v  p 

For  the  transmission  coefficient,  in  a  particular  case  of  a  flat  water-sediment  interface,  we 
have  W  —  l+V ,  where  V  is  the  seabed  reflection  coefficient. 


4  Incoherent  discrete  scattering  model 


Volume  scattering  in  the  sediment  can  be  attributed  to  two  different  kinds  of 
heterogeneity,  continuous  and  discrete.  Here,  in  this  paper,  we  consider  the  case  of 
scattering  from  discrete  inclusions  in  an  effective  fluid  sediment  and  assume  incoherent 
summation  of  the  scattered  intensities  for  different  scatterers  [6,7],  Some  considerations 
of  possible  coherent  effects  can  be  found  in  [10]  and  omitted  here.  In  the  frame  of 
incoherent  scattering  model,  the  scattering  cross  section  per  unit  sediment  volume  can  be 
presented  as  follows  [6] 

m  =  —  f  F  (k ,  a)a  ~2\ff  (a)cla  ■ 

An  y  ’ 

where 


V  v(a) 


aAVa 

v  A  a 


(4) 


is  a  dimensionless  volume  size  distribution  function,  Av  Iv  is  the  relative  part  of 


volume  occupied  by  particles  of  the  size  (equivalent  radius,  a )  within  the  correspondent 
interval  A  a  ,  k  is  the  wave  number  in  the  sediment  and  F  is  a  dimensionless  individual 
scattering  function  related  to  the  individual  scattering  cross  section  a  as  follows 

cr  =  a2F(k,a )  (5) 


Equations  ( 1  )-(5)  provide  a  solution  for  the  seabed  backscattering  coefficient  given 
an  arbitrary  depth-size  distribution  of  arbitrary  discrete  inclusions  in  the  sediment.  In  the 
calculations  below,  a  simple  case  is  considered.  Assume  that  inclusions  are  spherical, 
homogeneous  and  have  the  same  material  properties.  Also,  here  we  ignore  possible  shear 
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effects  within  inclusions  and  surrounding  sediments.  The  exact  solution  for  the 
backscattering  cross  section  and  the  scattering  function  in  this  case  can  be  found,  e.g.,  in 
[11].  A  reasonable  approximation  represents  a  solution  with  smoothed  interferential 
oscillations  at  ka>  1  [6]  and  used  in  the  following  sections  for  numerical  calculations  of 
the  integral  in  Eq.  (3). 


5  Size  distribution  at  SAX99 

An  environmental  data  set  obtained  at  SAX99  site  (near  Walton  Beach,  Florida)  [1] 
provides  depth-size  distribution  histograms  for  a  wide  range  of  grain  and  inclusions 
sizes.  Using  the  histograms  (courtesy  of  K.  Briggs)  the  volume  size  distribution  function 
\j/  defined  by  (4)  was  analyzed  and  results  are  presented  in  Figure  2. 


Figure  2.  Volume  size  distribution  for  the  SAX99  sediment  at  different  depths  (in-situ  data).  The 
solid  line  shows  the  best  multi-power  law  fit  to  the  environmental  in-situ  data.  The  dashed  line 
shows  the  distribution  that  would  provide  the  best  fit  for  the  acoustic  scattering  data. 

Different  symbols  represent  average  within  different  layers  of  2  cm  thickness  from 
the  cores  containing  top  26  cm  of  the  sediment.  It  is  seen  that  shape  of  the  size 
distribution  is  quite  different  for  different  depths  and  different  size  intervals.  For 
example,  concentration  of  shell  fraction  is  low  near  the  surface  at  0-2  cm  depth,  while 
there  is  a  remarkable  shell  layer  at  20-22  cm  depth.  The  solid  line  shows  the  best  multi¬ 
power  law  fit  to  the  environmental  in-situ  data  [6],  The  dashed  line  shows  the 
distribution  that  will  be  used  in  model/data  comparison  below. 

6  Volume  scattering  model/data  comparison 

Frequency  dependence  of  the  seabed  backscattering  strength,  10 log n?  ,  obtained  at 
SAX99  [3]  is  shown  in  Figure  3  at  a  fixed  grazing  angle,  %  =  35° ,  by  various  symbols. 

The  curves  show  results  of  calculations  for  the  model  of  discrete  scattering  described 
above,  using  Eqs.  (l)-(5).  Acoustic  parameters  of  the  sediment  in  the  model  correspond 
to  environmental  measurements  at  the  SAX99  site  [1],  The  grain  size  distribution  was 
taken  as  shown  in  Figure  2  (solid  curve)  and  with  changed  location  of  the  shell  layer 
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(dashed  curve)  to  demonstrate  an  effect  of  the  size  distribution  stratification.  Model/data 
comparison  shows  that  prediction  of  volume  scattering  using  real  size  distribution  is  a 
few  dB  lower  than  observed  bottom  scattering  at  SAX99  at  frequencies  below  about  100 
kHz.  However,  this  difference  can  be  easily  compensated,  if  the  shell  layer  is  located 
near  the  sediment  surface  (dashed  curve). 


SAX99  Backscattering,  /=35° 


Figure  3.  Frequency  dependence  of  the  bottom  backscattering  strength  at  SAX99  at  a  fixed  grazing 
angle  (35  degrees).  The  solid  line  shows  model  prediction  using  in-situ  size  distribution  data  (see 
Figure  1).  Comparison  of  the  solid  and  dashed  lines  shows  effect  of  the  shell  layer  location.  The 
shell  layer  is  placed  at  the  20-22  cm  depth  (solid)  and  0-2  cm  depth  (dashed)  respectively. 


SAX99  Backscattering,  y_=35° 


Figure  4.  Effect  of  a  thin  (3.5  mm)  top  homogeneous  (transition)  layer  at  high  frequencies.  The 
frequency  dependence  of  the  bottom  backscattering  strength  at  SAX99  shows  a  remarkable  roll-off 
above  1000  kHz.  The  solid  and  dashed  lines  show  results  for  a  model  of  discrete  scattering 
respectively  with  and  without  the  transition  layer. 
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In  Figures  4,  another  example  of  model/data  comparison  is  demonstrated,  showing 
that  the  model  provides  both  qualitative  and  quantitative  explanation  of  some  interesting 
experimental  results  obtained  at  SAX99  site  for  higher  frequencies  range  265-1850  kHz, 
including  backscattering  data  roll-off  above  1000  kHz  (see  Figure  4),  which  has  not  been 
explained  adequately  [12].  The  model  takes  into  account  the  fact  that  concentration  of 
shell  fraction  is  low  near  the  surface  at  0-2  cm  depth  (see  Figure  2)  which  allows  to 
consider  a  thin  (3.5  mm)  top  homogeneous  (transition)  layer  placed  on  a  sediment  half¬ 
space  with  depth  independent  size  distribution,  which  here  is  taken  as  shown  by  dashed 
curve  in  Figure  2.  The  result  is  shown  in  Figure  4  by  solid  curve.  The  dashed  line 
corresponds  to  the  case  without  the  transition  layer  and  perfectly  agrees  with  data  below 
300  kHz  but  fails  at  higher  frequencies. 

7  Volume  vs  roughness  scattering  at  SAX99 

Thus,  the  results  presented  above  show  that  the  model  of  scattering  from  discrete 
inclusions  in  the  sediment  can  be  a  good  descriptor  of  seabed  scattering  at  high 
frequency  (above  100  kHz)  and  grazing  angles  above  critical,  at  least  at  SAX99. 
Moreover,  a  simple  summation  of  the  two  mechanisms,  roughness  and  discrete  volume 
scatterers,  provides  a  good  model/data  comparison  for  backscattering  for  all  frequencies 
used  in  SAX99  at  grazing  angles  above  critical.  This  result  is  illustrated  in  Figure  5. 


SAX99  Backscattering,  x=35° 


Figure  5:  Model-data  comparison  for  frequency  dependence  of  SAX99  seabed  backscattering 
strength  at  a  fixed  grazing  angle  above  critical. 

However,  there  is  still  a  serious  problem  in  understanding  and  modeling  of 
backscattering  for  SAX99  at  high  frequencies  (about  200  kHz  and  higher)  at  sub-critical 
grazing  angles.  This  is  demonstrated  in  Figure  6.  Analogously  to  the  case  of  grazing 
angles  above  critical,  a  simple  summation  of  volume  and  roughness  scattering  was 
applied  assuming  no  interaction  between  these  two  mechanisms.  This  means,  that 
roughness  scattering  is  being  calculated  for  homogeneous  sediment  (with  no  inclusions) 
and  discrete  scattering  is  being  considered  assuming  a  flat  sediment  surface  (no 
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roughness).  This  simplified  approach,  used  successfully  at  higher  grazing  angles,  as 
demonstrated  in  Figure  5,  fails  in  attempts  to  explain  high  frequency  data  at  sub-critical 
grazing  angles  (see  Figure  6). 


SAX99  Backscattering,  x=20° 


Figure  6:  Model-data  comparison  for  frequency  dependence  of  SAX99  seabed  backscattering 
strength  at  a  fixed  grazing  angle  below  critical. 

The  discrepancy  can  be  due  to  ignoring  volume-roughness  interaction.  In  the  case  of 
flat  water-sediment  interface,  at  grazing  angles  below  critical,  volume  scattering  is 
significantly  reduced  because  of  small  depth  of  sound  penetration  into  the  sediment.  In 
the  case  of  randomly  rough  interface,  there  are  always  facets  with  local  grazing  angles 
above  critical,  which  causes  an  enhancement  of  sound  penetration  and  consequent 
enhancement  of  volume  scattering  in  the  sediment  [13].  The  effect  can  be  very 
significant  considering  the  fact  that  a  slope  of  roughness  at  sub-cm  scales  at  SAX99— 
SAX04  site  is  large  and  can  be  close  to  both  angle  of  repose  and  critical  angle  (about  30 
degrees). 

8  Discussion 

It  is  shown,  in  particular,  that  contribution  of  gravel  and  shell  inclusions  and  coarse  sand 
fraction  in  total  scattering  at  SAX99  site  can  be  dominating  (over  roughness)  at  very  high 
frequencies  (about  100  kHz  and  higher)  and  grazing  angles  above  critical  (about  30 
degrees)  while  roughness  is  likely  a  dominating  mechanism  of  bottom  scattering  at  lower 
frequencies  and  grazing  angles  below  critical.  A  simple  summation  of  the  two 
mechanisms,  roughness  and  discrete  volume  scatterers,  provides  a  good  model/data 
comparison  for  backscattering  for  all  frequencies  used  in  SAX99  at  grazing  angles  above 
critical. 

However,  there  is  a  serious  problem  in  modeling  of  backscattering  for  SAX99  at 
high  frequencies  (about  200  kHz  and  higher)  at  sub-critical  grazing  angles  where  there  is 
a  significant  model/data  discrepancy  which  can  be  due  to  ignoring  volume-roughness 
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interaction.  The  effects  of  such  interactions  can  be  very  significant  and  may  require 
further  theoretical  considerations  and  other  approaches. 

One  such  approach,  a  unified  approach  to  volume  and  roughness  scattering  [14,15], 
can  be  used  to  describe  volume-roughness  interactions  in  the  sediments.  It  is  exactly 
consistent  with  the  small  perturbation  method,  but  it  is  not  restricted  by  the  smallness 
requirement  for  roughness  height  and  slope,  which  makes  it  very  appropriate  tool  for 
considering  effects  of  non-small-slope  roughness  at  SAX99-SAX04  conditions. 
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A  model  of  discrete  scattering  in  marine  sediments  is  developed  which  permits,  using 
measurements  of  material  parameters  of  inclusions  and  their  size  distribution,  to  provide 
prediction  of  bottom  reverberation.  An  environmental  data  set  obtained  at  SAX99  site  (near 
Walton  Beach,  Florida),  including  the  size  distribution  and  sediment  acoustic  parameters, 
was  used  to  calculate  volume  component  of  the  seafloor  scattering.  Model/data  comparison 
shows,  in  particular,  that  discrete  scattering  from  gravel  and  shell  inclusions  and  coarse 
sand  fraction  can  be  a  dominating  mechanism  of  high  frequency  reverberation  at  grazing 
angles  above  critical  (about  30  degrees  for  sands).  Possibilities  for  inversion  of  various 
sediment  parameters  from  backs  cattering  data  are  discussed. 


1.  INTRODUCTION 

Acoustic  scattering  from  marine  sediments  with  discrete  inclusions,  such  as  shell 
fragments,  gravel  and  rocks,  can  be  described  in  terms  of  the  individual  scattering  functions 
of  discrete  targets  and  statistical  distributions  of  their  parameters  (size,  shape,  orientation, 
material,  etc)  [1,2].  An  environmental  data  set  obtained  at  SAX99  site  (near  Walton  Beach, 
Florida)  [3]  provides  the  necessary  data  on  the  size  distribution  for  a  wide  range  of  grain  and 
inclusions  sizes,  as  well  as  acoustic  parameters  of  the  sediments,  the  density,  sound  speed  and 
attenuation,  and  parameters  of  bottom  roughness.  Acoustic  measurements  of  the  seabed 
backscattering  strength  made  at  SAX99  [4]  along  with  environmental  data  provide  a  unique 
possibility  for  testing  various  models  of  scattering  and  for  evaluating  relative  contributions  of 
various  mechanisms  of  scattering. 

In  this  paper,  the  “inclusion  scattering”  model  is  demonstrated  as  a  good  descriptor  of 
seabed  backscattering  measured  at  SAX99  experiment  over  a  wide  frequency  range  (30  kHz 
to  300  kHz)  at  grazing  angles  above  critical  (about  30  degrees  for  sands).  Importantly,  the 
model  is  quite  simple  and  based  solely  on  environmental  measurements  of  a  small  number  of 
understandable  and  measurable  seabed  parameters.  Model/data  comparisons  and  analysis 
show  that  the  model  can  be  used  for  solution  of  both  forward  and  inverse  problems,  i.e., 


prediction  of  bottom  reverberation  given  sediment  properties  are  known  and  inversion  of  the 
sediment  parameters  from  multi-frequency  measurements  of  bottom  backscattering. 


2.  SEABED  BACKSCATTERING  COEFFICIENT 

Let  us  assume  that  the  seabed  scattering  is  due  to  volume  heterogeneity  of  the  sediment 
half-space  with  the  average  acoustical  parameters,  sound  speed  and  density,  independent 
from  the  depth.  In  this  case,  the  seabed  backscattering  coefficient,  or  scattering  cross  section 
per  unit  area  of  seabed  surface,  ms ,  is  related  to  the  volume  backscattering  coefficient  of  the 

sediment,  or  the  cross  section  per  unit  sediment  volume,  mv ,  as  follows  [1,5-7] 

ms  =|  W  |4  ju  2mvhp .  (1) 

Here,  W  is  the  sound  transmission  coefficient  of  the  water- sediment  interface,  ju  =  pi  pw  is 
the  sediment/water  density  ratio;  hp  is  the  depth  of  sound  penetration  into  the  sediment  at 
given  frequency  /  and  grazing  angle  in  water  %w , 


K  = 


R eV”2  -  cos2  z* 


2  2 
4klmJn  -cos 


2/?»0 


(2) 


where  k  =  2rf  /  c ,  (3  =  2k 5  and  n  =  no  (1  +  iS)  are  the  wave  number,  attenuation  coefficient 
and  complex  refraction  index  of  the  sediment,  5  is  the  loss  parameter,  no  =  cwl  c  is  the 

water/sediment  sound  speed  ration  (refraction  coefficient),  Re  and  Im  denote  real  and 
imaginary  parts  of  a  complex  value.  For  the  transmission  coefficient,  in  a  particular  case  of  a 
flat  water- sediment  interface,  we  have  W  =  1  +  V ,  where  V  is  the  seabed  reflection 
coefficient.  More  generally,  the  transmission  coefficient  can  take  into  account  the  seabed 
roughness  and/or  presence  of  a  stratified  top  transition  layer  in  the  sediment  [8]. 


3.  DISCRETE  INCOHERENT  SCATTERING  MODEL 

Volume  scattering  in  the  sediment  can  be  attributed  to  two  different  kinds  of 
heterogeneity,  continuous  and  discrete.  Here,  in  this  paper,  we  consider  the  case  of  scattering 
from  discrete  inclusions  in  an  effective  fluid  sediment  and  assume  incoherent  summation  of 
the  scattered  intensities  for  different  scatterers  [1].  Some  considerations  of  possible  coherent 
effects  can  be  found  in  [2]  and  omitted  here.  In  the  frame  of  incoherent  scattering  model,  the 
scattering  cross  section  per  unit  volume  of  sediment  can  be  presented  as  follows 

=  N<J^>  =  J  <?(“)'/' N(a)da ,  (3) 

where  a  is  the  scattering  cross  section  of  an  individual  scatterer  (inclusion)  and  <  a  >  is  its 
average,  the  ratio  N IV  is  the  average  number  of  scatterers  per  unit  sediment  volume,  a  is  a 


parameter  (or  a  set  of  parameters)  characterizing  individual  scatterers,  such  as  size,  shape, 
material  properties  and  others, 


VN(a)  = 


VA  a 


(4) 


is  the  probability  distribution  function  showing  the  number  of  scatterers  per  unit  volume  per 
unit  interval  of  the  parameter  a .  Here,  we  assume  that  inclusions  have  the  same  material 
properties,  and  a  is  considered  as  the  randomly  distributed  equivalent  radius  of  scatterers 
(radius  of  a  sphere  having  the  same  volume  as  a  non-spherical  inclusion).  In  this  case, 
relative  measurements  of  volume  and  weight  of  grains  are  equivalent.  Therefore,  using  a 
standard  fractional  weight  analysis,  so-called  volume  size  distribution  histograms  can  be 
obtained: 


AV 


,AVa  =  ^-m3ANa,AV  =  (\-P)V, 
3 

a 


(5) 


where  AFis  the  part  of  sediment  volume  occupied  by  inclusions  themselves  (grains),  P  is 
the  sediment  porosity.  As  y/  is  proportional  to  tj  at  small  // ,  we  will  be  using  a  normalized 

volume  size  distribution  function 


4  4  aAV 

y/Xa)  =  (l-P)y/Jri  =  -my/N(a)  =  -^-, 

which  is  independent  from  the  size  interval  used. 

Consider  a  simple  model  for  individual  scatterer  as  a  fluid  sphere  (representing  inclusions) 
with  acoustic  parameters,  sound  speed  and  density,  ca  and  pa ,  surrounded  by  a  fluid 
medium  (representing  surrounding  sediment)  with  respective  parameters  c  and  p  .  The  exact 
solution  for  the  backscattering  cross  section  in  this  case  can  be  found,  e.g.,  in  [9]  and 
presented  in  the  form 

<j  =  a2F(ka),  (7) 


where  F  is  a  dimensionless  scattering  function.  We  will  use  a  smoothed  function  of  the  form 


R1 

F(ka )  =  —^-(ka)4 
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where 
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PC2  V  Pg-P  R_PqCg  ~  PC 
PaCl J  Pa+Pt2’  PaCa+Pc' 


Using  Eqs.  (3,4)  and  (6,7),  we  obtain 


m  = 


4  n 


J F(ka)a  2y/v(a)da. 


(8) 


Eqs.  (1)  and  (8)  provide  a  solution  for  the  seabed  backscattering  coefficient  given  an  arbitrary 
size  distribution  of  discrete  inclusions  in  the  sediment  and  acoustic  parameters  of  these 
inclusions  and  surrounding  sediments. 


4.  SIZE  DISTRIBUTION  AT  SAX99 

An  environmental  data  set  obtained  at  SAX99  site  (near  Walton  Beach,  Florida)  [1] 
provides  the  size  distribution  histograms  for  a  wide  range  of  grain  and  inclusions  sizes.  Using 
these  histograms  (courtesy  of  K.  Briggs)  the  volume  size  distribution  function  defined  by  (6) 
was  analyzed  and  results  are  presented  in  Figure  1.  Different  symbols  represent  average 
within  different  layers  of  2  cm  thickness  from  the  cores  containing  top  26  cm  of  the  sediment. 


Fig.l:  Volume  size  distribution  of  the  sediment  grains  and  inclusions  at  SAX99  site 
( medium  and  coarse  sand  and  shell  fractions).  Solid  and  dashed  curves  correspond  to  a 
multi-power  law  approximation  with  slightly  different  parameters. 

It  is  seen  that  shape  of  the  size  distribution  is  quite  different  for  different  size  intervals.  We 
use  a  theoretical  approximation  for  the  volume  size  distribution  function  of  the  form 


y/Y{a)  =  A 


(1  +  (ala2)v-r) 


(1  +  (a/ a1)r) 


(9) 


This  function  is  of  a  multi-power  law  type  having  different  power  exponents  for  three 
different  intervals: 

y/v{a)~A,  a  <  av 

y/ v(a)  «  A(a / axyv ,  ax<a<a2, 

ipv(a)  ~  A(a2  / a1)-v (a/ a2YY ,  a  >  a2. 

The  solid  curve  in  Figure  1  represents  the  multi-power  law  function  with  the  best  fit  to 
measured  size  distribution  and  was  calculated  using  Eq.  (9)  with  the  following  values  of 
parameters:  A  =  0.5,  y  =  0.7,  v  =  5.873,  al  =  0.35mm,  a2  =  0.8mm.  Corresponding  values 
for  the  dashed  curve  are:  A  =  0.5,  y  =  0.7,  v  =  5.173,  ax  =  0.4mm,  a2  =  0.9mm ,  and  provide 
the  best  multi-power  law  lit  to  measured  backscattering  data  shown  in  Figure  2. 


5.  BACKSCATTERING  MODEL/DATA  COMPARISON 

Frequency  dependence  of  the  seabed  backscattering  strength,  1 0  log(mv) ,  obtained  at 
SAX99  [2]  is  shown  in  Figure  2  at  a  fixed  grazing  angle,  %w  =  35° ,  by  various  symbols. 


Frequency,  kHz 


Fig.2:  Frequency  dependence  of  the  seabed  backscattering  strength:  comparison  of  the 
discrete  scattering  model  and  experimental  data  obtained  at  SAX99.  The  grazing  angle  is  35 
degrees.  Solid  and  dashed  curves  correspond  to  a  multi-power  law  for  the  scatterers  size 

distribution  same  as  in  Fig.l. 

The  curves  show  results  of  calculations  for  the  model  of  discrete  scattering  described 
above,  using  Eqs.  (1,8).  The  grain  size  distribution  was  taken  as  shown  in  Fig.l  by  solid  and 


dashed  curves.  Acoustic  parameters  correspond  to  environmental  measurements  at  the 
SAX99  site  [1]  and  are  taken  as  follows:  p  =  2g/cm3,  c  =  1.7km/s,  8  =  0.01  (sediments); 
pa  =  2.7g/cm3 ,ca  =  5.7km/s,  8  =  0  (individual  grains);  and  pw  =  lg/cm3,civ  =  1.5km/s 
(water). 

The  model/data  comparison  represented  in  Figs.  1,2  and  its  analysis  show  that  the 
“inclusion  scattering”  model  is  a  good  descriptor  of  seabed  backscattering  measured  at 
SAX99  experiment  over  a  wide  frequency  range  (30  kHz  to  300  kHz)  at  grazing  angles  above 
critical  (about  30  degrees).  The  size  distribution  of  discrete  inclusions,  such  as  shells  and 
coarse  sand  particles,  is  the  critical  sediment  characteristic  in  this  model  and  an  important 
factor  determining  level  of  seabed  scattering.  Also  importantly,  the  model  is  quite  simple  and 
based  solely  on  environmental  measurements  of  a  small  number  of  understandable  and 
measurable  parameters.  The  model  can  be  used  for  solution  of  both  forward  and  inverse 
problems,  i.e.,  prediction  of  bottom  reverberation  given  sediment  properties  are  known  and 
inversion  of  the  sediment  parameters  from  multi-frequency  measurements  of  bottom 
backscattering.  For  example,  if  acoustic  parameters  of  the  sediment,  the  density,  sound  speed 
and  attenuation  are  known,  remote  acoustic  measurements  on  frequency  dependence  of 
bottom  reverberation  can  be  used  for  data  inversions  with  respect  to  size  distribution  of  the 
sediment  coarse  fractions.  From  the  other  hand,  the  size  distribution  can  be  carefully 
analyzed  in  laboratory  conditions,  and,  using  multi-frequency  measurements  of  bottom 
backscattering,  other  parameters  of  the  sediment  can  be  evaluated. 
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Bistatic  acoustic  measurements  offer  significant  advantages  in  detection  and  remote 
classification  of  natural  and  man-made  objects.  They  are  non-intrusive,  which  is 
especially  important  for  fragile  archeological  objects  buried  in  the  sediment.  The  results 
of  such  measurements  are  usually  expressed  through  the  scattering  cross  section  as  a 
function  of  frequency  and  directions  of  incident  and  scattered  waves.  While 
conventional  bistatic  techniques  are  performed  in  the  far-field,  with  accompanying  loss 
in  resolution,  we  propose  a  new  technique  employing  measurements  of  the  near-bottom 
Green's  function  and  its  exact  relationship  to  the  scattering  amplitude  (T-matrix).  Results 
of  numerical  experiments  are  presented  to  illustrate  the  new  technique  and  its 
advantages. 


1  Introduction 

Bistatic  acoustic  measurements  offer  significant  advantages  in  the  remote 
characterization  of  the  seabed,  and  in  detection  and  remote  classification  of  natural  and 
man-made  objects  in  marine  sediments.  They  are  non-intrusive,  which  is  especially 
important  for  fragile  archeological  objects  buried  in  the  sediment.  For  remote  acoustic 
characterization  of  the  objects,  in  many  cases,  rather  high  resolution  is  needed  to  capture 
classification  clues.  But  resolution  is  compromised  in  conventional  techniques  which 
require  the  measurements  to  be  made  far  enough  from  the  object  so  that  the  object  is  in 
the  far-field  region  of  the  source  and  the  receiver  is  in  the  far-field  region  of  the  object. 

In  this  paper,  a  new  technique  is  proposed  employing  near-bottom  bistatic 
measurements.  Ordinarily,  bistatic  measurements  are  extremely  difficult,  expensive  and 
time  consuming,  requiring  separate  source  and  receiver  platforms  and  accurate  position 
and  attitude  information  [1  ].  The  method  proposed  here  employs  a  single,  highly  portable 
platform  to  measure  the  near-bottom  Green's  function.  Given  this  Green’s  function, 
several  acoustically  important  quantities  can  be  determined,  e.g.,  the  scattering  cross 
section  and  the  reflection  coefficient,  which  usually  are  of  the  primary  practical  interest. 
The  basis  for  all  these  relationships  is  the  Fourier  transform  connection  between  the  near¬ 
bottom  Green’s  function  and  the  so-called  T-matrix.  The  T-matrix  provides  a  complete 
description  of  scattering,  giving  the  amplitude  of  the  plane  waves  scattered  in  all 
directions  for  any  given  incident  direction. 
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The  near-bottom  Green’s  function  can  be  measured  using  a  near-bottom  array.  The 
essential  point  is  that  the  near-bottom  array  provides  precise  measurements  that  can  be 
translated  into  practically  useful  estimates  of  scattering  and  reflection,  overcoming  the 
difficulties  and  uncertainties  of  previous  techniques  that  employ  transducers  placed 
higher  in  the  water  column.  In  this  paper,  theory  and  computer  simulations  are  used  to 
account  for  the  finite  aperture  and  discrete  spatial  sampling  of  the  actual  device. 

The  numerical  simulations  illustrate  the  new  technique,  and  its  advantages  in 
comparison  to  conventional  approaches.  In  particular,  the  proposed  technique  can 
increase  the  accuracy  and  angular  coverage  of  bistatic  measurements,  while  providing  a 
new  tool  for  scientific  measurements  of  acoustic  bottom  interaction  based  on  principles 
of  near-field  acoustic  holography. 

Two  simple  examples  are  considered.  The  first  example  is  a  flat  water- sediment 
interface  for  which  the  near-bottom  Green’s  function  was  calculated  via  wave-number 
integration  [2],  Then,  using  discrete  samplings  of  the  exact  Green’s  function,  an 
inversion  with  respect  to  the  reflection  coefficient  is  performed  using  the  discrete  Fourier 
transform.  The  result  is  compared  with  the  exact  expression  for  the  reflection  coefficient. 
This  procedure  allows  determination  of  the  required  array  size  and  the  element  spacing 
and  aids  in  assessing  the  trade-off  between  measurement  time  and  accuracy. 

As  another  example,  scattering  from  a  target  is  considered  using  a  point  scatterer 
model.  The  main  purpose  of  using  this  simple  model  is  to  provide  a  consistent  exact 
calculation  for  the  Greens  function  and  scattering  amplitude.  Again,  using  discrete 
sampling  of  the  exact  Green’s  function,  an  inversion  with  respect  to  the  scattering 
amplitude  (T-matrix)  is  performed  using  the  discrete  Fourier  transform.  The  result  is 
compared  with  the  exact  expression  for  the  T-matrix. 

An  important  advantage  of  our  approach  is  that  the  expressions  used,  in  particular, 
for  determining  the  reflection  coefficient  from  the  measured  field,  unlike  traditional 
approaches,  are  exact  and  independent  of  any  bottom  model.  Note  that  traditional 
approaches  of  measuring  the  reflection  coefficient  are  based  on  asymptotic  estimates  of 
the  above  mentioned  wave-number  integral  which  are  dependent  on  the  chosen  bottom 
model  and  measurement  geometry  and  do  not  permit  measurements  with  controlled 
accuracy.  Our  approach,  in  addition,  will  permit  us  to  evaluate  the  accuracy  and  validity 
of  traditional  approaches  simply  by  positioning  virtual  transducers  at  larger  distance  from 
the  bottom  and  to  compare  traditional  estimates  with  exact  expressions  for  the  reflection 
coefficient. 

2  Green’s  function  and  Scattering  Amplitude  (T-matrix) 

The  Green’s  function  of  the  fluid  medium  with  a  wave  number  k  obeys  the  equation 


V2G(r,r')  +  k2G(r,r')  -  S(r,r') 


and  can  be  presented  in  the  form 


G(r,r')  =  G0(r-r')  +  g(r,r'), 
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where  Q „  is  the  free-space  Green’s  function  describing  the  field  radiated  by  a  point 
source  and  g  describes  the  medium  response  (scattered  field),  r  =  (R,z)  and 
r'=(R',z')  are  position  vectors  of  the  receiver  and  the  source.  In  the  3D  case,  the  free- 
space  Green’s  function  can  be  presented  in  the  form 


G0(r)  =  - 


eik\r\ 

4  k  I  r 


In  the  2D  case,  the  free  space  Green’s  function  is 


G0(r)  =  -l-H^(k\r\y 


where  Hfp  is  the  outgoing  zero-order  Hankel  function. 

Let  all  the  scatterers  be  placed  in  a  lower  half-space,  while  the  source  and  receivers 
are  placed  in  the  upper  half-space,  i.e.,  z  >  0,  z'>  0  •  Below  the  source,  at  z  <  z\  the 

free-space  Green’s  function  can  be  represented  as  a  superposition  of  the  incident 
(downward-going)  plane  waves  as  follows 


2/G0(r-r')  = 


(2  nT 


tJ 


zk'-(r-r') 


r/K' 


where  k'=  (K',-V/(^1))  is  the  wave  vector  of  an  incident  plane  wave, 

v(K)  =  y[k*  —  K~  ,  n  =  2  or  n  =  3  in  the  2D  or  3D  case  respectively.  The 

corresponding  scattered  component  of  the  Green’s  function  also  can  be  presented  in  the 
form  of  a  plane  wave  integral  as  follows 


2/g(r,r')  = 


1 


JJ 


zk-r-z'k'-r' 


S(  K,K') 


dK  dK' 

_ » 

v(K )  v(K') 


(1) 


where  k  =  (K,u(^C))  is  the  wave  vector  of  a  scattered  going  upward  plane  wave, 
S(K,K')  is  the  so  called  scattering  amplitude  or  T-matrix  of  the  scattering  medium. 

The  T-matrix  provides  a  full  description  of  the  scattering  from  the  half-space  z  <  0 
regardless  of  its  complexity.  When  the  T-matrix  has  been  found,  all  the  basic 
characteristics  of  the  scattered  field  can  be  determined.  For  example,  the  reflection 
coefficient,  V (K) ,  for  an  arbitrary  plane-layered  medium  is  defined  by  the  following 

equation 


S(  K,K')  =  £(K-K>(iOR(K). 
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The  scattering  cross  section,  cr(K,K')>  for  an  arbitrary  3D  target  placed  in  the  lower 
half-space  can  be  determined  as  follows 

ct(K,K')  =  4tT  |  S(K,K')  |2. 

In  the  2D  case,  the  corresponding  equation  is 

o{K,K')  =  In  |  S(K,K')  |2  Ik. 


Equation  (1)  allows  determination  of  the  T-matrix,  if  the  Green’s  function  is  found, 
using  an  inverse  Fourier  transform 

S(  K,K')=  2i  v(K)v(f ]  JJ  e-lkr+lkW  g(r,  r'  )dRdR'  •  (2) 

{2n)n 

Equations  (1,2)  are  the  main  theoretical  results  of  this  paper. 

In  the  case  of  a  plane  layered  medium,  the  scattered  (reflected)  component  of  the 
Green’s  function  can  be  represented  as  follows 

g(r,r')  =  go(R-R';z  +  z') 


and  expressed  through  the  reflection  coefficient 


2ig0(R,z) 


1 

(2x)-' 


Je*KR+iV(^)zF(K) 


cK 

v{K) 


(3) 


The  reflection  coefficient  for  such  a  medium  can  be  obtained  as  follows 


F(K)  -  2  i  v(K)e~iv{K)z  j  e“/KRg0  (R,  z)dR  ■  (4) 

Equations  (3,  4)  were  used  in  a  numerical  example  of  inversion  of  the  reflection 
coefficient  below. 


3  Numerical  examples 

3. 1  Reflection  coefficient 

As  a  simple  numerical  example,  we  consider  a  2D  problem  of  reflection  from  a  flat 
interface  between  two  fluids  representing  water  and  the  sediment.  The  near-bottom 
Green’s  function  was  calculated  via  wave-number  integration  using  Eq.  (3)  for  a 
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sediment  having  density  ratio  2  and  refraction  index  0.8+i0.01.  In  Figure  1,  the  Green’s 
function  and  its  direct-path  and  reflected  components  are  shown. 


Figure  1.  An  example  of  the  2D  Green’s  function  and  its  components:  incident  field  (direct  path), 
reflected  and  total  fields  for  a  planar  interface  between  water  and  a  fluid  having  density 
ratio  2  and  refraction  index  0.8+i0.01.  A  discrete  sampling  of  the  total  field  was  then 
used  for  inversion  of  the  reflection  coefficient. 


grazing  angle,  deg 

Figure  2.  An  example  for  inversion  of  the  reflection  coefficient. 


If  the  Green’s  function  is  measured  by  a  horizontal  array.,  a  discrete  sampling  of  the 
exact  Green’s  function  is  provided.  Importantly,  complex  values  of  the  Green  function 
are  required,  including  its  magnitude  and  phase,  which  makes  this  procedure  similar  to 
that  used  in  near-field  holography.  Then,  inversions  with  respect  to  the  plane-wave 
reflection  coefficient  can  be  made  using  the  discrete  Fourier  transform  and  Eq.(4).  The 
results  of  inversion  for  the  magnitude  of  the  reflection  coefficient  and  its  comparison 
with  the  exact  expression  are  shown  in  Figure  2. 
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3.2  T-matrix  for  a  target  using  a  point  scattering  model 

As  another  simple  example,  we  use  the  point  scatterer  model  to  define  scattering  from  a 
2D  target.  This  simple  model  provides  a  consistent  exact  calculation  for  the  Green’s 
function  and  scattering  amplitude.  This,  again,  yields  a  measured  Green’s  function  using 
discrete  sampling,  and  then  allows  inversions  with  respect  to  the  scattering  amplitude  (T- 
matrix)  using  the  discrete  Fourier  transform  and  Eq.  (4). 


Figure  3.  An  example  the  Green’s  functions  for  a  2D  target,  discrete  samplings  of  which, 
imitating  the  measurements,  were  used  for  inversion  of  the  scattering  amplitude. 


L=1 51,  Ax=0.2>_,  /|nr=9Q:  L=5>„  Ax=0.51,  7|„,=90 


(a)  (b) 


Figure  4.  An  example  for  inversion  of  the  bistatic  scattering  amplitude  using  a  discrete  sampling 
of  the  Green’s  function,  corresponding  to  different  array  sizes,  15  (a)  and  5  (b),  and 
element  spacing,  0.2  (a)  and  0.5  (b),  with  all  dimensions  given  in  wavelengths.  The  solid 
line  shows  the  exact  solution. 
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An  example  of  this  procedure  is  illustrated  in  Figures  3  and  4.  In  Figure  3,  an  exact 
2D  Green’s  functions  calculated  for  a  target  comprised  of  point  2D  scatterers  is  shown. 
Its  discrete  sampling  imitates  the  measurements  with  an  array  and  then  is  used  for 
inversion  with  respect  to  the  scattering  amplitude  or  T-matrix.  The  result  of  such  an 
inversion  and  its  comparison  with  the  exact  expression  for  the  T-matrix  is  shown  in 
Figure  4.  This  type  of  numerical  procedure  can  be  very  useful,  in  particular,  in 
determination  of  the  required  array  size  and  the  element  spacing  and  can  aid  in  assessing 
the  trade-off  between  measurement  time  and  accuracy. 

4  Discussion 

In  this  paper,  a  new  technique  for  classification  of  the  sediment  and  buried  targets  is 
proposed  employing  measurements  of  the  near-bottom  Green's  function  and  its 
relationship  to  the  scattering  amplitude  (T-matrix).  Numerical  experiments  show  the 
applicability  of  the  new  technique  and  its  advantages  in  comparison  to  conventional 
approaches.  In  particular,  the  proposed  technique  can  increase  the  accuracy  and  angular 
coverage  of  bistatic  measurements,  while  providing  a  new  tool  for  scientific 
measurements  of  acoustic  bottom  interaction  based  on  principles  of  near-field  acoustic 
holography. 


a: — ? — 

N  N  /  / 

N  v  / 

\  — - — 

(b) 

Figure  5.  Sketch  of  near-bottom  bistatic  scattering  measurements.  The  linear  transmit 
and  receive  arrays  are  moved  as  shown  by  the  arrows  to  provide  2-D  synthetic 
apertures  covering  a  region  that  is  tens  of  wavelengths  in  dimension.  The 
actual  apparatus  would  have  a  drive  mechanism,  electronics,  computer,  and 
structure  designed  to  reduce  self-scattering. 
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Laboratory  measurements  with  free -field  and  buried  targets  are  the  next  logical  step 
in  the  development  of  this  technique.  Figure  5  shows  a  conceptual  picture  of  the  near¬ 
bottom  bistatic  measurement  apparatus  (a)  and  a  side  view  of  near-bottom  bistatic 
measurement  of  a  buried  target  (b).  Use  of  synthetic  apertures  for  transmission  and 
reception  simplifies  the  apparatus  and  reduces  the  level  of  unwanted  scattering  by  the 
array  elements. 
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